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PREFACE 


During recont years there has been an ever increasing interest in modern 
algebra not only of students in mathematics but also of those in physics, 
chemistry, psychology, economics, and statistics. Ny Modern Higher Alge¬ 
bra was intended, of course, to serve primarily the first of these groups, and 
its rather widespread use has assured me of the propriety of both its con¬ 
tents and its abstract mode of presentation. This assurance has been con¬ 
firmed by its successful use as a text, the sole prerequisite being the subject 
matter of L. E. Dickson’s First Course in the Theory of Equations. However, 

I am fully aware of the serious gap in mode of thought between the intuitive 
treatment of algebraic theory of the First Course and the rigorous abstract 
treatment of the Modern Higher Algebra , as well as the pedagogical difficulty 
which is a consequence. 

The publication recently of more abstract presentations of the theory of 
equations gives evidence of attempts to diminish this gap. Another such at¬ 
tempt has resulted in a supposedly less abstract treatise on modern algebra 
which is about to appear as these pages are being written. However, I have 
the feeling that neither of these compromises is desirable and that it would 
be far better to make the transition from the intuitive to the abstract by the 
addition of a new course in algebra to the undergraduate curriculum in 
mathematics a curriculum which contains at most two courses in algebra 
and these only partly algebraic in content. 

This book is a text for such a course. In fact, its only prerequisite ma¬ 
terial is a knowledge of that part of the theory of equations given as a chap¬ 
ter of the ordinary text in college algebra as well as a reasonably complete 
knowledge of the theory of determinants. Thus, it would actually be pos¬ 
sible for a student with adequate mathematical maturity, whose only train¬ 
ing in algebra is a course in college algebra, to grasp the contents. I have used 
the text in manuscript form in a class composed of third- and fourth-year 
undergraduate and beginning graduate students, and they all seemed to find 
the material easy to understand. I trust that it will find such use elsewhere 
and that it will serve also to satisfy the great interest in the theory of matrices 
which has been shown me repeatedly by students of the social sciences. 

I wish to express my deep appreciation of the fine critical assistance of 
Dr. Sam Perlis during the course of publication of this book. 


University of Chicago 
September 9, 1940 

v 


A. A. Albert 
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CHAPTER I 


POLYNOMIALS 

1. Polynomials in x. There are certain simple algebraic concepts with 
which the reader is probably well acquainted but not perhaps in the termi¬ 
nology and form desirable for the study of algebraic theories. We shall thus 
begin our exposition with a* discussion of these concepts. 

We shall speak of the familiar operations of addition, subtraction, and 
multiplication as the integral operations. A positive integral power is then 
best regarded as the result of a finite repetition of the operation of multi¬ 
plication. 

A polynomial fix) in x is any expression obtained as the result of the ap¬ 
plication of a finite number of integral operations to x and constants. If 
gix) is a second such expression and it is possible to carry out the operations 
indicated in the given formal expressions for f(x) and g(x) so as to obtain 
two identical expressions, then we shall regard f(x) and g(x) as being the 
same polynomial. This concept is frequently indicated by saying that fix) 
and g(x) are identically equal and by writing fix) s gix). However, we 
shall usually say merely that fix) and gix) are equal polynomials and write 
fix) = gix). We shall designate by 0 the polynomial which is the constant 
zero and shall call this polynomial the zero polynomial . Thus, in a discus¬ 
sion of polynomials fix) — 0 will mean that fix) is the zero polynomial. 
No confusion will arise from this usage for it will always be clear from the 
context that, in the consideration of a conditional equation fix) = 0 where 
we seek a constant solution c such that /(c) = 0, the polynomial fix) is not 
the zero polynomial. We observe that the zero polynomial has the 
properties 

0 • g(x ) = 0 , 0 + g(x) = g(x) 

for every polynomial gix). 

Our definition of a polynomial includes the use of the familiar term con¬ 
stant. By this term we shall mean any complex number or function inde¬ 
pendent of x. Later on in our algebraic study we shall be much more ex¬ 
plicit about the meaning of this term. For the present, however, we shall 
merely make the unprecise assumption that our constants have the usual 
properties postulated in elementary algebra. In particular, we shall assume 
the properties that if a and b are constants such that ab = 0 then either 
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a or b is zero; and if a is a nonzero constant then a has a constant inverse 
a~ l such that acr 1 = 1. 

If fix) is the label we assign to a particular formal expression of a poly¬ 
nomial and we replace x wherever it occurs infix) by a constant c, we ob¬ 
tain a corresponding expression in c which is the constant we designate by 
/(c). Suppose now that g{x) is any different formal expression of a poly¬ 
nomial in x and that fix) = g(x) in the sense defined above. Then it is 
evident that/(c) = g(c). Thus, in particular, if h(x), q{x), r(x) are poly¬ 
nomials in x such that fix) = hix)qix) + r(x) then /(c) = hic)qic) + r(c) 
for any c. For example, we have fix) = x 3 — 2x 2 + 3x, hix) = x — 1, 
qix) = x 2 — x, r(x) = 2x, and are stating that for any c we have c 3 — 
2c 2 + 3c = (c — l)(c 2 - c) + 2c. 

If the indicated integral operations in any given expression of a poly¬ 
nomial fix) be carried out, we may express fix) as a sum of a finite num¬ 
ber of terms of the form ax k . Here k is a non-negative integer and a is a 
constant called the coefficient of x k . The terms with the same exponent k 
may be combined into a single term whose coefficient is the sum of all their 
coefficients, and we may then write 

(1) fix) = a 0 x n + air*" 1 + . . . + a n _ X x + a n . 

The constants a* are called the coefficients of fix) and may be zero, but 
unless fix) is the zero polynomial, we may always take a 0 ^ 0. The ex¬ 
pression (1) of fix) with a 0 0 is most important since, if gix) is a second 
polynomial and we wTite gix) in the corresponding form 

(2) gix) = b 0 x m + b x x m ~ l + . . . + b m 

with bo 5 * 0, then/(x) and g(x) are equal if and only if m = a, a, = hi for 
i = 0, . . . , n. In other words, we may say that the expression (1) of a 
polynomial is unique, that is, two polynomials are equal if and only if their 
expressions (1) are identical. 

The integer n of any expression (1) of fix) is calk'd the virtual degree of 
the expression (1). If a 0 0 we call n the degree * of fix). Thus, either any 
fix) has a positive integral degree, or fix) = a n is a constant and will be 
called a constant polynomial in x. If, then, a n ^ 0 we say that the con¬ 
stant polynomial fix) has degree zero. But if a„ = 0, so that/(x) is the 
zero polynomial, we shall assign to it the degree minus infinity. This will 

* Clearly any polynomial of degree n ( , may be written as an expression of the form (1) 
of virtual degree any integer n ^ n 0 . We may thus speak of any such n as a virtual de¬ 
gree of f(x). 
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be done so as to imply that certain simple theorems on polynomials shall 
hold without exception. 

The coefficient ao in (1) will be called the virtual leading coefficient of this 
expression of f(x) and will be called the leading coefficient of /(x) if and only 
if it is not zero. We shall call/(x) a monic polynomial if a 0 = 1. We then 
have the elementary results referred to above, whose almost trivial verifica¬ 
tion we leave to the reader. 

Lemma 1. The degree of a product of two polynomials f(x) and g(x) is the 
sum of the degrees of f(x) and g(x). The leading coefficient of f(x) • g(x) is 
the product of the leading coefficients of f(x) and g(x), and thus y if f(x) and 
g(x) are monic , so is f(x) • g(x). 

Lemma 2. A product of two nonzero polynomials is nonzero aiid is a con¬ 
stant if and only if both factors are constants. 

Lemma 3. Let f(x) be nonzero and such that f(x)g(x) = f(x)h(x). Then 
g(x) = h(x). 

Lemma 4. The degree of f(x) + g(x) is at most the larger of the two degrees 
of f(x) and g(x). 

EXERCISES* 

1. State the condition that the degree of f{x) -f g(x) be less than the degree of 
either/(x) or g(x). 

2. What can one say about the degree of/Or) + g(x) if f(x) and g(x) have posi¬ 
tive leading coefficients? 

3. What can one say about the degree of / 2 , of / 3 , of f k for/ = /Or) a polynomial, 
k a positive integer? 

4. State a result about the degree and leading coefficient of any polynomial 
s(x) = f \ +...+/? for t > 1, ft = ffx) a polynomial in x with real coefficients. 

5. Make a corresponding statement about g(x)s(x) where g(x) has odd degree 
and real coefficients, s(x) as in Ex. 4. 

6. State the relation between the term of least degree in f(x)g(x) and those of 
least degree in f(x) and g(x). 

7. State why it is true that if x is not a factor of f(x) or g(x) then x is not a fac¬ 
tor of f(x)g(x). 

8. Use Ex. 7 to prove that if k is a positive integer then x is a factor of [/Or)]* if 
and only if x is a factor of /Or). 

9. Let / and g be polynomials in x such that the following equations are satisfied 
(identically). Show, then, that both/ and g are zero. Hint: Verify first that other- 

* The early exercises in our sets should normally be taken up orally. The author’s 
choice of oral exercises will be indicated by the language employed. 
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wise both / and g are not zero. Express each equation in the form a(x) — b(x) and 
apply Ex. 3. In parts (c) and (d) complete the squares. 

a) P + = 0 c) P + 2xpg z + (* 2 - x)g* = 0 

b) P - xY = 0 d) p + 2 xjg -xg 2 = 0 

10. Use Ex. 8 to give another proof of (a), (6), and (d) of Ex. 9. Hint: Show that 
if / and g are nonzero polynomial solutions of these equations of least possible de¬ 
grees, then x divides / = x}\ as well as g = xgi. But then/i and gi are also solutions 
—a contradiction. 

11. Use Ex. 4 to show that i if, g } and h are polynomials in x with real coefficients 
satisfying the following equations (identically), then they are all zero: 

a) p — xg 2 = xh 2 

b) P -xg 2 + h 2 = 0 

c) P + g* + Cr + 2)ft* = 0 

12. Find solutions of the equations of Ex. 11 for polynomials /, g t h with complex 
coefficients and not all zero. 

2. The division algorithm. The result of the application of the process 
ordinarily called long division to polynomials is a theorem which we shall 
call the Division Algorithm for polynomials and shall state as 

Theorem 1. Let f(x) and g(x) be polynomials of respective degrees n and rn, 
g(x) A 0. Then there exist unique polynomials q(x) and r(x) such that r(x) 
has virtual degree m — 1, q(x) is either zero or has degree n — m, and 

(3) f(x) = q(x)g(x) + r(x) . 

For let/Or) and g(x) be defined respectively by (1) and (2) with b 0 A 0. 
Then, either n < m and we have (3) with q(x) = 0, r(x) = fix), ora 0 A 0, 
n > m. If c k is the virtual leading coefficient of a polynomial h(x) of virtual 
degree m + k > m, a virtual degree of h(x) — ba l CkX k g(x) is m + k — 1. 
Thus a virtual degree of f(x) — bo l a 0 x n ^ m g(x) is n — 1, and a finite repetition 
of this process yields a polynomial r(x) = f(x) — b' 0 l (a {) x n '™ + . . .)g(x) of 
virtual degree m — 1, and hence (3) for q(x) of degree n — m and leading 
coefficient ajb^ 1 A 0. If also f(x) = qo(x)g(x) + r 0 (x) for r Q (x) of virtual 
degree m — 1, then a virtual degree of s(x) = r 0 (x) — r(x) is m — 1. But 
Lemma 1 states that if t(x) = q{x) — q 0 (x) A 0 the degree of s(x) = 
t{x)g(x) is the sum of rn and the degree of t(x). This is impossible; and 
t(x) = 0, q(x) = q 0 (x), r(z) = r 0 (x). 

The Remainder Theorem of Algebra states that if we use the Division 
Algorithm to write 


/(x) = q(x)(x - c) + r(x) 
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so that g{x) = x — c has degree one and r = r(x) is necessarily a constant, 
then r = /(c). The obvious proof of this result is the use of the remark in 
the fifth paragraph of Section 1 to obtain /(c) = q(c)(c — c) + r, /(c) = r 
as desired. It is for this application that we made the remark. 

The Division Algorithm and Remainder Theorem imply the Factor Theorem 
—a result obtained and used frequently in the study of polynomial equa¬ 
tions. We shall leave the statements of that theorem, and the subsequent 
definitions and theorems on the roots and corresponding factorizations of 
polynomials* with real or complex coefficients, to the reader. 

* If f(x ) is a polynomial in x and c is a constant such that /(c) = 0 then we shall 
call c a root not only of the equation f{x) = 0 but also of the polynomial/(x). 

EXERCISES 

1. Show by formal differentiation that if c is a root of multiplicity m of /(x) = 
(x — c) m q(x) then c is a root of multiplicity ra — 1 of the derivative /'(x) of /(x). 
What then is a necessary and sufficient condition that /(x) have multiple roots? 

2. Let c be a root of a polynomial /(x) of degree n and ordinary integral coeffi¬ 
cients. Use the Division Algorithm to show that any polynomial h(c) with rational 
coefficients may be expressed in the form bo + b\C + . . . -f h„_ic n “ 1 for rational 
numbers b 0 , . . . , i. Hint: Write h(x) = q(x)f(x) + r(x) and replace x by c. 

3. Let f(x) = x 3 + 3x 2 + 4 in Ex. 2. Compute the corresponding for each of 
the polynomials 

a) c 6 + 10c 4 + 25c 2 c) c 8 — 2c 4 + c 2 

b) c 4 + 4c 3 + Oc 2 + 4c + 1 d) (2c 2 + 3)(c 3 + 3c) 

3. Polynomial divisibility. Let /(x) and g(x) ^ 0 be polynomials. Then 
by the statement that g(x) divides f(x) we mean that there exists a poly¬ 
nomial q(x) such that/(x) = q(x)g{x). Thus, g(x) j* 0 divides/(x) if and 
only if the polynomial r(x) of (3) is the zero polynomial, and we shall say 
in this case that f(x) has g(x) as a factor , g(x) is a factor of f(x). 

We shall call two nonzero polynomials/(x) and ^(x) associated polynomials 
if /(x) divides g(x) and g(x) divides /(x). Then /(x) = q(x)g(x), g(x) = 
h(x)f(x), so that/(x) = q(x)h(x)f(x). Applying Lemmas 3 and 2, we have 
q(x)h(x) = 1, q(x) and h(x) are nonzero constants. Thus f(x) andg(x)are 
associated if and only if each is a nonzero constant multiple of the other . 

It is clear that every nonzero polynomial is associated with a monic poly¬ 
nomial. Observe thus that the familiar process of dividing out the leading 
coefficient in a conditional equation /(x) = 0 is that used to replace this 
equation by the equation g(x) = 0 where <j(x) is the monic polynomial 
associated with /(x). 
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Two associated monic polynomials are equal. We see from this that if 
g(x) divides f(x) every polynomial associated with g(x) divides f(x) and 
that one possible way to distinguish a member of the set of all associates 
of g(x) is to assume the associate to be monic. We shall use this property 
ater when we discuss the existence of a unique greatest common divisor 
(abbreviated, g.c.d.) of polynomials in x. 

In our discussion of the g.c.d. of polynomials we shall obtain a property 
which may best be described in terms of the concept of rational function. 
It will thus be desirable to arrange our exposition so as to precede the study 
of greatest common divisors by a discussion of the elements of the theory of 
polynomials and rational functions of several variables, and we shall do so. 

EXERCISES 

1. Let / = f(x) be a polynomial in x and define m(f) = x m f(l/x) for every posi¬ 
tive integer m. Show that m(f) is a polynomial in x of virtual degree m if and only 
if m is a virtual degree of f(x). 

2. Show that w*(0) = 0, m[m(f)] = /. 

3. Define/ = 0 if/ = 0 , and/ = n(/) if/is any nonzero polynomial of degree n. 
Show that m(J) = x m ~ n j for every m> n and that, if / 5 ^ 0 , x is not a factor of /. 

4. Let g be a factor of /. Prove that g is a factor of m(f) for every m which is at 
least the degree of /. 

4. Polynomials in several variables. Some of our results on polynomials 
in x may be extended easily to polynomials in several variables. We define 
a polynomial/ = f(x 1 , . . . , x q ) in x\ 9 . . . , x q to be any expression obtained 
as the result of a finite number of integral operations on Xi, . . . , x q and 
constants. As in Section 1 we may express f(x 1 , . . . , x q ) as the sum of a 
finite number of terms of the form 

(4) ax\ l xfy . . . x . 

We call a the coefficient of the term (4) and define the virtual degree in 
X\, ... 9 x q of such a term to be k\ + . . . + k g , the virtual degree of a par¬ 
ticular expression of / as a sum of terms of the form (4) to be the largest of 
the virtual degrees of its terms (4). If two terms of / have the same set of 
exponents fci, . . . , k q , we may combine them by adding their coefficients 
and thus write/as the unique sum, that is, the sum with unique coefficients, 

/ = /(*».*«) = ^ a kl . . . . . . x k „Q . 

k j » 0, 1. . . . . n ; 


(5) 
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Here the coefficients a kx . . . kq are constants and n, is the degree of f(x i, 

. . . , x q ) considered as a polynomial in x, alone. Also / is the zero poly¬ 
nomial if and only if all its coefficients are zero. If / is a nonzero poly¬ 
nomial, then some a kl . . . kq 9 ^ 0, and the degree of / is defined to be the 
maximum sum fci + . . . + k q for a kl . . . kq ^ 0. As before we assign the 
degree minus infinity to the zero polynomial and have the property that 
nonzero constant polynomials have degree zero. Note now that a poly¬ 
nomial may have several different terms of the same degree and that con¬ 
sequently the usual definition of leading term and coefficient do not apply. 
However, some of the most important simple properties of polynomials in 
x hold also for polynomials in several x i} and we shall proceed to their 
derivation. 

We observe that a polynomial / in X\ f . . . , x q may be regarded as a 
polynomial (1) of degree n = n 7 in x = x q with its coefficients a 0 , . . . , a„ all 
polynomials in x\ y . . . , x q -i and a 0 not zero. If, similarly, g be given by 
(2) with b Q not zero, then a virtual degree in x q of fg is m + n , and a virtual 
leading coefficient of fg is a Q bo. If q = 2, then a 0 and b 0 are nonzero poly¬ 
nomials in xi and a 0 fe 0 ^ 0 by Lemma 2. Then we have proved that the 
product fg of two nonzero polynomials / and g in X\> x 2 is not zero. If we 
prove similarly that the product of two nonzero polynomials in x h . . . , x q -\ 
is not zero, we apply the proof above to obtain a 0 bo 0 and hence have 
proved that the product fg of two nonzero polynomials in X\ y . . . , x q is not 
zero. We have thus completed the proof of 

Theorem 2. The product of any tiro nonzero polynomials in Xi, . . . , Xq 
is not zero. 

We have the immediate consequence 

Theorem 3. Let f, g, h be polynomials in xi, . . . , x q and f be nonzero , 
fg = fh. Then g = h. 

To continue our discussion wo shall need to consider an important special 
type of polynomial. Thus we shall call f(x i, . . . , x q ) a homogeneous poly¬ 
nomial or a form in Xi, . . . , x q if all terms of (5) have the same degree 
A; = A'i + . . . + k q . Then, if / is given by (5) and we replace Xi in (5) by 
yxi , we see that each power product . . . x q Q is replaced by y*»+ • • • +*«;r{i 
. . . x q i and thus that the polynomial f{yx u . .. , yx q ) = y k f(x i,.. . , x q ) identi¬ 
cally in y y x i, . . . , x q if and only if f(x i, . . . , x q ) is a form of degree k 
in x q . 

The product of t wo forms / and g of respective degrees n and m in the 
same X\, . . . , x q is clearly a form of degree m + n and, by Theorem 2, is 
nonzero if and only if / and g are nonzero. We now use this result to obtain 
the second of the properties we desire. It is a generalization of Lemma 1. 
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Observe first that all the terms of the same degree in a nonzero poly¬ 
nomial (5) may be grouped together into a form of this degree and then 
we may express (5) uniquely as the sum 

(6) / = f(x 1 , . . . , X q ) = fo 4" fn ) 

where/ 0 is a nonzero form of the same degree n as the polynomial/and/,- is 
a form of degree n — i. If also 

(7) g = g(x i, . . . , x q ) = g 0 + . . . + g m , 

for forms g> of degree m — i and such that go 0, then clearly 

(8) fg = h 0 4" • • • + ^mfn , 

where the 6, are forms of degree m + n — i and h 0 = fogo. By Theorem 2 
ho 9 * 0. Thus if we call / 0 the leading form of /, we clearly have 

Theorem 4. Let f and g he 'polynomials in xi, . . . , Xq. Then the degree 
of fg is the sum of the degrees of f and g and the leading form of fg is the prod¬ 
uct of the leading forms of f and g. 

The result above is evidently fundamental for the study of polynomials 
in several variables—a study which we shall discuss only briefly in these 
pages. 


6. Rational functions. The integral operations together with the opera¬ 
tion of division by a nonzero quantity form a set of what are called the 
rational operations . A rational function of x\, . . . , x q is now defined to be 
any function obtained as the result of a finite number of rational operations 
on Zi, . . . , x q and constants. The postulates of elementary algebra were 
seen by the reader in his earliest algebraic study to imply that every rational 
function of zi, . . . , x q may be expressed as a quotient 


( 9 ) 


f = • t %q) 

J b(x 1 , . . . , Xq) ’ 


for polynomials a(xi, . . . , x q ) and b(x i, . . . , x q ) 0. The coefficients of 
a(x i, . . . , x q ) and b(x i, . . . , x q ) are then called coefficients of /. Let us ob¬ 
serve then that the set of all rational functions in x\, . .•. , x q with complex 
coefficients has a property which we describe by saying that the set is 
closed with respect to rational operations. By this we mean that every rational 
function of the elements in this set is in the set. This may be seen to be 
due to the definitions a/6 + c/d = (ad + bc)/bd } (a/6) • (c/d) = (ac)/(bd). 
Here 6 and d are necessarily not zero, and we may use Theorem 2 to obtain 
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bd 7 * 0. Observe, then, that the set of rational functions satisfies the prop¬ 
erties we assumed in Section 1 for our constants, that is, fg = 0 if and only 
if / = 0 or g = 0, while if / ^ 0 then f~ l exists such that f*f~ l = 1. 

6. A greatest common divisor process. The existence of a g.c.d. of two 
polynomials and the method of its computation are essential in the study 
of what are called Sturm’s functions and so are well known to the reader 
who has studied the Theory of Equations. We shall repeat this material here 
because of its importance for algebraic theories. 

We define the g.c.d. of polynomials f\ (x), . . . , f t (x) not all zero to be any 
monic polynomial d(x) which divides all the fi(x), and is such that if g(x) di¬ 
vides every fi(x) then g(x) divides d(x). If d 0 (x) is a second such polynomial, 
then d(x) and d 0 (x) divide each other, d(x) and d 0 (x) are associated monic 
polynomials and are equal. Hence, according to our definition, the g.c.d. 
of/i(x), . . . } ft(x) is a unique polynomial. 

If g(x) divides all the fi(x), then g{x) divides d(x ), and hence the degree 
of d(x) is at least that of g{x). Thus the g.c.d. d(x) is a common divisor 
of the fi(x) of largest possible degree and is clearly the unique monic com¬ 
mon divisor of this degree. 

If dj{x) is the g.c.d. of /i(x), . . . , fj(x) and d Q (x) is the g.c.d. of d,(x) and 
fi+ i(x), then d 0 (x) is the g.c.d. of /i(x), . . . ,/,+i(x). For every common 
divisor h(x) of fi(x) y . . . ,/,+ i(x) divides /i(x), . . . ,/,(x), and hence both 
dj{x) and fj+i{x), h{x) divides d 0 (x). Moreover, d 0 (x) divides /,+i(x) and 
the divisor d y (x) of/i(x), . . . ,/,(x), d 0 (x) divides /i(x), . . . ,/ ;+1 (x). 

The result above evidently reduces the problems of the existence and 
construction of a g.c.d. of any number of polynomials in x not all zero to 
the case of two nonzero polynomials. We shall now study this latter prob¬ 
lem and state the result we shall prove as 

Theorem 5. Let f(x) and g(x) be polynomials not both zero . Then there 
exist polynomials a(x) and b(x) such that 

(10) d(x) = a(x)f(x) + b(x)g(x) 

is a monic common divisor of f(x) and g(x). Moreover , d(x) is then the unique 
g.c.d. of f(x) and g(x). 

Forif/(x) = 0, then d(x) is associated with g(x),a(x) = 1 and 6(x) = 
is a solution of (10) if g(x) is given by (2). Hence, there is no loss of generality 
if we assume that both /(x) and g{x) are nonzero and that the degree of 
g(x) is not greater than the degree of /(x). For consistency of notation 
we put 

(11) 


K{x) = f{x) , /ii(x) = g(x) . 
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By Theorem 1 

(12) h 0 (x) = qi(x)hi(x) + h 2 (x) , 

where the degree of h 2 {x) is less than the degree of hi(x). If h 2 (x) ^ 0, we 
may apply Theorem 1 to obtain 

(13) fh(x) = q 2 (x)h 2 (x) + h z (x) , 

where the degree of h z (x) is less than that of h 2 (x). Thus our division process 
yields a sequence of equations of the form 

(14) hi- 2 (x) = g i _ 1 (x)/ii_i(x) + hi(x) , 

where if r?, is the degree of hiix) then rq > n 2 > . . . , while > 0 unless 
hiix) = 0. We conclude that our sequence must terminate with 

(15) h r — 2 (x) = q r — i(x)/z r —l(#) h r ix) 

and 

(16) hr(x) 7 ^ 0 , hr- i(x) = q r (x)h r (x) 
for r > 1. 

filiation (16) implies that (15) may be replaced by h r ~ 2 ix) = 
[q r -\ix)q r ix) + l]h r (x). Thus h r (x) divides both h r ~ 2 {x) and h r -\(x). If we 
assume that h r (x) divides hi(x ) and hi- i(x), then (14) implies that h r (x) 
divides hi- 2 ix). An evident proof by induction shows that h r (x) divides 
both fto(x) = f{x) and hi(x) = g(x). 

Equation (12) implies that h 2 (x) — a a (x)/(x) + b 2 ix)gix) with a 2 (x) = 1, 
b 2 (x) = —qi(x). Clearly also hi(x) = a x (x)f(x) + bi(x)g{x) with adx) = 0, 
b x (x) = 1. If, now, hi- 2 (x) - ai_ 2 (x)/(x) + bi- 2 ix)g{x) and hi^i(x) = 
ai-i(x)f(x) + bi-iix)g(x) then (14) implies that hiix) = [a,_ 2 (x) — 
7 i_i(x)a»_i(x)]/(x) + [bi- 2 (x) - qi-iix)bi-iix)]g(x) = a*(x)/(x) + biix)gix). 
Thus we obtain ft r (x) = a r (x)/(x) + 6 r (x)g(x). The polynomial Kix) is a 
common divisor of /(x) and gix) and is associated with a monic common 
divisor d(x) = cKix ). Then d(x) has the form (10) for a (x) = ca r (x),6(x) = 
cbrix). We have already shown that d(x) is unique. 

The process used above was first discovered by Euclid, who utilized it in 
his geometric formulation of the analogous result on the g.c.d. of integers. 
It is therefore usually called Euclid's process. We observe that it not only 
enables us to prove the existence of d(x) but gives us a finite process by 
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means of which dix) may be computed. Notice finally that d{x) is computed 
by a repetition of the Division Algorithm on /(x), g(x) and polynomials se¬ 
cured from fix) and g(x) as remainders in the application of the Division 
Algorithm. But this implies the result we state as 

Theorem 6‘. The 'polynomials a(x), b(x), and hence the greatest common 
divisor d(x) of Theorem 5 all have coefficients which are rational functions 
with rational number coefficients of the coefficients of f(x) and g(x). 

We thus have the 

Corollary. Let the coefficients of f(x) and g(x) be rational numbers. Then 
the coefficients of their g.c.d. are rational numbers. 

If the only common divisors of f(x) and g(x ) are constants, then d(x) = 1 
and we shall call/(x) and g(x) relatively prime polynomials . We shall also 
indicate this at times by saying that fix) is prime to g(x) and hence also 
that g(x) is prime to/(x). When/U) and g(x) are relatively prime, we use 
(10) to obtain polynomials a(x) and b(x) such that 

(17) a(x)f(x) + g(x)b(x) = 1 . 

It is interesting to observe that the polynomials a(x) and b(x) in (17) are 
not unique and that it is possible to define a certain unique pair and then 
determine all others in terms of this pair. To do this we first prove the 

Lemma 5. Let f(x), g(x), and h(x) be nonzero polynomials such that f(x) 
is prime to g(x) and divides g(x)h(x). Then f(x) divides h(x). 

For we may write g(x)h(x) = fix)q(x) and use (17) to obtain [a(x)f(x) + 
b(x)g(x)]h(x) = [a(x)/i(x) + b(x)g{x)]f(x) = h(x) as desired. 

We now obtain 

Theorem 7. Let f(x) of degree n and g(x) of degree in be relatively prime. 
Then there exist unique polynomials ao(x) of degree at most m — 1 and b 0 (x) 
of degree at most n — 1 such that a 0 (x)f(x) + b 0 (x)g(x) = 1. Every pair of 
polynomials a(x) and b(x) satisfying (17) has the form 

(18) a(x) = a 0 (x) + c(x)g(x) , b(x) = b 0 (x) - c(x)f(x) 
for a polynomial c(x). 

For, if a(x) is any solution of (17), we apply Theorem 1 to obtain the 
first equation of (18) with ao(x) the remainder on division of a(x) by g(x). 
Then a 0 (x) has degree at most m — 1, a(x)/(x) + b(x)g(x) = a 0 (x)/(x) + 

[ b(x) + c(x)f(x)]g(x) = 1. We define b 0 (x) = b{x) + c(x)/tx) and see that 
b 0 (x)g(x) = — a 0 (x)/(x) + 1 has degree at most m + n — 1. By Lemma 1 
the degree of 6 0 (x) is at most n — 1, a 0 (x)/(x) + b 0 (x)g(x) = 1 as desired. 
If now ai(x) has virtual degree m — 1, 6i(x) virtual degree n — 1 and 
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«i(x)/(x) + b x (x)g(x) = a 0 (x)/(x) + b 0 (x)g(x), then fix) clearly divides 
[fco(x) ~ bi(x)]g(x). By Lemma 5 the polynomial b 0 (x) — bi(x) of virtual 
degree n — 1 is divisible by fix) of degree n and must be zero. Hence, 
&o(x) = &i(x),so that ai(x)/(x) = ao(x)/(x), ai(x) = a 0 (x). This proves a 0 (x) 
and b 0 (x) unique. But the definition above of a 0 (x) as the remainder on 
division of a(x) by gix) shows that then (18) holds. 

There is also a result which is somewhat analogous to Theorem 7 for the 
case where/(x) and ^(x) are not relatively prime. We state it as 

Theorem 8. Let f(x) ^ 0 and g(x) 0 have respective degrees n and m. 
Then polynomials a(x) ^ 0 of degree at rnost m — 1 and b(x) 0 of degree 
at most n — 1 such that 

(19) a(x)f(x) + b(x)g(x) = 0 

exist if and only if f(x) and g(x) are not relatively prime. 

For if the g.c.d. of/(x) and g{x) is a nonconstant polynomial d(x), wo 
have f(x) = fi(x)d(x),g(x) = gi{x)d(x) f gi(x)f(x) + [~Mx)g{x)] = 0 whore 
gi(x) has degree less than m and /i(x) has degree less than n. Conversely, 
let (19) hold. If f(x) and g(x) are relatively prime, we have a 0 (x)/(x) + 
bo{x)g(x) = 1, a(x) = a 0 (x)a(x)/(x) + a(x)6 0 (x)^(x) = g(x)[a{x)b Q (x) - 
a 0 (x)6(x)]. But then ^(x) of degree m divides a(x) ^ 0 of degree at most 
m — 1 which is impossible. 


EXERCISES 

1. Extend Theorems 5, 6, and the corollary to a set of polynomials /i(x), . . . , 
/<(*)• 

2. Let /i(x), . . . ,f t (x) be all polynomials of the first degree. State their pos¬ 
sible g.c.d.^ and the conditions on the /»(x) for each such possible g.c.d. 

3. State the results corresponding to those above for polynomials of virtual 
degree two. 

4. Prove that the g.c.d. of f(x) and g(x) is the monic polynomial of least possible 
degree of the form (10). Hint: Show that if d(x) is this polynomial then f(x) — 
q(x)d(x) + r(x), r(x) has the form (10) as well as degree less than that of d(x) and 
so must be zero. 

5. A polynomial f(x) is called rationally irreducible if f(x) has rational coefficients 
and is not the product of two nonconstant polynomials with rational coefficients. 
What are the possible g.c.d.’s of a set of rationally irreducible/»(x) of Ex. 1? 

6. Let f(x) = 0 be rationally irreducible, gix) have rational coefficients. Show 
that/(x) either divides gix) or is prime to {/(x). Thus, fix) is prime to gix) if the 
degree of g(x) is less than that of fix). 
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7. Use Ex. 1 of Section 2 together with the results above to show that a rational¬ 
ly irreducible polynomial has no multiple roots. 

8. Find the g.c.d. of each of the following sets of polynomials as well as of all 
possible pairs of polynomials in each case: 

a) /i = 2x 5 - x 3 - 2x 2 - 6x + 4 
/ 2 = x 4 + x 3 - x 2 - 2x - 2 

b) fi = 3x 4 + 8x 2 - 3 
U = x 3 + 2x 2 + 3x + 6 

c) f, = x 4 - 2x 3 - 2x 2 - 2x - 3 
/ 2 = x 3 + fix 2 + llz + 6 
fz = x 4 - 8x 2 - 5x + 6 

d) f\ — x 3 + 4x 2 + x — 6 
/ 2 = x 2 - 3x + 2 
/ 3 = x 5 — 3x 3 + x 2 — 4x — 4 

9. Let /(x) be a rationally irreducible polynomial and c be a complex root of 
f(x) — 0. Show that, if g(x) is a polynomial in x with rational coefficients and g(c) ^ 
0, there then exists a polynomial h(z) of degree less than that of fix) and with 
rational coefficients such that g(c)h{c) = 1. 

10. Let /(x) be a rationally irreducible quadratic polynomial and c be a complex 
root of /(x) = 0. Show that every rational function of c with rational coefficients 
is uniquely expressible in the form a + be with a and b rational numbers. 

11. Let/i, ...,/«be polynomials in x of virtual degree n and/i ^ 0. Use Ex. 4 
of Section 3 to show that if d(x) is the g.c.d. of/i, . . . , f t then the g.c.d. of/i, . . . ,/* 
is (7. Thus, show that the g.c.d. of n(/0, .. ., n{f t ) has the form x k d for an integer 

k > 0. 


e) fi = x R + 2x 4 - x 3 - 5x 2 - 6x - 3 
U = x 3 + x 2 + 3x + 3 
fz — X 4 + x 3 — x — 1 
/) /i = X 3 + 2x 2 - 3x - G 
fi = x 2 + x — 2 
fz = x 4 + 2x 3 + x + 2 


7. Forms. A polynomial of degree n is frequently spoken of as an n-ic 
polynomial. The reader is already familiar with the terms linear , quadratic , 
cubic t quartic , and quintic polynomial in the respective cases n = 1 , 2, 3,4, 5. 

In a similar fashion a polynomial in Xi, . . . , x q is called a q-ary poly¬ 
nomial. As above, we specialize the terminology in the cases q = 1 , 2, 3, 
4, 5 to be unary, binary , ternary , quaternary , and quinary. 

The terminology just described is used much more frequently in connec¬ 
tion with theorems on forms than in the study of arbitrary polynomials. 
In particular, we shall find that our principal interest is in n-ary quadratic 
forms. 

There are certain special forms which are quadratic in a set of variables 
x. 1 , . . . , x m , ?/i, . . . , y n and which have special importance because they 
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are linear in both a*i, . . . , x m and y\ y . . . , y ny separately. We shall call 
such forms bilinear forms. They may be expressed as forms 

j -1..... n 

(20) f = ^ ZiCt,,?/,-, 

* = 1, . . . , m 

so that we may thus write 

n m 

(21) / = ^asjj, a, = 5> w , 

i — l f * l 

and see that / may be regarded as a linear form in y Xy . . . , y n whose coeffi¬ 
cients are linear forms in x Xy . . . , x m . 

A bilinear form / is called symmetric if it is unaltered by the interchange 
of correspondingly labeled members of its two sets of variables. This state¬ 
ment clearly has meaning only if m = n; and / is symmetric if and only if 
/ = Hxianyj = ^y^ijXj. But / = and hence/ is symmetric if and 

only if m = n y 

(22) an = an (i, j = 1, . . . , n) . 

A quadratic form / is evidently a sum of terms of the type a { x; as well as 
the type CijXiXj for i ^ j. We may write an = a if « i; - = a ;i = £**,-,■ for i ^ j 
and have + anXjX ly so that 

n 

(23) / = ^ -r.a./X, 

t, J = 1 

(fl;y = Oji] i y j = 1, . . . , ?i) . 

We compare this with (22) and conclude that a quadratic form may be re¬ 
garded as the result of replacing the variables y u . . . , y n in a symmetric 
bilinear form in x iy ... y x n and y iy ... y y n by x Xy ... y x ny respectively. 
Later we shall obtain a theory of equivalence of quadratic forms and shall 
use the result just derived to obtain a parallel theory of symmetric bilinear 
forms. 

A final type of form of considerable interest is the skew bilinear form. 
Here again m = n y and we call a bilinear form / skew if f = f(xi y ... , x n ; 
2 /i, ... , 7/ n ) = —f(yi, . . . , ynl X\ y ... y x n ). Thus skew bilinear forms arc 
forms of the type 


(24) 


f = y' WaVi, 
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where 

(25) dja = —an (i, j = 1, . . . , n) . 

It follows that an + an = 0, that is 

(26) an = 0 (i = 1, . . . , n) . 

Hence / is a sum of terms a tJ (xi?/ ; — Xjy t ) for i j, i — 1, . . . , n — 1, 
j = 2, . . . , n. It is also evident that if we replace the y, by corresponding 
Xj, then the new quadratic form/(xi, . . . , x n ; Xi, . . . , x n ) is the zero poly¬ 
nomial. It is important for the reader to observe thus that while (22) may 
be associated with both quadratic and symmetric bilinear forms we must 
associate (25) only with skew bilinear forms. 

ORAL EXERCISES 

1. Use the language above to describe the following forms: 

a) x 3 -f 3 xy 2 -f z 3 d) x\ 4* 2 z x y x 

b) x\ + y\ e) xiy 2 - x 2 yi 

c) 2x\iji + x 2 yi + xiy 2 

2. Express the following quadratic forms as sums of the kind given by (23): 

a) 2x\ + 3x,x 2 

b) x? - x\ + 2xjx 3 + 2x\ - 4x 2 .T3 

8. Linear forms. A linear form is expressible as a sum 

(27) f a^Xi 4" • • • -}- a n x n . 

We shall call (27) a linear cotnbination of x\, . . . , x n with coefficients 
a 1 , ... f a n . The concept of linear combination has already been used with¬ 
out the name in several instances. Thus any polynomial in x is a linear 
combination of a finite number of non-negative integral powers of x with 
constant coefficients, a polynomial in xi, . . . , x q is a linear combination 
of a finite number of power products xji . . . xj<? with constant coefficients, 
the g.c.d. of fix) and g(x) is a linear combination (10) of f(x) and g(x) with 
polynomials in x as coefficients. 

The form (27) with ai = a 2 = ... = a n = 0 is the zero form. If g is 
a second form, 

(28) g = 61X1 4- . . . + b n x n , 

with constant coefficients 61, , b n) we sec that 

(29) / + g = («i 4- bi)xi +... + (<*» + b n )x n 
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Also if c is any constant, we have 

(30) cf = (cai)xi + . . . + (ca n )x n . 

We define —/ to be the form such that / + (—/) = 0 and see that 

(31) -f = -1 • / = (-ai)xi + . . . + (-a n )x n . 

Then / = g if and only if / — g = / + ( — g) — 0, that is a» = (t = 1, 

. . . , n). 

The properties just set down are only trivial consequences of the usual 
properties of polynomials and, as such, may seem to be relatively unimpor¬ 
tant. They may be formulated abstractly, however, as properties of se¬ 
quences of constants (which may be thought of, if we so desire, as the coeffi¬ 
cients of linear forms) and in this formulation in Chapter IV will be very 
important for all algebraic theory. The reader is already familiar with these 
properties which he has used in the computation of determinants by opera¬ 
tions on its rows and columns. 

Let, then, u be a sequence 

(32) u = (cti, . . . , n n ) 

of n constants a» called the elements of the sequence u. If a is any constant, 
we define 

(33) au = ua = (aai, . . . , aa n ) 

and call au the scalar product of u by a. Wo now consider a second sequence, 

(34) v = (< b u . . . , b n ) , 
and define the sum of u and v by 

(35) u + v = («i + 6i, . . . , a n + b n ) . 

Then the linear combination 

(36) au + bv = (aai + bbi, . . . , aa n + bb n ) 

has been uniquely defined for all constants a and b and all sequences 
u and v. 

The sequence all of whose elements are zero will be called the zero 
sequence and designated by 0. It is clearly the unique sequence z with the 
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property that u + z = u for every sequence u. Evidently if a is the zero 
constant au = 0 for every u. 

We define .the negative -u of a sequence u to be the sequence v such 
that u + v — 0. Evidently, then, — u is the unique sequence 

(37) — u = -1 • u = (-ai, . . . , -a n ) , 

and we see that the unique solution of the equation u + x = v is the se¬ 
quence u + ( — u). We evidently call this sequence 

(38) v - u = (bi - a,, . . . , b n - a n ) . 

The reader should observe now that the definitions and properties derived 
for linear combinations of sequences are precisely those which hold for the 
sequences of coefficients of corresponding linear combinations of linear 
forms and that the usual laws of algebra for addition and multiplication 
hold for addition of sequences and multiplication of sequences by constants. 

9. Equivalence of forms. If / = f(x i, . . . , x q ) is a form of degree n in 
X\ 9 ... 9 x q and we replace every x t in / by a corresponding linear form 

(39) x 9 = anyi + . . . + a ir y r (i = 1, . . . , q) , 

we obtain a form g = < 7 ( 1 / 1 , ... 9 y r ) of the same degree n in .t/i, . . . , y T . 
Then we shall say that / is carried into g (or that g is obtained from/) by 
the linear mapping (39). If q = r and the determinant 


( 40 ) 


an 

ai2 

. . . 

a\ q 

a q \ 

a q 2 

. . . 

a qq 


is not zero, we shall say that (39) is nonsingular. In this case it is easily 
seen that we may solve (39) for 7 / 1 , . . . , y q as linear forms in xi 9 ... 9 x q and 
obtain a linear mapping which we may call the inverse of (39). This termi¬ 
nology is justified by the fact that the equation f(x 1 , . . . , x g ) = g(yi 9 
. . . , y q ) is an identity, and thus if we replace y u . . . , y q in <7(1/1, . . . , y q ) 
by the corresponding linear forms in x\ 9 . . . , x q we obtain the original form 

fi.X 1| . . . j Xq). 

We now consider two forms / = f(x 1 , . . . , x q ) and g = g(x 1 , . . . , x q ) of 
the same degree n. Then we shall say that / is eqxdvalent to g if / is carried 
into 0 ( 2 / 1 , . . . , y q ) by a nonsingular linear mapping. The statements above 
imply that if / is equivalent to g then g is also equivalent to /. Thus, we 
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shall usually say simply that / and g are equivalent . We shall not study the 
equivalence of forms of arbitrary degree but only of the special kinds of 
forms described in Section 7, and even of those forms only under restricted 
types of linear mappings. 

We have now obtained the background needed for a clear understanding 
of matrix theory and shall proceed to its development. 

EXERCISES 

1 . The linear mapping (39) of the form x» = yi for i - 1 , . . . , q is called the 
identical mapping. What is its effect on any form// 

2 . Apply a nonsingular linear mapping to carry each of the following forms to an 
expression of the type a\y\ + a 2 y\. Hint: Write / = ai(xi + cx 2 ) 2 + ... by com¬ 
pleting the square on the term in x\ and put X\ + cx 2 = 2 / 1 , x 2 = 2 / 2 . 

a) 2x'i — 4xix 2 + 3x3 d) 2x\ — XiX 2 

b) x\ + 14xiX 2 + e) 3x? + 2xix 2 — x\ 

c ) 3x? + 18xix 2 + 24x5 

3. Find the inverses of the following linear mappings: 

x f 2 xi + x 2 = 2 A f x x + x 2 = 2/1 

' \3xi + 2x 2 = 2/2 ; 1 -2xi + x 2 = 2/2 

4. Apply the linear mappings of Ex. 3 to the following forms / tq obtain equiva¬ 
lent forms g and their inverses to g to obtain /. 


a) f = x\ + x] 


b) f - 4x? — 4xjx 2 + 3x5 



CHAPTER II 


RECTANGULAR MATRICES AND ELEMENTARY 
TRANSFORMATIONS 

1. The matrix of a system of linear equations. The concept of a rec¬ 
tangular matrix may be thought of as arising first in connection with the 
study of the solution of a system 

dnyi + a n y 2 + . . . + a ]n y n = fci , 

G212/1 + do 2 y 2 + . . . + d 2 n y n = k 2 , 


dmiyi ~l“ a m2 y 2 + . . . + a mn y n k m 

of m linear equations in n unknowns y h . . . , y n , with constant coefficients 
a<y. The array of coefficients arranged as they occur in (1) has the form 




and is called the coefficient matrix of the system (1). We shall henceforth 
speak of the coefficients a^ and ki in (1) as scalars and shall derive our theo¬ 
rems with the understanding that they are constants (with respect to the 
variables y lt . . . , y n ) according to the usual definitions and hence satisfy 
the properties usually assumed in algebra for rational operations. In a later 
chapter we shall make a completely explicit statement about the nature 
of these quantities. 

It is not only true that the concept of a matrix arises as above in the 
study of systems of linear equations, but many matrix properties are ob¬ 
tainable by observing the effect, on the matrix of a system, of certain natu¬ 
ral manipulations on the equations themselves with which the reader is 
very familiar. We shall devote this beginning chapter on matrices to that 
study. 

Let us now recall some terminology with which the reader is undoubtedly 
familiar. The line 

XI i = (dil, • • • y din) 
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of coefficients in the ith equation of (1) occurs in (2) as its ith horizontal 
line. Thus, it is natural to call u» the ith row of the matrix A. Similarly, 
the coefficients of the unknowns in (1) form a vertical line 

I “• 1 

(4) u<*' 


\ «m, 


which we call the jth column of A. 

We may now speak of A as a matrix of m rows and n columns, as an 
ra-rowed and n-columned matrix or, briefly, as an m by n matrix. Then 
the rows of A are 1 by n matrices and its columns are m by 1 matrices. We 
shall speak of the scalars as the elements of A, and they may be regarded 
as one by one matrices. The notation a t; which we adopt for the element 
of A in its ith row and jth column will be used consistently, and this usage 
will be of some importance in the clarity of our exposition. To avoid bulky 
displayed equations we shall usually not use the notation (2) for a matrix 
but shall write instead 

(5) A = (a,j) (t = 1, . . . , m; j = 1, . . . , n) . 

If rn = n then A is a square matrix , and we shall speak of A simply as an 
/tt-rowed square matrix. This, too, is a concept and terminology which we 
shall use very frequently. 


ORAL EXERCISES 

1. Read off the elements an, ai 2 , au, a^ in the following matrices 


a) 


c) 



2. Read off the second row and the third column in each of the matrices of Ex. 1. 

3. Read off the systems of equations (1) with constants all zero and matrices 
of coefficients as in Ex. 1. 
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2. Submatrices. In solving the system (1) by the usual methods the 
reader is led to study subsystems of s < m equations in certain t < n of 
the unknowns. The corresponding coefficient matrix has s rows and t col¬ 
umns, and its elements lie in certain $ of the rows and t of the columns of A. 
We call such a matrix an s by t submatrix B of A. If s < m and t < n, 
the elements in the remaining m — s rows and n — t columns form an 
m — s by n — t submatrix C of A and we shall call C the complementary 
sub matrix of B. Clearly, then, B is the complementary submatrix of C. 

It will be desirable from time to time to regard a matrix as being made 
up of certain of its submatrices. Thus we write 

(6) A = (A a) (i = 1 , . . . , j ’ = 1, . . . , t) , 

where now the symbols An themselves represent rectangular matrices. We 
assume that for any fixed i the matrices A a, A i2 ,. . . , A it all have the same 
number of rows, and for fixed k the matrices An, A 2 *, . . . , A, k have the 
same number of columns. It is then clear how each row of A is a 1 by t 
matrix whose elements are rows of An, , A<i in adjacent positions and 
similarly for columns. We have thus accomplished what we shall call the 
partitioning ((>) of A by what amounts to drawing lines mentally parallel to 
the rows and columns of A and between them and designating the arrays 
of elements in the smallest rectangles so formed by A <; . Our principal use 
of ((>) will be the use of the case where we shall regard A as a two by two 
matrix 



whose elements A i, A 2 , A 3 , A 4 are themselves rectangular matrices. Then 
Ai and A 2 have the same number of rows, A 3 and A 4 have the same num¬ 
ber of rows, and every row of A consists partially of a row of Ai and of a 
corresponding row of A 2 or of a row of A 3 and a corresponding row of A 4 . 

Note our usage in ( 2 ), (5), ( 6 ), (7) of the symbol of equality for matrices. 
We shall always mean that two matrices are equal if and only if they are 
identical, that is, have the same size and equal corresponding elements. 

EXERCISES 

1 . State how the columns of A of (7) are connected with the columns of A 1 , A 2 , 
A 3 , and A 4 . 

2. Introduce a notation of an arbitrary six-rowed square matrix A and partition 
A into a three-rowed square matrix whose elements are two-rowed square matrices. 
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Also partition A into a two-rowed square matrix whose elements are three-rowed 
square matrices. 

3. Write out all submatrices of the matrix 

/2 —1 3 4 5\ 

1 0 2 - 1-21 
0 12 3 6 

\7 3 2 1 0/ 

and, if they exist, the complementary submatrices. 

4. Which of the submatrices in Ex. 3 occur in some partitioning of A as a matrix 
of submatrices? 

5. Partition the following matrices so that they become three-rowed square mat¬ 
rices whose elements are two-rowed square matrices and state the results in the 
notation (6). 


2 

-1 

3 

4 

1 

2 \ 

l 

1 

1 

-1 

0 

0 

1 

O' 

0 

1 

1 

-1 

2 

ll 


1 

1 

0 

0 

0 

1 

1 

-2 

1 

0 

0 

0 

b) 


-2 

0 

2 

0 

3 

1 

0 

0 

1 

1 

-1 

3 


0 

2 

0 

2 

-1 

2 

0 

0 

4 

2 

1 

o 


1 

0 

0 

-1 

1 

-4 

0 

10 

0 

-1 

3 

0 

1 

] 

i 

0 

0 

-1 

-1 

0 

-4/ 


6 . Partition the matrices of Ex. 5 into two-rowed square matrices whose elements 
are three-rowed square matrices. 

7. Partition the matrices of Ex. 5 into the form (7) such that .1 1 is a two by three 
matrix; a one by six matrix; a two by two matrix. Read olT A?, A h and A a and 
state their sizes. 

3. Transposition. The theory of determinants arose in connection with 
the solution of the system (1). The reader will recall that many of the prop¬ 
erties of determinants were only proved as properties of the rows of a de¬ 
terminant, and then the corresponding column properties were merely 
stated as results obtained by the process of interchanging rows and columns. 
We call the induced process transposition and define it as follows for mat¬ 
rices. Let A be an m by n matrix, a notation for which is given by (5), 
and define the matrix 

(8) A' = ( g,i) (ga = dij-J = 1, ... ,n;i = 1, ... ,m) , 

which we shall call the transpose of A. It is an n by rn matrix obtained from 
A by interchanging its rows and columns. Thus, the element a», in the tth 
row and jth column of A occurs in A' as the element in its jth row and ith 
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column. Note then that in accordance with our conventions (8) could have 
been written as 

(9) A' = ((ffljj)j-,) O' = 1, • • • , n; i = 1, . . . , m). 

We also observe the evident theorem which we state simply as 

(10) {Ay = A . 

The operation of transposition may be regarded as the result of a certain 
rigid motion of the matrix which we shall now describe. If A is our m by n 
matrix and m < n we put q — in and write 

( 11 ) A = (A u A ,) , 


where A\ is a 7 -rowed square matrix. On the other hand, if n < m, we put 
q = n and have 


( 12 ) 



where the matrix Ai is again a q- rowed square matrix. The line of elements 
an, a 2 2 , . . . , a qq of A and hence of Ay is called the principal diagonal of A 
or, simply, the diagonal of A. It is a diagonal of the square matrix Ay and 
is its principal diagonal. We shall call the an the diagonal elements of A. 
Notice now that A' is obtained from A by using the diagonal of A as an 
axis for a rigid rotation of A so that each row of A becomes a column of A f . 
We should also observe that if .4 has been partitioned so that it has the 
form ((>) then A' is the t by s matrix of matrices given by 


(13) 


A ' — (On) (On = A'ij-yj = 1 , . . . , t; i = 1 , . . . , s). 


We have now given some simple concepts in the theory of matrices and 
shall pass on to a study of certain fundamental operations. 


EXERCISES 

1 . Let A have the form (7). Give the corresponding notation for A'. Give also 
A 9 if A is any matrix of Ex. 1 of Section 1. 

2. In Ex. 1 assume that A = A'. What then is the form (7) of A ? Obtain the 
analogous result if A — —A', where — A is the matrix whose elements are the nega¬ 
tives of those of A. 

3. Let A be a three-rowed square matrix. Find the form (2) of A if A = A' and 
also if A = —A'. 
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4. Prove that the determinant of every three-rowed square matrix A with the 
property that A' = —A is zero. 

5. Solve the system (1) with matrix 



3 

for V t, 2 / 2 , 2/3 in terms of k it A* 3 . Write the results as y i-2 m , and thus eom- 

>» l 

pute the matrix B = (&»•/). Do this also for the system (1) with matrix A ' and com¬ 
pare the results. 

4. Elementary transformations. The system (1) may be solved by the 
method of elimination, and the reader is familiar with the operations on 
equations which are permitted in this method and which yield systems 
said to be equivalent to (1). The resulting operations on the row's of the 
matrix A of the system and corresponding operations on the columns of A 
are called elementary transformations on A and will turn out to be very useful 
tools in the theory of matrices. 

The first of our transformations is the result on the rows of A of (he inter¬ 
change of tw'o equations of the defining system. We define this and the 
corresponding column transformation in the 

Definition 1. Let i r and B be the matrix obtained from A by inter¬ 
changing its ith and r th rows (columns). Then B is said to be obtained from A 
by an elementary rote (column) transformation of type 1. 

The rows (columns) of an m by n matrix are sequences of n (of m) ele¬ 
ments, and the operations of addition and scalar multiplication (i.e., mul¬ 
tiplication by a scalar) of such sequences were defined in Section 1.8.* The 
left members of (1) are linear forms. The addition of a scalar multiple of 
one equation of (1) to another results in the addition of a corresponding 
multiple of a corresponding linear form to another and hence to a corre¬ 
sponding result on the rows of A. Thus we make the following 

Definition 2. Let i and r be distinct integers , c be a scalar , and B be the 
matrix obtained by the addition to the ith row (column) of A of the multiple 
by c of its rth row (column). Then B is said to be obtained from A by an ele¬ 
mentary row (column) transformation of type 2. 

Our final type of transformation is induced by the multiplication of an 

* We shall use a corresjionding notation henceforth when we make references any¬ 
where in our text to results in previous chapters. Thus, for example, by Section 4.7, 
Theorem 4.8, Lemma 4.9, equation (4.10) we shall mean Section 7, Theorem 8, Lemma 9, 
equation (10) in Chapter IV. However, if the prefix is omitted, as, for example, Theorem 8, 
we shall mean that theorem of the chapter in which the reference is made. 



RECTANGULAR MATRICES 


25 


equation of the system (1) by a nonzero scalar a. The restriction o^O 
is made so that A will be obtainable from B by the like transformation 
for cr l . Later we shall discuss matrices whose elements are polynomials 
in x and use elementary transformations with polynomial scalars a. We 
shall then evidently require a to be a polynomial with a polynomial inverse 
and hence to be a constant not zero. In view of this fact we shall phrase 
the definition in our present environment so as to be usable in this other 
situation and hence state it as 

Definition 3. Let the scalar a possess an inverse a -1 and the matrix B be 
obtained as the result of the multiplication of the ith row (column) of A by a. 
Then B is said to be obtained from A by an elementary row (column*) trans¬ 
formation of type 3. 

The fundamental theorems in the theory of matrices are connected with 
the study of the matrices obtained from a given matrix A by the applica¬ 
tion of a finite sequence of elementary transformations, restricted by the 
particular results desired, to A. Thus, it is of basic importance to study 
first what occurs if we make no restriction whatever on the elementary 
transformations allowed. For convenience in our discussion we first make 
the 

Definition. Let A and Bbembyn matrices and let B be obtainable from A 
by the successive application of finitely many arbitrary elementary transforma¬ 
tions. Then we shall say that A is rationally equivalent to B and indicate this 
by writing A ^ B. 

We now observe some simple consequences of our definition. First, we 
see that, if A is rationally equivalent to B and B is rationally equivalent 
to C, the combination of the elementary transformations which carry A to B 
with those which carry B to C will carry A to C. Then A is rationally 
equivalent to C. Observe next that every m by n matrix A is rationally 
equivalent to itself. For the elementary transformations of type 2 with 
c = 0 and of type 3 with a — 1 are identical transformations leaving all 
matrices unaltered. 

Finally, we see that if an elementary transformation carries A to B there 
is an inverse transformation of the same type carrying B to A. In fact, the 
inverse of any transformation of type 2 defined for c is that defined for — c, 
of type 3 defined by a is that defined for a -1 , of type 1 is itself. But then A 
is rationally equivalent to B if and only if B is rationally equivalent to A . 

• The reader should verify the fact that, if we apply any elementary row transforma¬ 
tion to A and then any column transformation to the result, the matrix obtained is the 
same as that which we obtain by applying first the column transformation and then the 
row transformation. 
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Thus we may and shall replace the terminology A is rationally equivalent 
to B in the definition above by A and B are rationally equivalent. 

We have now shown that in order to prove that A and B are rationally 
equivalent it suffices to prove A and B both rationally equivalent to the 
same matrix C. As a tool in such proofs we then prove the following 
Lemma 1 . Let r < m, s < n, and A and B be m by n matrices of the form 



for r by s rationally equivalent matrices Ai and Bi and m — r by n — s ra¬ 
tionally equivalent matrices A* and B*. Then A and B are rationally equivalent. 

For it is clear that any elementary transformation on the first r rows and s 
columns of A induces a corresponding transformation on A \ and leaves A 2 
and the zero matrices bordering it above unaltered. Clearly the sequence 
of such transformations induced by the transformations carrying A\ to B\ 
will replace A by the matrix 



We similarly follow this sequence of elementary transformations by ele¬ 
mentary transformations on the last m — r rows and n — s columns of A 0 
which carry A 2 to B 2 and obtain B. 

It is important also to observe that we may arbitrarily permute the rows 
of A by a sequence of elementary row transformations of type 2, and simi¬ 
larly we may permute its columns. For any permutation results from some 
properly chosen sequence of interchanges. 

Before continuing further with the study of rational equivalence we shall 
introduce the familiar properties of determinants in the language of matrix 
theory and shall also define some important special types of matrices. We 
shall then discuss another result used for the types of proofs mentioned 
above. 

5. Determinants. Let B be the square matrix 
( 14 ) £ = (bn) d, j = 

The corresponding symbol 




bn .. 

bn 

( 15 ) 

D = 

621 • • 

bu 



bn . . 

btt 
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is called a t-rowed determinant or determinant of order t. It is defined as the 
sum of the t! terms of the form 

(16) (-1 ) i bu l b 2 i i ...b tit , 

where the sequence of subscripts i\, . . . , i t ranges over all permutations 
of 1, 2, . . . , t and the permutation i\ y . . . , i t may be carried into 1,2 ,..., 
t by i interchanges. That the sign ( —1)‘ is unique is proved in L. E.. Dick¬ 
son ’s First Course in the Theory of Equations , and we shall assume this re¬ 
sult as well as all the consequent properties of determinants derived there. 

The determinant D will be spoken of here as the determinant of the matrix 
B and we shall indicate this by writing 

(17) D = |fi| 

(read D equals determinant B ). Nonsquare matrices A do not have de¬ 
terminants, but their square submatrices have determinants called the 
minors of A. If A is a square matrix of n > t rows, the complementary sub¬ 
matrix of any J-rowed square submatrix B is an (n — t)-ro wed square ma¬ 
trix whose determinant and that of B are minors of A called complementary 
minors . In particular, every element a i; of a matrix A defines a one-rowed 
square submatrix of A whose determinant is the element itself. Thus we 
have seen that the elements of a matrix may be regarded either as its one-rowed 
square submatrices or as its one-rowed minors. We now pass to a statement 
of some of the most important results on determinants. 

The result on the interchange of rows and columns of determinants men¬ 
tioned in Section 3 may now be stated as 
Lemma 2. Let A be a square matrix. Then 

(18) | A'| = | A | . 

The next three properties of determinants are those frequently used in 
the computation of determinants, and we shall state them now in the lan¬ 
guage we have just introduced. 

Lemma 3. Let B be the matrix obtained from a square matrix A by an ele¬ 
mentary transformation of type 1. Then |B| = — |A|. 

Lemma 4. Let B be the matrix obtained from a square matrix A by an ele¬ 
mentary transformation of type 2. Then | B | = | A |. 

Lemma 5. Let B be the matrix obtained from a square matrix A by an ele¬ 
mentary transformation of type 8 defined for a scalar a. Then |B| =a - |A|. 

The reader will recall that Lemma 3 may be used in a simple fashion 
to obtain 
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Lemma 6 . If a square matrix has two equal rows or columns , its determi¬ 
nant is zero. 

Another result of this type is 

Lemma 7. If a square matrix has a zero row or column , its determinant is 
zero. 

Finally, we have 

Lemma 8. Let A, B, C be n -rowed square matrices such that the ith row 
(column) of C is the sum of the ith row (column) of A and that of B while all 
other rows (columns) of B and C are the same as the corresponding rows (col¬ 
umns) of A. Then 

( 19 ) |C| - IA | + | B | . 

There are, of course, many other properties of determinants, and of these 
we shall use only very few. Those we shall use are, of course, also well 
known to the reader. Of particular importance is that result which might 
be used to define determinants by an induction on order and which does 
yield the actual process ordinarily used in the expansion of a determinant. 
We let A be an n-rowed square matrix A = (a»,) and define to be the 
complementary minor of a,-,. Then the result we refer to states that if we 
define ca = (— then 

n n 

(20) | A | = aikCki — CjkCLkj (i 9 j } = 1, • • • y n) . 

k - 1 k - 1 

Thus, the determinant of A is obtainable as the sum of the products of the 
elements a,-,- in any row (column) of A by their cofactors cthat is, the 
properly signed and labeled minors (— l) <+y d„. 

The result (20) is of fundamental importance in our theory of matrices 
and will be applied presently together with the following parallel result. 
Let B be the matrix obtained from a square matrix A by replacing the ith 
row of A by its gth row. Then B has two equal rows and by Lemma 6 
| B | =0. We expand B as above according to the elements of its ith row 
and obtain as its vanishing determinant the sum of the products of all ele¬ 
ments in the gth row of A by the cofactors of the elements in the ith row 
of A. Combining this result with the corresponding property about col¬ 
umns we have 

n n 

QikCkq — C.kdk) “ 0 
k - 1 k - 1 

(i * q,s* j ; i, j, q, s = 1 , . . . , n) . 


( 21 ) 
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The equations (20) and (21) exhibit certain relations between the arbitrary 
square matrix A and the matrix we define as 

(22) adj A = ( ca) (i, j = 1,. . ., n) . 

These relations will have important later consequences. We call the matrix 
(22) the adjoint of A and see that if A = (a*/) is an n-rowed square matrix 
its adjoint is the n-rowed square matrix with the cofactor of the element 
which appears in the jth row and ith column of A as the element in its 
own ith row and jth column. Clearly, if A = 0 then adj A = 0. 


EXERCISES 

1. Compute the adjoint of each of the matrices 


°) (i i) b) ( 


-1 0 \ 
4 2 


3 0 -1 
-1 2 1 
-3 6 3 


/I 1 -2 

d) (l 1 -3 

\2 3 —4, 


2. Expand the determinants below and verify the following instances of Lem¬ 
ma 8: 


3 

2 

-1 

3 

-3 

-1 3 

-1 

-1 

1 

2 

0 

+ 1 

-1 

0 = 1 

1 

0 

0 

-1 

3 

0 

0 

3 0 

-1 

3 

1 

-1 

4 

-1 

1 

-3 0 

0 

1 

1 

-1 

1 

+ 1 

-1 

1 = 1 

-1 

1 

2 

3 

4 

2 

3 

4 2 

3 

4 


6. Special matrices. There are certain square matrices which have spe¬ 
cial forms but which occur so frequently in the theory of matrices that they 
have been given special names. The most general of these is the triangular 
matrix, that is, a square matrix having the property that either all its ele¬ 
ments to the right or all to the left of its diagonal are zero. Thus a square 
matrix A — (a< 7 ) is triangular if it is true that either a a = 0 for all j > i or 
that a i; = 0 for all j < i. It is clear that A is triangular if and only if A' is 
triangular; and, moreover, we have 

Theorem 1. The determinant of a triangular matrix is the product ana *2 
. . . a nn of its diagonal elements . 

The result above is clearly true if n = 1 so that A — (an), |A| = an. 
We assume it true for square matrices of order n — 1 and complete our 
induction by expanding | A | according to the elements of its first row or 
first column in the respective cases above. 

A matrix A = (a*,) is called a diagonal matrix if it is a square matrix, 
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and a*, = 0 for all i j. Clearly, a diagonal matrix A is triangular so 
that its determinant is the product of its diagonal elements. 

If all the diagonal elements of a diagonal matrix are equal, we call the 
matrix a scalar matrix and have \A | = a^. The scalar matrix for which 
an = 1 is called the n -rowed identity matrix and will usually be designated 
by I. If there be some question as to the order of I or if we are discussing 
several identity matrices of different orders, we shall indicate the order by 
a subscript and thus shall write either /„ or / as is convenient for the n- 
rowed identity matrix. 

Any scalar matrix may be indicated by 

(23) al , 

where a = an is the common value of the diagonal elements of the matrix. 
We shall discuss the implications of this notation later. 

It is natural to call any m by n matrix all of whose elements are zeros 
a zero matrix. In any discussion of matrices we shall use the notation 0 to 
represent not only the number zero but any zero matrix. The reader will 
find that this usage will cause neither difficulty nor confusion. 

We shall frequently feel it desirable to consider square matrices of either 
of the forms 



where A i is a square matrix. Then (24) implies that A 4 is necessarily square, 
and the reader should verify the fact that the Laplace expansion of de¬ 
terminants implies that 

|A| = |A.| = 1^1 • \A<\ . 

The property above and that of Theorem 1 are special instances of a 
more general situation. We let A be a square matrix and partition it as in 
(6) with 8 = t and the submatrices An all square matrices. Then the La¬ 
place expansion clearly implies that if all the A „ are zero matrices for either 
all i > j or all i < j, then \A \ = | A n \ . . . | A u \ . Evidently Theorem 1 
is the case where the A are one-rowed square matrices and the result con¬ 
sidered in (24) the case where t = 2. 

In connection with the discussion just completed we shall define a nota¬ 
tion which is quite useful. Let A be a square matrix partitioned as in (6) 
and suppose that s = t, the An are all square matrices, and every A</ = 0 
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for i ^ j. Then A is composed of zero matrices and matrices A< = An 
which are what we may call its diagonal blocks , and we shall indicate this 
by writing 

(25) A = diag [A u . . . , A t \ . 

As above, the determinant of A is the product of the determinants of its 
submatrices Ai, . . . , A t . 

In closing we note the following result which we referred to at the close 
of Section 4. 

Lemma 9. Every nonzero m by n matrix A is rationally equivalent to a 
matrix 



where I is an identity matrix. 

For by elementary transformations of type 1 we may carry any element 
a pq 0 of A into the element b n of a matrix B which is rationally equivalent 
to A. By an elementary transformation of type 3 defined for a = b ft 1 we 
replace B by a rationally equivalent matrix C with cn = 1. We then apply 
elementary row transformations of type 2 with c — —c r i to replace C by 
the rationally equivalent matrix D = (d„) such that dn = 1, d r \ = 0 for 
r > 1, and then use elementary column transformations of type 2 with 
c = —dir to replace D by the matrix 



Now A i is an m — 1 by n — 1 matrix, and I\ is the identity matrix of one 
row. Clearly, A and A 0 are rationally equivalent. Moreover, either A\ = 0 
and we have (26) for I = h, or our proof shows that Ai is rationally equiv¬ 
alent to a matrix 


(28) 



But then, by Lemma 1, A is rationally equivalent to a matrix 
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After finitely many such steps we obtain (26). 

We shall show later that the number of rows in the matrix I of (26) is 
uniquely determined by A . 


EXERCISES 

1. Carry the following matrices into rationally equivalent matrices of the 
form (26) 


a) 


3 5 

0 1 
-1 -1 


10 4 
2 1 
■ 2 0 


b) 


1 1 2 
2 3 5 
1 2 3 


11 

0 

1 

2 \ 

/-i 

1 

1 

n 

0 

2 

3 1 

« 4 

-3 

-3 

_3 i 

\0 

0 

0 

0 

2 

2 

2 

\l 

0 

1 

1/ 

\ o 

0 

1 

0 / 


c) 


2. Apply elementary row transformations only and carry each of the following 
matrices into a matrix of the form (26) 


a) 


2 1 —2 
1 1 —3 

3 2 —4, 


b) 


10 10 ^ 
2 2-50 
- 21-52 
0 0 0 1 } 


c ) 


3 2 0 0 

2 4 0 0 

4 5 0 0, 


3. Apply elementary column transformations only and carry each of the follow¬ 
ing matrices into a matrix of the form (26). 


/ 3 2 
a) 2 3 
\5 4 





1 0 3 -1\ 
-1 2 4 5/ 



0 10 2 
0 5 13 


) 


4. Show that if the determinant of a square matrix A is not zero then A can be 
carried into the identity matrix by elementary row* transformations alone. Hint : 
The property | A | ^ 0 is preserved by elementary transformations. Some element 
in the first column of A must not be zero, and by row transformations we may 
carry A into a matrix with ones on the diagonal and zeros below and then into 7. 


7. Rational equivalence of rectangular matrices. The largest order of any 
nonvanishing minor of a rectangular matrix A is called the rank of A . The 
result of (20) states that every (l + 1)-rowed minor of A is a sum of nu- 

* Not every matrix may be carried into the form (26) by row transformations only, 
e.g., take A = (1 1). 
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merical multiples of its /-rowed minors, and if the latter are all zero so are 
the former. We thus clearly have 

Lemma 10. Let all (r + 1 )-rowed minors of A vanish. Then the rank of A 
is at most r. 

We may also state this result as 

Lemma 11. Let A have a nonzero r -rowed minor and let all (r + 1 )-rowed 
minors of A vanish. Then r is the rank of A. 

Note that we are assigning zero as the rank of the zero matrix and that 
the rank of any nonzero matrix is at least 1 . 

The problem of computing the rank of a matrix A would seem from our 
definition and lemmas to involve as a minimum requirement the computa¬ 
tion of at least one r-rowed minor of A and all (r + l)-rowed minors. The 
number of determinants to be computed would then normally be rather 
large, and the computations themselves generally quite complicated. How¬ 
ever, the problem may be tremendously simplified by the application of 
elementary transformations. We are thus led to study the effect of such 
transformations on the rank of a matrix. 

Let then A o result from the application of an elementary row transforma¬ 
tion of either type 1 or type 3 to A. By Lemmas 3 and 5 every Crowed 
minor of A 0 is the product by a nonzero scalar of a uniquely corresponding 
/-rowed minor of A, and it follows that A and A 0 have the same rank. If 
A 0 results when we add to the ith row of A the product by c 5 * 0 of its 0 th 
row and B is a /-rowed square submatrix of A, the correspondingly placed 
submatrix B 0 of A 0 is equal to B if no row of B is a part of the ith row of A. 
If, however, a row of B is in the ith row of A and a row of B is in the 0 th 
row of B } then by Lemma 2.4 we have | B 0 1 = | B |. If, finally, a row of B 
is in the ith row of A but no row is in the 0 th row of A, then by Lemma 8 
\Bq\ = | B | + c | C i, where C is a /-rowed square matrix all but one of 
whose rows coincide with those of B, and this remaining row is obtained by 
replacing the elements of B in the ith row of A by the correspondingly 
columned elements in its 0 th row. But then it is easy to see that ±\C\ is 
a minor of A as well as of A 0 . If A has rank r we put t = r + 1 and see 
that |B| = |(?| = 0 , |B 0 1 = 0 for every (r + l)-rowed minor |B 0 | of A 0 . 
Also there exists an r-rowed minor | B | 5 * 0 in A, and our proof shows the 
existence of a corresponding minor | B 0 1 = | B | or | B 0 1 = | B | + c | C | in 
A 0 . But then |B 0 | =0 implies that | C | = — c - ' 1 1B | 5 ^ 0 , and A 0 has a 
nonzero r-rowed minor ± |C|, A and A 0 have the same rank. 

We observe, finally, that if an elementary row transformation be ap¬ 
plied to the transpose A' of A to obtain Ai and the corresponding column 
transformation be applied to A to obtain A 0 , then A' 0 = A\. By the above 
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proof A'and Ai have the same rank, by Lemma 2 A and A' have the same 
minors and hence the same rank, so that Ai and A[ = A 0 have the same 
rank. Hence, A and A 0 have the same rank. Thus, any two rationally 
equivalent matrices have the same rank. 

Conversely, let the rank r of two m by n matrices A and B be the same 
and use Lemma 9 to carry A into a rationally equivalent matrix (26). It 
is clear that the rank of the matrix in (26) is the number of rows in the ma¬ 
trix /. By the above proof A and this matrix have the same rank, I = / r , 
A is rationally equivalent to 



Similarly, B is rationally equivalent to (30) and to A. We have thus proved 
what we regard as the principal result of this chapter. 

Theorem 2. Two m by n matrices are rationally equivalent if and only if 
they have the same rank. 

We have also the consequent 

Corollary. Every m by n matrix of rank r is rationally equivalent to an 
m by n matrix (30). 

A matrix is called nonsingular if it is a square matrix and its determinant 
is not zero. But then Theorem 2 implies, as in Ex. 4 of Section 6, 

Theorem 3. Every n-rowed nonsingular matrix is rationally equivahnl to 
the n-rowed identity matrix. 

In closing let us observe a result of the application to a matrix of either 
row transformations only or column transformations only. We shall prove 

Theorem 4. Every m by n matrix of rank r > 0 may be carried into a 
matrix of the forms 

(31) (q), (H 0), 

respectively , by a sequence of elementary row or column transformations only, 
where G is an r-rowed matrix , II is an r -columned matrix , and both (1 and II 
have rank r. 

For A is equivalent to (30) by a sequence of elementary row and column 
transformations. Clearly, we may obtain (30) by first applying all the row 
transformations and then all the column transformations. If we then apply 
the inverses of the column transformations in reverse order to (30), we ob¬ 
tain the result of the application of the row transformations alone to A. 
But column transformations applied to (30) clearly carry this matrix into 
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a matrix of the form given by the first matrix of (31). Moreover, it is evi¬ 
dent that the rank of this matrix is that of G, G has rank r. The result for 
column transformations is obtained similarly. 


EXERCISES 

1 . Compute the rank r of the following matrices by using elementary transforma¬ 
tions to carry each into a matrix with all but r rows (or columns) of zeros and an 
obvious r-rowed nonzero minor. 


a) 




r) 


1 

4 


-3 4 \ 
-12 16 
- 9 12 / 



-84 

-48 

-12 

-24 


2. Carry the first of each of the following pairs of matrices into the second by 
elementary transformations. Hint: If necessary carry A and B into the form (30) 
and then apply the inverses of those transformations which carry B into (30) to 
carry (30) into B. 


a) A = 

b) A = 

c) A — 


(? 


-1 

-1 


3 ) 

0 / 



-2 

-4 

-6 



-G 



0 

0 

0 

0 

0 

1 

1 

0 




CHAPTER III 

EQUIVALENCE OF MATRICES AND OF FORMS 


1. Multiplication of matrices. If xi, ... , x m and y\, ... ,y n are variables 
related by a system of linear equations 

n 

(1) Xi = (t = 1, . . . , m), 

> - l 

this system was said in Section 1.9 to define a linear mapping carrying the 
Xi to the ifi . We call the m by n matrix A = (an) the matrix of the map - 
ping (1). 

Suppose now that z if . . . , z q are variables related to y u ... , y n by a 
second linear mapping 

Q 

(2) Vi = ^ b ik z k 0=1 , . . . ,n), 

*-l 

with n by q matrix B = (b,*), carrying the y, to the z k . Then, if we sub¬ 
stitute (2) in (1) we obtain a third linear mapping 

Q 

(3) Xi = ^ chZic O' = 1, . . . ,m) , 

k - 1 


with m by q matrix C = (c»jk), and it is easily verified by substitution that 

n 

(4) Cik — 

3 - 1 

~ 1 | • • • j m,) k = l t ... f g) . 

The linear mapping (3) is usually called the product of the mappings ( 1 ) 
and ( 2 ), and we shall also write C = AB and call the matrix C the product 
of the matrix A by the matrix B. 

We have now defined the product AB of an m by n matrix A and an 
n by q matrix B to be a certain m by q matrix C. Moreover, we have de¬ 
fined C so that the element c»* in its ith row and feth column is obtained 

36 
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as the sum c»* = a»ibu + ai 2 b 2 k + . . . + a» n b»* of the products of the 
elements aa in the ith row 

(5) (u*i, d%2i . . . , Uin) 

of A, by the corresponding elements &,* in the kth column 

' bik \ 
b 2 k 

( 6 ) ; 

\ b n k j 

of B. Thus we have stated what we shall speak of as the row by column rule 
for multiplying matrices. 

Observe that if either A or B is a zero matrix the product AB is also. 
Moreover, if A is an m by n matrix, then we have 

(7) 1mA = A/ n = A , 


where I r represents the r-rowed identity matrix defined for every r as in 
Section 2.6. Observe also that I m is the matrix of the linear transforma¬ 
tion Xi = yi for i = l,...,???, and thus in this case the product of (1) 
by (2) is immediately (2); hence (7) is trivially true. 

We have not defined and shall not define the product AB of two matrices 
in which the number of columns in A is not the same as the number of rows 
in B. Then, evidently, the fact that AB is defined need not imply that BA 
is defined. But when both are defined, they are generally not equal and 
may not even be matrices of the same size. This latter fact is clearly so if, 
for example, A is m by ?i, B is n by ???, and ??i n. 

If A and B are n-rowed square matrices, we shall say that A and B are 
commutative if AB = BA. Note the examples of noncommutative square 
matrices in the exercises below. 

Finally, let us observe the following 

Theorem 1. The transpose of a product of two matrices is the product of 
their transposes in reverse order. 

In symbols we state this result as 

(8) (AB)' = B'A'. 


Here A is an m by n matrix, B is an n by q matrix, (AB)' is a g by m ma¬ 
trix, which we state is the product of the q by n matrix B' and the n by m 
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matrix A'. We leave the direct application of the row by column rule to 
prove this result as an exercise for the reader. 


EXERCISES 

1. Compute (3) and hence the product C = AB for the following linear changes 
of variables (mappings) and afterward compute C by the row by column rule. 


f *1 = 2 IJy -b 2/2 

a) j x 2 = 3 j/i - 2/2 

1*3 = 2/1 + 2 2/2 


/ 2/1 = — Z-1+ 2a 

{ 2/2 = -221 + 22 - 323 


I *1 = 2yi + 32/2 — 2/3 
6) ^ x 2 — 2/i 4“ 2/2 ^2/3 

1*3 = 2/2 - 92/3 


f 2/1 = — 132i -f 2 () 2 2 -f 13s 3 

j 2/2 = 9^1 - 1822 - 923 

1 '2/3 = 2i — 22 2 — 23 


2. Compute the following matrix products AB. Compute also BA in the cases 
where the latter product is defined. 


4 3 2\ /—l —2 -3\ 
3 2 1 1 1 3 0 
\2 1 -1/ \-l -2 4/ 



c) 



" 3 \ 

5 /3-2 1 3\ 

-1 [l 14 0 / 
- 2 / 



d) 



2 3 2) 


-1 


3 4) 


3. Let the symbol Ea represent the three-rowed square matrix with unity in the 
ith row and jth column and zeros elsewhere. Verify by explicit computation that 
EijEjk - E ik and that if j j* g then FJijEgk ~ 0. 


2. The associative law. It is important to know that matrix multiplica¬ 
tion has the property 

(9) ( AB)C = A{BC) 

for every m by n matrix A, n by q matrix B, q by s matrix C. This result 
is known as the associative law for matrix multiplication and it may be 
shown as a consequence that no matter how we group the factors in forming 
a product A\ ... A t the result is the same. In particular, the powers A 1 of 
any square matrix are unique. We shall assume these two consequences o f 
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(9) without further proof and refer the reader to treatises on the founda¬ 
tions of mathematics for general discussions of such questions. 

To prove (9) we write A = (a„), B = (ft/*), C = ( c k i ), where in all 
cases i = m; j = 1 , . . . , n; k = 1 , . . . , q; l = 1 , . . . , s. Then it 

is clear that the element in the *th row and li\\ column of G — (AB)C is 

(10) g,i = , 

while that in the same position in H = A(BC) is 

(11) = Vn,/ ’S^bjkCkA . 

Each of these expressions is a sum of nq terms which are respectively of the 
form (■ a l} b ik )cu and a tj (b jk c k i). But those terms in the respective sums with 
the same sets of subscripts are equal, since we have already assumed that 
the elements of our matrices satisfy the associative law a (be) = iab)c. 
Hence, g tl = hu for all i and /, G = H , and (9) is proved. 


EXERCISES 

Compute the products (AB)C and A(BC) in the following cases. 



3. Products by diagonal and scalar matrices. Let A = (a„) be an m by n 

matrix and B be a diagonal matrix, and designate the ith diagonal element 
of B by bi. Then our definition of product implies that if B is w-rowed so 
that BA is defined, then 

BA = (6>«o) 

0 = 1 ,..., m ; j = !,...,«): 


(12) 
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while if B is n-rowed, then 

(13) AB = (anbj) 

(* = 1, . . ., m; j = 1,. . . , n). 

Thus the product of a matrix A by a diagonal matrix B on the left is ob¬ 
tained as the result of multiplying the rows of A in turn by the correspond¬ 
ing diagonal elements of B; the product of A by a diagonal matrix on the 
right is the result of multiplying the columns of A in turn by the corre¬ 
sponding diagonal elements of B. 

Now, let m = n, A be a square matrix. Then from (12) and (13) we see 
that AB = BA if and only if 

(14) (6* - b,)a<, = 0 (i, j = 1, . . . , n) . 

As an immediate consequence of the case bi = . . . = 6 n we have the very 
simple 

Theorem 2. Every n-rowed scalar matrix is commutative with all n -rowed 
square mairices. 

We next see that, if i j* j and bi ^ b,-, then (14) implies that an = 0. 
This gives the result we shall state as 

Theorem 3. Let the diagonal elements of an n-rowed diagonal matrix B be 
all distinct. Then the only n-rowed square mairices commutative with B are 
the n-rowed diagonal matrices. 

We may now prove the converse of Theorem 2—a result which is the 
inspiration of the name scalar matrix. 

Theorem 4. The only n-rowed square matrices which are commutative with 
every n-rowed square matrix are the scalar mairices . 

For let 

(15) B = (bi,) (i, j, = 1, . . . , n) 

and suppose that B is commutative with every n-rowed square matrix A. 
We shall select A in various ways to obtain our theorem. First, we let Ej be 
the diagonal matrix with unity in its jth row and column and zeros else¬ 
where and put BEj = EjB. Equations (12) and (13) imply that the jth 
row of EjB is the same as that of B and the jth column of BEj is the 
same as that of B } while all other columns of BEj are zero. Thus if i ^ j, 
the elements in the ith column of EjB must be zero. Since is in the 
ith column, we have bn = 0 for j 5^ i, and B is a diagonal matrix. If D, 
is the matrix with 1 in its first row and jth column and zeros elsewhere, the 
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product BDj has bn in its first row and jth column and is equal to the 
matrix DjB which has in this same place. Hence, = b n = b for j = 2, 
. . . , n, and B is the scalar matrix bl n . 

Let us now observe that if A is any mby n matrix and a is any scalar, 
then 

(16) (aI m )A = A(al n ) 

is an m by n matrix whose element in the ith row and jth column is the 
product by a of the corresponding element of A. This is then a type of 
product 

(17) aA = Aa 


like that defined in Chapter I for sequences, and we shall call such a product 
the scalar product of a by A. However, we have defined (17) as the in¬ 
stances (16) of our matrix product (4). 


EXERCISES 

1. Compute the products AB and BA by the use of (12) and (13) as well as the 
row by column rule if 


(l 0 0\ 

а) A = 0 -2 0 , 

\0 0 3/ 

/! 0 0 \ 

б ) A = 0 -1 0 , 

\0 0 0 / 


B = 


B = 


2 4 — 5\ 

-4 0 1 

-1 -3 -2/ 

-1 0 — 2 \ 

2 1 41 

-2 3 -3/ 


c) A = 


d) A = 


0 

0 

°\ 

/° 

0 ■ 

-1 

°\ 

2 

0 

o' 

1 D o 

0 

0 

-l\ 

0 

-2 

0| 

B ~ 1 1 

0 

0 

°) 

0 

0 

-2/ 

\o 

1 

0 

0/ 


^3 0 0 
0 2 0 
\0 0 1 


B = 


0 1 0 

0 0 1 
2 0 0 , 


2. Find all three-rowed square matrices B such that BA = AB if 

/I 0 0\ /0 1 0\ 

a) A = l 0-1 0) 6) A = (0 0 1) 

\0 0 2/ \l 0 0/ 
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3. Prove by direct multiplication that BA — —AB if t 2 = —1 and 



4. Let w be a primitive cube root of unity. Prove that BA = wAB if 

/I 0 0 \ /0 1 0\ 

A = / 0 w 0 ) , B = (0 0 1) 

\0 0 W \a 0 0/ 

5. Show that the matrix B of Ex. 4 has the property B 3 = a/ and that the 
matrix A has the property A 3 = /. Obtain similar results for the matrices of Ex. 3. 

6. Compute BAB if 



where c and d are any ordinary complex numbers, c and d are their conjugates. 

4. Elementary transformation matrices. We shall show that the matrix 
which is the result of the application of any elementary row (column) trans¬ 
formation to an m by n matrix A is a corresponding product EA (the prod¬ 
uct AE ) where E is a uniquely determined square matrix. We might of 
course give a formula for each E and verify the statement, but it is simpler 
to describe E and to obtain our results by a device which is a consequence 
of the following 

Theorem 6. Let A be an m by n matrix , B be an n by q matrix so that 
C = AB is an m by q matrix. Apply an elementary row ( column ) transforma¬ 
tion to A (to B) resulting in what we shall designate by A 0 (by B (0) ), and then 
the same elementary transformation to C resulting in C 0 (in C (0) ). Then 

(18) Co = A 0 B , C<°> = AB (0) . 

For proof we see that if we replace by a $i in the right members of (4) 
we get c,k . Thus we obtain the result of an elementary row 7 transformation 
of type 1 on AB by applying it to A. Our definitions also imply that for 
elementary row transformations of type 3 the result stated as (18) follows 
from 

n n 

(19) (aaii)bjk = a a% } bjk ~ ac%k 

j - 1 3 » 1 

(i = 1, . . . , m; k = 1, . . . , q) . 
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Finally, we see that for type 2 they follow from 


(20) ^(a i; + ca 9i )bj k = 

> -1 


^ Q'ijbjk^ "F C ^ — Cik “I" CCgk 

(i = 1, . . . , n; k = 1, . . . , q) , 


and we have proved the first equation in (18). The corresponding column 
result C (0) = AP (0) has an obvious parallel proof or may be thought of as 
being obtained by the process of transposition. It is surely unnecessary to 
supply further details. 

We now write A = IA where 7 is the m-rowed identity matrix and apply 
Theorem 5 to obtain A 0 = 7oA, in which 7 0 is the matrix called E above. 
Then we see that to apply any elementary row transformation to A we 
simply multiply A on the left by either 77,/, 7 J „(c), or Ri(a). Here w r e define 
En to be the matrix obtained by interchanging the ith and jth rows of 
7, Pa(c) by adding c times the jth row of 7 to its ith row, R<(a) by multiply¬ 
ing the ith row of 7 by a 0. We shall call £\„ Pn(c), and Ri{a) elementary 
transformation matrices of types 1, 2, and 3, respectively. 

Observe that 


(21) - E„ = E u , Eifia = 7 , 

so that, if we now assume that E^ is n-rowed, the product AE „ is the result 
of an elementary column transformation of type 1 on A. Similarly, 

(22) [Pa(c)Y = Pn(e) , P, y (-c)P,/(c) = PiMPiA-c) = I , 

and if P,,(c) is w-rowed, then AP Jt (c) is the result of an elementary column 
transformation of type 2 on A. Finally, 

(23) [Ri(a)Y = Ri(a) , Ri(a~ l )Ri(a) = R^R^) = 7, 


and if 7£»(a) is n-rowed, then APj(a) is the result of an elementary column 
transformation of type 3 on A. Thus the elementary column transforma¬ 
tions give rise to exactly the same set of elementary transformation mat¬ 
rices as were obtained from the row transformations. 

We shall now interpret the results of Section 2.7 in terms of elementary 
transformation matrices. First of all, we may interpret Theorem 2.2 as the 
following 
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Lemma 1. Let A and B be m by n matrices. Then there exist elementary 
transformation matrices Pi, ... , P s , Qi, . . . , Q t such that 

B = (Pi. . . P.)A(Qi. . . Q t ) 

if and only if A and B have the same rank. 

Theorem 2.3 is the case of Lemma 1 where A is the n-rowed identity 
matrix, and consequently B = Pi. . . P.Qi. . . Q t . Thus we obtain 

Lemma 2. Every nonsingular matrix is a product of elementary transforma¬ 
tion matrices. 

We shall close the results with an important consequence of Theorem 2.4. 

Theorem 6. The rank of a product of two matrices does not exceed the rank 
of either factor. 

For, by Theorem 2.4, if the rank of A is r, there exists a sequence of ele¬ 
mentary transformations carrying A into an m by n matrix Aq whose bot¬ 
tom m — r rows are all zero. By Theorem 5, if we apply these transforma¬ 
tions to C = AB, we obtain a rationally equivalent matrix C 0 = A Q B. Then 
the bottom m — r rows of the ra-rowed matrix Co are all zero, and the 
rank of C 0 is the rank of C and is at most r, as desired. The corresponding 
result on the relation between the ranks of C and B is obtained similarly. 

EXERCISES 

1. Express the following as products of elementary transformation matrices. 



2. Find the elementary transformation matrices corresponding to the elementary 
row transformations used in Ex. 2 of Section 2.6 and carry the matrices of that exer¬ 
cise into the form (2.30) by matrix multiplication. 

6. The determinant of a product. In the theory of determinants it is 
shown that the symbol for the product of two determinants may be com¬ 
puted as the row by column product of the symbols for its factors. In our 
present terminology this result may be stated as 

Theorem 7. The determinant of a product of two square matrices is the 
product of the determinants of the factors , that is, 

(24) I AB | = | A | • | B | . 

The usual proof in determinant theory of the result is quite complicated, 
and it is interesting to note that it is possible to derive the theorem as a 
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simple consequence of our theorems which were obtained independently 
and which we should have wished to derive even had Theorem 7 been as¬ 
sumed. We shall, therefore, give such a derivation. We thus let A and B 
be n-rowed square matrices and see that Theorem 6 states that if \A | =0 
then | AB | =0. Hence, let A be nonsingular so that, by Lemma 2, 

A = . . . P ., 

where the P% are elementary transformation matrices. From our defini¬ 
tions we see that, if E is an n-rowed elementary transformation matrix of 
type 1, 2, or 3, then \E\ = — 1, 1, or a, respectively. Thus, if G is an 
n-rowed square matrix and Go = EG, then Lemmas 2.3, 2.4, and 2.5 imply 
that |G 0 | = — |Gj, | G |, a |G|, respectively, and hence |G 0 | = \E\ • |G|. 
It follows clearly that, if Ei, . . ., E t are any elementary transformation 
matrices, then \E t E t -i . . . EiG\ = \E t \ ... |Ei| • |G|. We apply this re¬ 
sult first to A to obtain \A\ = |Pi| ... \P g \ and then to AB to obtain 
\AB\ = |Pi . . . P.B | = |Pi| . . . |P.| • |P| = |A| • |B| as desired. 

6 . Nonsingular matrices. An n-rowed square matrix A is said to have an 
inverse if there exists a matrix B such that AB = BA = J is the n-rowed 
identity matrix. Clearly, B is an n-rowed square matrix which we shall 
designate by A~ l and thus write 

(25) AA~ l = A~ l A = I . 

Moreover, we have 

Theorem 8. A square matrix A has an inverse if and only if A is non- 
singular . 

For if (25) holds, we apply Theorem 7 to obtain |A| • |A -1 | = \I\ = 1, 
| A | 5 ^ 0. The converse may be shown to follow from (21), (22), (23), and 
Lemma 2; but we shall prove instead the result that if | A | ^ 0 then A -1 is 
uniquely determined as the matrix 

(26) A" 1 = | A | - 1 • adj A . 

This formula follows from the matrix equation 

(27) A(adj A) - (adj A)A - | A | • I 

by multiplication by the scalar | A I” 1 , and we observe that (27) is the in¬ 
terpretation of (2.20) and (2.21) in terms of matrix product, where adj A is 
our symbol for the adjoint matrix defined in (2.22). 
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We now prove A" 1 unique by showing that if either AB = 7 or BA = 7 
then B is necessarily the matrix A -1 of (26). This is clear since in either 
case \A \ 0, A -1 of (26) exists, A" 1 = A~ l I = A~ l (AB) = (A -1 A)B = 

IB = B, and similarly if BA = 7. 

We note also that, if A and B are nonsingular, 

(28) (AB)" 1 = B~ l A ~ l . 

For (AB)(B" 1 A“ l ) = A(BB -1 )A" 1 = A A -1 = 7. Finally, if A is nonsingu¬ 
lar we have 

(29) (A" 1 )' = (A')“ l . 

For 7' = 7 = (AA- 1 )' = (A~ l )'A', (A -1 )' is the inverse of A'. 

A linear mapping (1) with m = n has an inverse mapping (2) with 
n = q if the product (3) is the identical mapping, that is, if C = AB is the 
identity matrix. But this is then possible if and only if | A | ^ 0, that is, 
(1) is what we called a nonsingular linear mapping in Section 1.9. We shall 
use this concept later in studying the equivalence of forms. 


EXERCISES 

1. Show that if A is an n-rowed square matrix of rank n - 1 ^ 0, then adj A 
has rank 1. Hint: By (27) we have A (adj A) = 0, PAQ((f~ l • adj A) = Ofor 



Then the first n — 1 rows of Q~ l • adj A must be zero, adj A has rank 1. 

2. Use the result above and (28) to prove that adj (adj A) — \A | n- M if n > 2 
and adj (adj A) — A if n = 2 and | A | ^ 0. 

3. Use Ex. 1 and (27) to show that |adj A| = | A | n ~ 


4. Compute the inverses of the following matrices by the use of (26). 



be — 1 
c 


/0 0 2 \ 
b) ( 1 0 0 ) 
\0 1 0/ 



d) 


2 3 4 
1 2 6 
0 0 1 


(! i 5) 

\l 3 6/ 


( 1-10 2 \ 
0 1 1 -l\ 

2 12 1 I 

3-2 1 6/ 


( 2-11 6 \ 
0 1 2 l\ 

-1 1 1 -2 
1 0 2 3/ 



2 

3 

-1 

0 


-3 

-1 

0 

1 
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5. Let A be the matrix of ( d ) above and B be the matrix of (e). Compute C = 
by (28) and verify that (AB)C — I by direct multiplication. 

7. Equivalence of rectangular matrices. Our considerations thus far have 
been devised as relatively simple steps toward a goal which we may now 
attain. We first make the 

Definition. Two m by n matrices A and B are called equivalent if there 
exist nonsingular matrices P and Q such that 

(30) PAQ = B . 

Observe that P and Q are necessarily square matrices of m and n row T s, 
respectively. By Lemma 2 both P and Q are products of elementary trans¬ 
formation matrices and therefore A and B are equivalent if and only if A and 
B are rationally equivalent. The reader should notice that the definition of 
rational equivalence is to be regarded here as simply another form of the 
definition of equivalence given above and, while the previous definition is 
more useful for proofs, that above is the one which has always been given 
in previous expositions of matrix theory. We may now apply Lemma 1 and 
have the principal result of the present chapter. 

Theorem 9. Two m by n matrices are equivalent if and only if they have 
the same rank. 

We emphasize in closing that, if A and B are equivalent, the proof of 
Theorem 2.2 shows that the elements of P and Q in (30) may be taken to 
be rational functions, with rational coefficients, of the elements of A and B. 

EXERCISES 

1. Compute matrices P and Q for each of the matrices A of Ex. 1 of Section 2.6 
such that PAQ has the form (2.30). Hint: If A is m by n t we may obtain P by ap¬ 
plying those elementary row transformations used in that exercise to I m and simi¬ 
larly for Q. (The details of an instance of this method are given in the illustrative 
example at the end of Section 8.) 

2. Show that the product AB of any three-rowed square matrices -4 and B of 
rank 2 is not zero. Hint: There exist matrices.P and Q such that A 0 = PAQ has the 
form (2.30) for r = 2. Then, if AB = 0, we have AqB 0 = 0 where B 0 = Q~'B has 
the same rank as B and may be shown to have two rows with elements all zero. 

3. Compute the ranks of A, B t AB for the following matrices. Hint: Carry A 
into a simpler matrix Aq — PA by row transformations alone, B into B 0 = BQ by 
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column transformations alone, and thus compute the rank of AoBq = P(AB)Q in¬ 
stead of that of AB. 



/l 2 3\ 

a) A = 

2 3 4 - 


\1 3 5/ 


/-I 2 0\ 

b) A = 

1 1 0 ■ 


\—4 5 0/ 


B = 


3 4 5 
0 2 4 
-1 1 3 


B = 


1 3 2\ 

2 6 4) 

-1 -4 1/ 


e) 



2 

1 

0 

3 



B = 


( -4 10 2 6\ 

-2 6 2 4 \ 

3-4 2-1 
-12 0 1 / 


8. Bilinear forms. As we indicated at the close of Chapter I, the problem 
of determining the conditions for the equivalence of two forms of the same 
restricted type is customarily modified by the imposition of corresponding 
restrictions on the linear mappings which are allowed. We now precede 
the introduction of those restrictions which are made for the case of bi¬ 
linear forms by the presentation of certain notations which will simplify our 
discussion. 

The one-rowed matrices 


(31) x’ = X m ) , y' = (i y u . . . , y n ) 


have one-columned matrices x and y as their respective transposes. We 
let A be the m by n matrix of the system of equations (1) and see that this 
system may be expressed as either of the matrix equations 

(32) x = Ay , x' = y'A'. 

We have called (1) a nonsingular linear mapping if m = n and A is non¬ 
singular. But then the solution of (1) for y if . . . , y n as linear forms in 
xi, . . . , x n is the solution 

(33) y « x , y' = x'(A')* 1 - x'(A -1 )' 

of (32) for y in terms of x (or y' in terms of x'). We shall again consider 
m by n matrices A and variables x< and yj and shall now introduce the no¬ 
tation 


( 34 ) 


x = P'u 
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for a nonsingular linear mapping carrying the x < to new variables u k , for i, 
k = 1, . . . , ra, so that the transpose P of P* is a nonsingular m-rowed 
square matrix and u' = (u h . . . , u m ). Similarly, we write 

(35) y = Qv 

for a nonsingular n-rowed square matrix Q and v' = (vi, . . . , v n ). We now 
return to the study of bilinear forms. 

A bilinear form / = 2 Xia^y^ for i = 1, . . . , m and j = 1, . . . , n, is a 
scalar which may be regarded as a one-rowed square matrix and is then the 
matrix product 

(36) / = x’Ay . 

Here x and y are given by (31), and we call the m by n matrix A the matrix 
of/, its rankle rank off. Also let g — x'By be a bilinear form in xi, . . . , x m 
and 2 / 1 , . . . , y n with m by n matrix B. Then we shall say that / and g are 
equivalent if there exist nonsingular linear mappings (34) and (35) such 
that the matrix of the form in Ui, . . . , u m and v h . . . , v n into which / is 
carried by these mappings is B. But if (34) holds, then x ' = u'P and 

/= (u'P)A(Qr) = u'(PAQ)v, 

so that B = PAQ and A are equivalent. Thus, two bilinear forms f and g 
are equivalent if and only if their matrices are equivalent . By Theorem 9 we 
see that two bilinear forms are equivalent if and only if they have the same 
rank. It follows also that every bilinear form of rank r is equivalent to 
the form 

(37) Xit/i + . . . + x r y r . 

These results complete the study of the equivalence of bilinear forms. 
ILLUSTRATIVE EXAMPLE 

We shall find nonsingular linear mappings (34) and (35) which carry the form 
/ = 2xiyi - 3*12/2 + xiy 9 - xty x + bx^y z - 6x 5 yi + Zx z y 2 + I9x z y z 
into a form of the type (37). The matrix of / is 


A = 


2 -3 V 
-10 5 

-6 3 19, 
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We interchange the first and second rows of A f add twice the new first row to the 
new second row, add —6 times the new first row to the third row, and obtain 



which evidently has rank 2. We then add the second row to the third row, multiply 
the first row by — 1, the second row by — J, and obtain 



The matrix P is obtained by performing the transformations above on the three- 
rowed identity matrix, and hence 


P = 




We continue and carry PA into PAQ of the form (2.30) for r = 2 by adding five 
times the first column and V times the second column of PA to its third column. 
Then 


1 0 

0 1 

,0 0 



The corresponding linear mappings (34) and (35) are given respectively by 


I Xi = ~iu2 + Ut 
x 2 = —ui — Jt h — 4 u, , 
x 3 = u 2 


1 2 /i = vi + 5 v t 
2 /i = v 2 + • 

2/i = v z 


We verify by direct substitution that 


/ = (“}ui + Ui)(2vi + lOvj — 3t>2 — lit;, + Vi) 

+ (ui + §u 2 + 4ui)(v\ + 5t;, - 5t>,) 

4 u,(— 6 vi — 30 t;, -f- 3 vi -f* lit;, 4 19 t;,) 

= i + w$2 + 2uiV\ — 3uiV 2 + UiVi 4 iUiVi 4 4u,t;i — 6u,t;i 4 3utv 2 

— UiVi 4 Utf ) 2 , 


as desired. 


EXERCISES 

1. Use the method above to find nonsingular linear mappings which carry the 
following bilinear forms into forms of the type (37). 

a) 2x x yi 4 3x4/* — x#i — 2xjy 2 

b) xii/i - Xij/j 4 xijji + 4x22/i - 3x22/a 4 2 x 22/1 - x, 2 /i + 2x z y 2 — 3 x,y« 
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c) 2x x y x + 3x x y 2 — x x y z -f bx 2 y x + 2x 2 y 2 + x 2 y z -f 3x z y x — x z y 2 + 2x z y z 

d) 2x x yi — x x y 2 + 3x x y z 4- x x y 4 — 2x 2 y\ + 3x 2 y 2 — \3x 2 y z + 5 x 2 y A + x z y x 

- x z y z + 2x z y A 4* 4x 4 2/i — x K y 2 4- x<y z -f 5 x A y A 

2. Use the method above to find a nonsingular linear mapping on x Xf x 2 , x z such 
that it and the identity mapping on y it y 2i y z carry the following forms into forms 
of the type (37). 

a) —3xi2/i + 5xi2/ 2 4- x 2 y x — 2x 2 y 2 

b) 3xi y x + x x y 2 -f 4x 2 t/i 4- 2 x 2 y 2 

c) x x y x — 4xi2/ 2 4- x x y z — x 2 y x + x 2 y 2 -f 2 x z y 2 — x z y z 

d) 3x x y 2 4- 5xi y z -|- x 2 y 2 4- 2x*y z + x z y x - 2 x z y 2 + x z y z 

e ) -x x y 2 4- x x y z + x x y 4 4* x 2 y 2 4- 2x 2 y z - 3 x 2 y A 4- x z y x — 2x z y 2 4- 3x z y z 

4- 4 x z y K 4- 2 x 4 2/2 - x^ z - 3 x K y A 

3. Find a nonsingular linear mapping on y X) y 2} y z such that it and the identical 
mapping on x Xt x 2 , x 3 carry the forms of Ex. 2 into forms of the type (37). Hint: The 
matrices A of the forms of Ex. 2 are nonsingular so that the corresponding matrices 
P have the property PA = /. Then P = A~ l , AQ = I has the unique solution 

Q - P- 


4. Use elementary row transformations as above to compute the inverses of the 
matrices of Section 6, Ex. 4. 


5. Use elementary column transformations to compute the inverses of the fol¬ 
lowing matrices. 



c) 





9. Congruence of square matrices. There are some important problems 
regarding special types of bilinear forms as well as the theory of equivalence 
of quadratic forms which arise when we restrict the linear mappings we use. 
We let m = n so that the matrices A of forms/ = x f Ay are square matrices. 
Then (34) and (35) are called cogredient mappings if Q = P '. Then / and g 
are clearly equivalent under cogredient mappings if and only if 

(38) B = PAP f . 

We shall call A and B congruent if there exists a nonsingular matrix P 
satisfying (38). 
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We shall not study the complicated question as to the conditions that 
two arbitrary square matrices be congruent but shall restrict our attention 
to two special cases. 

Before passing to this study we observe that Lemma 2 and Section 7 
imply that A and B are congruent if and only if A may be carried into B by 
a sequence of operations each of which consists of an elementary row trans¬ 
formation followed by the corresponding column transformation. Thus 
P = P a . . . Pi, where the P* are elementary transformation matrices, B = 
P. . . . P X AP[ . . . P.' = P. . . . P 2 (PiAP1)P' 2 . . . P.', and so forth as de¬ 
sired. We shall speak of such operations on a matrix A as cogredient ele¬ 
mentary transformations of the three types and shall use them in our study 
of the congruence of matrices. 

10. Skew matrices and skew bilinear forms. A square matrix A is called 
symmetric if A' = A, and skew if A' = —A. If B = PAP' is congruent to 
A then B ' = PA'P', B is symmetric if and only if A is symmetric, B is 
skew if and only if A is skew. 

If A = (an) is a skew matrix, then a,, = — a»y for all i and j. Hence 
an = — an and consequently every diagonal element of A is zero. We use 
this result in the proof of 

Theorem 10. Two n-rowed skew matrices are congruent if and only if they 
the same rank r. Moreover r is an even integer 2t, and every skew matrix is 
thus congruent to a matrix 

/0 ~I t 0\ 

(39) It 0 0 . 

\0 0 0/ 

For either A = 0 = PAP' for every P and our result is trivial, or some 
an 0. We may interchange the ith row and first row, the jth and sec¬ 
ond row and thus also the corresponding columns by cogredient elementary 
transformations of type 1. We thus obtain a skew matrix H = (A„) con¬ 
gruent to A and with an = h\ 2 0, h 2 \ = — h\ 2 , hn = A 22 = 0. Multiply 

the first row and column of H by hi i 1 and obtain a skew matrix C = (c» 7 ) 
congruent to A and with c i2 = —1, c 2 \ = 1, c n = C 22 = 0. We now apply 
a sequence of cogredient elementary transformations of type 2 where we 
add — cn times the second row of C to its jth row as well as cy* times the 
first row of C to its jth row for j = 3, . . . , n, and thus obtain the skew 
matrix 
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The matrix A 0 is congruent to A, Ai is necessarily a skew matrix of n — 2 
rows, and any cogredient elementary transformation on the last n — 2 rows 
and columns of A 0 induces corresponding transformations on Ai, carrying 
it to a congruent skew matrix H\ and carrying A 0 to a congruent skew 
matrix 



It follows that, after a finite number of such steps, we may replace A by 
a congruent skew matrix 

(41) G = diag { E\ , . . . , E t} 0 , . . . , 0 } 

with each E k = E. Clearly, G and A have rank 2 1 . If B also is a skew ma¬ 
trix of rank 2 1 , then B is congruent to (41) and to A, both A and B are con¬ 
gruent to (39). 

If A is a skew matrix, the corresponding bilinear form x'Ay is a skew 
bilinear form. Hence, two such forms are equivalent under cogredient non¬ 
singular linear mappings if and only if they have the same rank. Moreover, 
if / is a skew bilinear form of rank 2 1 it is equivalent under cogredient non¬ 
singular linear mappings to (xi y 2 - x 2 yi) + . . . + {x 2 t -iy 2i - x 2 t y 2 t -i). 

EXERCISES 

Use a method analogous to that of the exercises of Section 8 to find cogredient 
linear mappings carrying the following skew forms to forms of the type above. 

a) x<y 2 - x 2 yi + 2xiy 3 - 2x 3 ?/i + 3x^3 — 3 x 3 y 2 

b ) 2x21/1 - 2x x y 2 - X1I/3 + x 3 y x + 2x 2 y 3 - 2x 3 y 2 

c) xiy 2 - x 2 i/i + Sx x y 3 - Sx 3 yi - 4x<yi + 4xi+ 4x^3 - 4 x 3 i/2 

+ 8x 2 i/ 4 - 8x4j/2 + Sx 3 y A — 8x4/3 

11. Symmetric matrices and quadratic forms. The theory of symmetric 
matrices is considerably more extensive and complicated than that of skew 
matrices, and we shall obtain only some of its most elementary results. Our 
principal conclusion may be stated as 

Theorem 11. Every symmetric matrix A is congruent to a diagonal matrix 
of the same rank as A. 

We may evidently assume that A = A f ^ 0, and we shall prove first 
that A is congruent to a symmetric matrix H = (ha) with some diagonal 
element ha 5 * 0 . This is true for A = H if some diagonal element of A is 
not zero. Otherwise there is some a % , 5 ^ 0 , a,< = a*„ and an = a,-,- = 0 . We 
then obtain H as the result of adding the jth row of A to its ith row and 
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the jth column of A to its ith column, ha = a,< + an = 2a u ^ 0 as de¬ 
sired. We now permute the rows and corresponding columns of H to ob¬ 
tain a congruent symmetric matrix C with Cn = ha ^ 0. Then we add 
—CnCkx times the first row to its fcth row, follow with the corresponding col¬ 
umn transformation, and obtain a symmetric matrix 



congruent to A. Clearly, A\ is a symmetric matrix with n —* 1 rows. As 
in the proof of Theorem 10 we carry out a finite number of such steps, and 
it is clear that we ultimately obtain a diagonal matrix. It is a matrix equiv¬ 
alent to A and must have the same rank. 

The result above may be applied to obtain a corresponding result on 
symmetric bilinear forms of rank r, that is, bilinear forms / = x'Ay defined 
by a symmetric matrix A of rank r. Theorem 11 then states that/is equiva¬ 
lent under cogredient transformations to a form 

(43) aiXiyi + . . . + a r x r y r • 

The results of Theorem 11 may also be applied to quadratic forms / = 
f(x i, . . . , x n ). As we saw in Section 1.7, / is the one-rowed square matrix 
product 

n 

(44) / = x'Ax = ^ XjajjXj 

i,3 - 1 

for a symmetric matrix A. We call the uniquely determined symmetric 
matrix A the matrix of f and its rank the rank of f. Now in Section 1.9 we 
defined the equivalence of any quadratic form / and any second quadratic 
form g = y'By . We may then use the notations developed above and see 
that / and g are equivalent if and bnly if A and B are congruent. For if 
our nonsingular linear mapping is represented by the matrix equation x = 
P'u, then x f = u'P is a consequence, and / = u'(I>AP')u,f and g are equiv¬ 
alent if and only if B = PAP'. Thus Theorem 11 states that every quad¬ 
ratic form of rank r is equivalent to 

(45) a x x\ + . . . + a r x \. 

The form (45) is of course to be regarded as a form in x\ t ... t x n with 
matrix diag {<zi, . . . , a f , 0, . . . , 0}. However, it may be regarded as a 
form in £i, . . . , x r with nonsingular matrix diag [ai f ... f a r \. We shall 
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call a quadratic form / = x'Ax in Xi, . . . , x n with nonsingular matrix A 
a nonsingular form, and we have shown that every quadratic form of rank r in 
n variables may be written as a nonsingular form in r variables whose matrix 
is a diagonal matrix . 


EXERCISES 

1. What is the symmetric matrix A of the following quadratic forms? 

a ) 3z? + x &2 + 2xiX 3 — 3x 2 x 4 + x\ 

b) 2 x\X 2 — 3xjX 4 + x\ 

c ) x\ — 3xix 2 + 2x\x 4 — 3x\ 


2 . Find a nonsingular matrix P with rational elements for each of the following 
matrices A such that PAP’ is a diagonal matrix. Determine P by writing A = 
IAP and by applying cogredient elementary transformations. 


a) 




b) 


2 4 
4 9 
1 4 


c) 




/ 10 -3 3\ 

d) ( - 3 1-1 

V 3 -1 1/ 


e) 


/ 0 

0 

2 

3 \ 


/ 7 

f> 

0 

— 1 \ 

1° 

0 

-2 

_3 i 

/) 

1 6 

15 

0 

1 

U 

-2 

0 

-1 

0 

0 

2 

3 


[3 -3 -1 —3/ 


,-l 1 3 


g) 






3. Write the symmetric bilinear forms whose matrices are those of (a), (6), (c), 
and ( d ) of Ex. 2 and use the cogredient linear mappings obtained from that exercise 
to obtain equivalent forms with diagonal matrices. 

4. Apply the process of Ex. 3 to (e), (/), ( 0 ), and ( h ) of Ex. 2 for quadratic forms. 

5. Which of the matrices of Ex. 2 are congruent if we allow any complex num¬ 
bers as elements of Pf 

6. Show that the forms / = x\ + x\ and g « x\ — x\ are not equivalent under 
linear mappings with real coefficients. Hint: Consider the possible signs of values 
of / and of 0. 
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12. Nonmodular fields. In our discussion of the congruence and the 
equivalence of two matrices A and B the elements of the transformation 
matrices P and Q have thus far always been rational functions, with rational 
coefficients, of the elements of A and B. While we have mentioned this 
fact before, it has not, until now, been necessary to emphasize it. But the 
reader will observe that we have not, as yet, given conditions that two sym¬ 
metric matrices be congruent, and our reason is that it is not possible to do 
so without some statement as to the nature of the quantities which we allow 
as elements of the transformation matrices P. We shall thus introduce an 
algebraic concept which is one of the most fundamental concepts of our 
subject—the concept of a field. 

A field of complex numbers is a set F of at least two distinct complex 
numbers a, 6, ... , such thht a + b, ab, a — b f and a/c are in F for every 
a, b } c 5 * 0 in F. Examples of such fields are, then, the set of all real num¬ 
bers, the set of all complex numbers, and the set of all rational functions 
with rational coefficients of any fixed complex number c. 

If F is any field of complex numbers, the set K = F(x) of all rational 
functions in x with coefficients in F is a mathematical system having prop¬ 
erties, with respect to rational operations, just like those of F. Now it is 
true that even if one were interested only in the study of matrices whose 
elements are ordinary complex numbers there would be a stage of this study 
where one would be forced to consider also matrices whose elements are 
rational functions of x. Thus we shall find it desirable to define the concept 
of a field in such a general way as to include systems like the field K defined 
above. We shall do so and shall assume henceforth that what we called con¬ 
stants in Chapter I and scalars thereafter are elements of a fixed field F. 

The fields we have already mentioned all contain the complex number 
unity and are closed with respect to rational operations. But it is clearly 
possible to obtain every rational number by the application to unity of a 
finite number of rational operations. Thus all our fields contain the field 
of all rational numbers and are what are called nonmodular fields. The fields 
called modular fields will be defined in Chapter VI. We now make the fol¬ 
lowing brief 

Definition. A set of elements F is said to form a nonmodular field if F 
contains the set of all rational numbers and is closed with respect to rational 
operations such that the following properties hold : 

I. (a + b) + c = a + (b + c) , (ab)c = a(bc) 

II. a(b + c) = ab + ac ; 

III. a + b = b + a , ab = ba 
for every a, b, c of F. 
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The difference a — 6 is always defined in elementary algebra to be a solu¬ 
tion x in F of the equation 

(46) x + 6 = a , 


and the quotient a/6 to be a solution y in F of the equation* 

(47) yb = a . 

Thus our hypothesis that F is closed with respect to rational operations 
should be interpreted to mean that any two elements a and 6 of F determine 
a unique sum a + b and a unique product ab in F such that (46) has a 
solution in F and (47) has a solution in F if 6 5 * 0. In the author’s Modern 
Higher Algebra it is shown that the solutions of (46) and (47) are unique. 
In fact, it may be concluded that the rational numbers 0 and 1 have the 
properties 

(48) a + 0 = al=a, aO = 0, 

for every a of F and that there exists a unique solution x = — a of x + a = 
0 and a unique solution y = b~ l of yb = 1 for 6^0. Then the solutions 
of (46) and (47) are uniquely determined by x = a + ( — 6 ), y = a6 _1 . 

We also see that the rational number —1 is defined so that ( —1) + 
1 = 0 , and thus (— 1 + l)a = 0 • a = 0 , whereas ( —1 + l)a = — 1 • a + 
1 • a = — 1 • a + a. Hence —a = — 1 • a. It is also true that —( — a) = 
a, ( — a) ( — 6) = ab for every a and 6 of a field F. 


EXERCISES 


1. Let a, 6, and c range over the set of all rational numbers and F consist of all 
matrices of the following types. Prove that F is a field. (Use the definition of ad¬ 
dition of matrices in (52).) 



»(•:* r-) 



2. Show that if a, 6, c, and d range over all rational numbers and i 2 = — 1, the 
set of all matrices of the following kind form a quasi-field which is not a field. 

/a -}- bi 3(c -f- di)\ 

\c — di a — bi / 

* Note that in view of III the property II implies that (6 + c)a = ba + ca, and the 
existence of a solution of (46) is equivalent to that of b + x = a, of (47) to that of 
by =* o. But there are mathematical systems called quasi-fields in which the law ab — ba 
does not hold, and for these systems the properties just mentioned must be made as 
additional assumptions. 
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13. Summary of results. The theory completed thus far on polynomials 
with constant coefficients and matrices with scalar elements may now be 
clarified by restating our principal results in terms of the concept of a field. 
We observe first that if f(x) and g{x) are nonzero polynomials in x with 
coefficients in a field F then they have a greatest common divisor d(x) with 
coefficients in F, and d{x) = a(x)f(x) + b(x)g(x) for polynomials a(x) y b(x) 
with coefficients in F . 

We next assume that A and B are two m by n matrices with elements in 
a field F and say that A and B are equivalent in F if there exist nonsingular 
matrices P and Q with elements in F such that PAQ = B. Then A and B 
are equivalent in F if and only if they have the same rank. Moreover, corre¬ 
spondingly, two bilinear forms with coefficients in F are equivalent in F if 
and only if they have the same rank. Since the rank of a matrix (and of a 
corresponding bilinear form) is defined without reference to the nature of 
the field F containing its elements, the particular F chosen to contain the 
elements of the matrices is relatively unimportant for the theory. 

If A and B are square matrices with elements in a field F, then we call A 
and B congruent in F if there exists a nonsingular matrix P with elements 
in F such that PAP ' = B. Similarly, we say that the bilinear forms x'Ay 
and x'By are equivalent in F under cogredient transformations* if A and B 
are congruent in F. When A' = — A the matrix A is skew, and every ma¬ 
trix B congruent in F to A is skew, two skew matrices with elements in F 
are congruent in F if and only if they have the same rank. Hence, the pre¬ 
cise nature of F is again unimportant. 

Let A = A' be a symmetric matrix with elements in F so that any ma¬ 
trix B congruent in F to A also is a symmetric matrix with elements in F. 
Then two corresponding quadratic forms x f Ax and x'Bx are equivalent in 
F if and only if A and B are congruent in F. Moreover, we have shown 
that every symmetric matrix of rank r and elements in F is congruent in F 
to a diagonal matrix diag {<zi, . . . , a„ 0, . . . , 0} with 0 in F and 
that correspondingly every quadratic form x'Ax is equivalent in F to 
0\Xi "j" • • • 4“ 

The problem of finding necessary and sufficient conditions for two quad¬ 
ratic forms with coefficients in a field F to be equivalent in F is one in¬ 
volving the nature of F in a fundamental way, and no simple solution of this 
problem exists for F an arbitrary field. In fact, we can obtain results only 
after rather complete specialization of F, and these results may be seen to 
vary as we change our assumptions on F. 

• We leave to the reader the explicit formulation of the definitions of equivalence in 
F of two forms, of two bilinear forms, and of two bilinear forms under cogredient linear 
mappings, where all the forms considered have elements in F. 
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The simplest conditions are those given in 

Theorem 12. Let F be a field with the property that for every a of F there 
exists a quantity b such that b 2 = a. Then two symmetric matrices with ele¬ 
ments in F are congruent in F if and only if they have the same rank . 

For every A = A' of rank r is congruent to A 0 = diag {ai, . . . , a r , 0, 

. . . , 0} for a» t* 0 and a< = h\ for hi j* 0 in F. Then if P = diag {6J" 1 , . . . , 
&r l , 0, . . . , 0}, PA 0 P' = diag {/ r , 0}. If also B = B’ has rank r, then B is 
also congruent in F to diag { I r , 0} and hence to A . The converse follows 
from Theorem 2.2. 

We then have the obvious.consequences. 

Corollary I. Two symmetric matrices whose elements are complex num¬ 
bers are congruent in the field of all complex numbers if and only if they have 
the same rank . 

Corollary II. Let F be the field of either Theorem 12 or Corollary I. 
Then two quadratic forms with coefficients in F are equivalent in F if and only 
if they have the same rank. Hence every such form of rank r is equivalent 
in F to 

(49) x 2 + . . . + x 2 . 

14. Addition of matrices. There is one other result on symmetric matrices 
over an arbitrary field which will be seen to have evident interest when we 
state it. Its proof involves the computation of the product of two matrices 



which have been partitioned into two-rowed square matrices whose ele¬ 
ments are rectangular matrices. If these matrices were one-rowed square 
matrices, that is to say in F } we should have the formula 


(51) 


/AiBi + AtBz AiBz + A2Bi\ 

\A 3 B 1 + A 4 B 3 AzBt + A 4 B 4 ) 


But it is also true that if the partitioning of any A and B is carried out so 
that the products in (51) have meaning and if we define the sum of two 
matrices appropriately, then (51) will still hold. Thus (51) will have major 
importance as a formula for representing matrix computations. 

We now let A = (a*,) and B = (6 t; ), where i = 1, . . . , m and j = 1, 
. . . , n, so that A and B are m by n matrices. Then we define 

(52) S = A + B = ( sh ), sn = an + 

(i - 1, . . . , m ;j = 1, . . . , n) . 
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We have thus defined addition for any two matrices of the same shape such 
that A + B is the matrix whose elements are the sums of correspondingly 
placed elements in A and B. 

The elements of our matrices are in a field F, and a w + &», = &»/ + a*/. 
If C = (c„) is also an m by n matrix, we have (a»/ + bn) + c< ; = a*, + 
(bij + dj). Hence we have the properties 

(53) A + B = B + A , (A + B) + C = A + (B + C) . 

We observe also that if 0 is the m by n zero matrix, then A + 0 = A, 
A + (-1 • A) = 0. 

Now let C = (Cjk) be any n by q matrix. Then 

(54) (A + B)C = AC + BC . 

For this equation is clearly a consequence of the corresponding property 
of F, that is, of 

(55) {an + bij)cjk = aijCjk + . 

We have thus seen that addition of matrices has the properties (53) al¬ 
ways assumed for addition of the elements of our matrices and that the 
law (54), which we call the distributive law for matrix addition and multipli¬ 
cation, also holds. Clearly if D is a matrix such that DA is defined, then 
similarly we have 

(56) D(A +B) - DA + DB . 

Observe, however, that if n > 1 and A and B are n-rowed square mat¬ 
rices, then |A + B\ and \A \ + |B| are usually not equal. For example, 
if A and B are equal, then |A| 4- |B| = 2|A| while | A + B | = 12 A | = 
2* | A |. 

Let us now apply our definitions to derive (51). We let A = (a*,-) be an 
m by n matrix, B = (&,•*) be an n by q matrix, and Ai be an s by t matrix, 
so that Ai = (a<,*) but now with i = 1 , . . . , s and j = 1 , . . . , t. Then 
(51) has meaning only if B\ has t rows, and we thus assume that B\ is a 
matrix of t rows and g columns. Our partitioning is now completely de¬ 
termined and necessarily A 2 has s rows and n — t columns, A 3 has m — s 
rows and t columns, A 4 has m — s rows and n — t columns, B 2 has t rows 
and q — g columns, B z has n — t rows and g columns, Bi has n — t rows 
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and q — g columns. The element in the ith row and fcth column of A B may 
clearly be expressed as 

n t n 

y £ a abjk — ^y^ aubjk + ^y^ Qijbjk 
i - l j - i ; -1 +1 

(i = 1, . . ., to; k = 1,. . . , q) . 

But this equation is equivalent to the matrix equation given by 

(57) A = (A A), B = (J) , AB = AA + , 

where we define A, A, A, A by 

(58) A = , A = (£), Ei = (5i, A), A = (A, A) . 


Moreover, we may obtain (51) from (57) by simply using the ranges 1, 2, 
. . . , s and s + 1, . . . , m for i separately, as well as 1, 2, . . . , g and 
g + 1, . . . , q for j separately. In matrix language we have used (58) and 
computed 


(59) 


(1) <«" *> - (J* IS) ■ (!>• - (IS as) 


as the result of partition of matrices and then have used (57) and addition 
of matrices in (59) to give (51). 

We shall now apply the process above to prove the following theorem on 
symmetric matrices mentioned above. 

Theorem 13. Let Ai and Bi be r -rowed nonsingular symmetric matrices 
with elements in F, and A and B be the corresponding n-rowed symmetric 
matrices 



of rank r. Then A and B are congruent in F if and only if Ai and Bi are 
congruent in F. 

For if A i and B\ are congruent in F there exists a nonsingular matrix P\ 
such that P\A X P[ = Si. Then P = diag {Pi, / n -r} is nonsingular, and com- 
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putation by the use of (51) gives PAP' = B. Conversely, if PAP' = B for 
a nonsingular matrix P we may write 



where Pi is an r-rowed square matrix, and we shall have 



PAP' = P 


a* 


_ (p,a x p[ 

o ) ~ V*ViiP', 


PiWa 

p^p'J 


But then B x = P X A X P' X) and P x must be nonsingular since \B X \ = 
|P(I * |^i | • I^Pil 9* 0. Hence B x and A x are congruent in F. 

The result above thus states that, if / and g are quadratic forms with 
coefficients in F so that / and g are equivalent only if they have the same 
rank r, then/ may be written as a nonsingular form / 0 in r variables, g may be 
written as a nonsingular form g 0 in r variables, and finally the original forms 
/ and g are equivalent in F if and only if / 0 and g 0 are equivalent in F. 


ORAL EXERCISES 


1. Compute A -f B if 



B = 



/-3 4 5\ /-3 5 3\ 

b) A = 5-1 2 , B = I 4-1 1 I 

\ 3 11/ \ 5 2 1/ 

2. Verify that ( A + BY = A' + B' { or any m by n matrices A and B. 

3. Show that every n-rowed square matrix A is expressible uniquely as the sum 
of a symmetric matrix B and a skew matrix C. Hint: Put A = B -f C with B = B' t 
— C = C', compute A', and solve. 


16. Real quadratic forms. We shall close our study of symmetric mat¬ 
rices and hence of quadratic forms with coefficients in F as well by consider¬ 
ing the case where F is the field of all real numbers. Let then / = x'Ax 
have rank r so that we may take / = a x x\ +...+, a r xj for real 5 * 0. 
We now call the number of positive a» the index t of / and prove 

Theorem 14. Two quadratic forms with real coefficients are equivalent in 
the field of all real numbers if and only if they have both the same rank and the 
same index . 
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For proof we observe, first, that by Section 14 there is no loss of generality 
if we assume that / = x'Ax where A is a nonsingular diagonal matrix, that 
is, r = n. Moreover, there is clearly no loss of generality if we take A = 
diag {di, . . . , d t , —d t + 1 , . . . , — d r \ for positive d*. Then there exist real 
numbers 9 * 0 such that p\ = d { , and if P is the nonsingular matrix diag 
{pF 1 ; • . . , P7 1 } then PAP ' is the matrix diag {1, . . . , 1, — 1, . . . , — 1} 
with t elements 1 and r — t elements — 1 . Thus, / is equivalent in F to 

(61) (x\ + . . . + x?) - (x ? +1 + . . . + x*) . 

Now, if g has the same rank and index as/, it is equivalent in F to (61) and 
hence to/. Conversely, let g have rank r = n and index s so that g is equiva¬ 
lent in F to 

(62) (x\ + . . . + x?) - (zj +1 + . . . + x?) . 

We propose to show that s = t. There is clearly no loss of generality if we 
assume that s ^ t and show that if s > t we arrive at a contradiction. 

Hence, let s > t. Our hypothesis that the form / defined by (61) and 
the form g defined by (62) are equivalent in the real field implies that there 
exist real numbers d i; such that if we substitute the linear forms 

(63) x% = diii/i • • • "h di n y n (i — 1, . . . , w) , 

in / = f(x u . . . , x n ), we obtain as a result (y\ + . . . + yl) - (t/J +1 + 

. . . + yl). Put z x = x 2 = . . . = x t = 0 and y a+l = . . . = y n = 0 in (63) 

and consider the resulting t equations 

(64) dnyi + . . . + di 9 y 9 = 0 (i = 1, . . . , t) . 

These are t linear homogeneous equations in 5 > t unknowns, and there 
exist real numbers v iy ... f v a not all zero and satisfying these equations. 
The remaining n — t equations of (63) then determine the values of x,- as 
certain numbers w, for j = t + 1 , . . . , n; and we have the result h = 

/(0, 0, . . . , 0 , u t+l , . . . , u n ) = rj + . . . + v* > 0. But clearly h = 

— (u ? +1 + . . . + ul) ^ 0 , a contradiction. 

We have now shown that two quadratic forms with real coefficients and 
the same rank are equivalent in the field of all complex numbers but that, 
if their indices are distinct, they are inequivalent in the field of all real 
numbers. We shall next study in some detail the important special case 
t = r of our discussion. 
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A symmetric matrix A and the corresponding quadratic form / = x'Ax 
are called semidefinite of rank r if A is congruent in F to a matrix 

fal r 0 \ 

\o o; * 

for a 0 in F. Thus / is semidefinite if it is equivalent to a form a(x\ + 

. . . + a?J). We call A and / definite if r = n, that is, A is both semidefinite 
and nonsingular. 

If F is the field of all real numbers, we may take a = 1 and call A and / 
positive or take a = — 1 and call A and / negative. Then A and / are nega¬ 
tive if and only if —A and —/ are positive. Thus we may and shall re¬ 
strict our attention to positive symmetric matrices and positive quadratic 
forms without loss of generality. 

If f(x i, . . . , x n ) is any real quadratic form, we have seen that there 
exists a nonsingular transformation (63) with real da such that f = y\ + 

• ••+!/?— ( 2/?+1 + • . . + yl). If Ci, ...» c n are any real numbers and 
if we put Xi = d in (63), there exist unique solutions t/, = d, of the result¬ 
ing system of linear equations, and the d, may readily be seen to be all zero 
if and only if the c, are all zero. Now if t < r, we have/ < 0 for y x = ... = 
y t = 0, y t +1 = l,/(ci, . . . , c n ) < 0. Conversely, if/(ci, . . . , c n ) < 0, then 
t < r. For otherwise / - y\ + . . . + 2/?, /(ci, . . . , c n ) = d\ + . . . + 
di ^ 0. If t = r < n, then we put y r +i = 1 and all other t/ ; = 0 and have 
Ci, . . . , c n not all zero such that /(ci, . . . , c n ) = 0. Hence, if /(ci, . . . , 
c n ) > 0 for all real Ci not all zero, the form / is positive definite. Conversely, 
if / is positive definite, we have / = y\ + . . . + yl, /(ci, ...,c*) = dj + 

. . . + di > 0 for all d, not all zero and hence for all c< not all zero. We 
have proved 

Theorem 16. A real quadratic form f(xi, . . . , x n ) is positive semidefinite 
if and only if f(ci, . . . , c B ) 0 for all real Ci, is positive definite if and only 
if f(ci, ...» Co) > 0 for all real c* not all zero . 

As a consequence of this result we shall prove 

Theorem 16. Every principal submatrix of a positive semidefinite matrix 
is positive semidefinite, every principal submatrix of a positive definite matrix 
is positive definite. 

For a principal submatrix B of a symmetric matrix A is defined as any 
m-rowed symmetric submatrix whose rows are in the iith, . . . , i m th rows 
of A and whose corresponding columns are in the corresponding columns 
of A. Put Xq = {Xi x , . . . , Xi m ) so that g — x' 0 Bx 0 is the quadratic form with 
B as matrix, and we obtain g from / by putting x, = 0 in / for j 7 * u. 
Clearly, if / ^ 0 for all x< = a, then g ^ 0 for all values of the x ik , and 
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hence B is positive semidefinite by Theorem 15. If A is positive definite 
and B is singular, then g = 0 for x lk not all zero, and hence / = 0 for the 
xj above all zero and for the Xi k not all zero, a contradiction. 

The converse of Theorem 16 is also true, and we refer the reader to the 
author’s Modern Higher Algebra for its proof. We shall use the result just 
obtained to prove 

Theorem 17. Let A be an m by n matrix of rank r and with real elements. 
Then AA' is a positive semidefinite real symmetric matrix of rank r. 

For we may write 



where P and Q are nonsingular matrices of m and n rows, respectively. 
Then QQ' is a positive definite symmetric matrix, and we partition QQ' so 
that Qi is an r-rowed principal submatrix of QQ'. By Theorem 16 the ma¬ 
trix Q i is positive definite, 


AA' = P 




= PCP' 


is congruent to the positive semidefinite matrix C of rank r and hence has 
the property of our theorem. 


EXERCISE 

What are the ranks and indices of the real symmetric matrices of Ex. 2, Sec¬ 
tion 11? 



CHAPTER IV 


LINEAR SPACES 

1. Linear spaces over a field. The set V n of all sequences 

(1) u = (Cl, . . . , c n ) 

may be thought of as a geometric n-dimensional space. We assume the 
laws of combination of such sequences of Section 1.8 and call u a point or 
vector , of F„. We suppose also that the quantities Ci are in a fixed field F 
and call c* the ith coordinate of u, the quantities Ci, . . . , c n the coordinates 
of u. The entire set V n will then be called the n -dimensional linear space 
over F. 

The properties of a field F may be easily seen to imply that 

(2) (u + v) + w = u + (y + w) , u + v = v + u 

(3) a(bu) = (i ab)u , (a + b)u = au + bu , a(u + v) = au + av 

for all a, 6 in F and all vectors u, v } w of V n . The vector which we designate 
by 0 is that vector all of whose coordinates are zero, and we have 

(4) u + 0 = u, u + (— u) = 0 , 
where —u = (— a, . . . , — c n ). Then 

(5) 0 • u = 0 , 1 • u = u , —u= — l • u . 

Note that the first 0 of (5) is the quantity 0 of F, and the second zero is the 
zero vector. We shall use the properties just noted somewhat later in an 
abstract definition of the mathematical concept of linear space, and we leave 
the verification of the properties (2), (3), (4), and (5) of V n to the reader. 

2. Linear subspaces. A subset L of V n is called a linear subspace of V n if 
au + bv is in L for every a and b of F f every u and v of L. Then it is clear 
that L contains all linear combinations 

(6) u = aiUi + . . . + a m u m , 
for a% in F , and u» in L. 
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We observe that the set of all linear combinations (6) of any finite num¬ 
ber m of given vectors, u h . . . , u m is a linear subspace L of V n according 
to this definition. If now L is so defined, we shall say that u h . . . , u m 
span the space L and shall write 

(7) L = [ut, , u m } . 

It is clear that, if e is the vector whose jth coordinate is unity and whose 
other coordinates are zero, then e if ... f e n span V n . 

The space spanned by the zero vector consists of the zero vector alone 
and may be called the zero space and designated by L = {0). In what fol¬ 
lows we shall restrict our attention to the nonzero subspaces L of V ni and, 
for the time being, we shall indicate, when we call L a linear space over F, 
that L is a linear subspace over F of some V n over F. 

3. Linear independence. Our definition of a linear space L which is a 
subspace of V n implies that every subspace L contains the zero vector. If 
L = {ui, . . . , u m }, then 0 = Oui + . . . + 0u m . Hence the zero vector of 
L is expressible in the form (6) in at least this one way. We shall say that 
Ui, . . . , Um are linearly independent in F or that u h ... } u m are a set of 
linearly independent vectors (of V n over F), if there is no other such expres¬ 
sion of 0 in the form (6). Thus, Ui, . . . , u m are linearly independent if it 
is true that a linear combination a^Ui + . . . + a^Um = 0 if and only if 
ai = a 2 = . . . = a m = 0. If . . . , u m are not linearly independent in F, 
we shall say that they are linearly dependent in F. 

A set of vectors u\, ... , Um are now seen to be linearly independent in F 
if and only if the expression of every u of L = {ui, . . . , u m } in the form ( 6 ) 
is unique. For this property clearly implies linear independence as the spe¬ 
cial case u = 0. Conversely, if u h . . . , u m are linearly independent and 

u = a\U\ + . . . + a m Um.= Mi + . - . + b m ii m , then 0 = («i — &i )u x + 

. . . + (a m — b m )u m , a* — bi = 0, a» = as desired. We now make the 

Definition. Let L = {ui,..., u ro } over F and Ui,..., u m be linearly inde¬ 
pendent in F. Then we shall call Ui, . . . , u m a basis over F of L and indicate 
this by writing 

(8) L = UiF + . . . + u m F . 

It is evident that V n = eiF + . . . + e n F. But we may actually show 

that every subspace L spanned by a finite number of vectors of V n has a 
basis in the above sense. Observe, first, that the definition of linear inde- 
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pendence in case m = 1 is that au = 0 only if a = 0 and thus that u^O. 
Then let Ui, . . . , u r be linearly independent vectors of V n and m ^ 0 be 
another vector. Then either iti, . . . , u r , m are linearly independent in F or 
OLiUi + . . . + drUr + a Q u = 0 for a* not all zero. If a 0 = 0, then aiUi + 

. . . + a T u r = 0, from which a x = . . . = a r = 0, a contradiction. But then 
u = (-ao 1 ^)^ + . . . + ( — ai‘ 1 a n )w r is in {u ly . . . , w r ). It follows that, 
if Mi, . . . , u m are any m distinct nonzero vectors, we may choose some 
largest number r of vectors in this set which are linearly independent in F 
and we will then have the property that all remaining vectors in the set are 
linear combinations with coefficients in F of these r. We state this result as 
Theorem 1. Let L be spanned by m distinct nonzero vectors Ui, . . . , u m . 
Then L has a basis consisting of certain r of these vectors , 1 < r < m. 


EXERCISES 

1. Determine which of the following sets of three vectors form a basis of the 
subspace they span. Hint: It is easy to see whether some two of the vectors, say 
ui and m 2 are linearly independent. To see if u i} u 2 , m 3 are linearly dependent we 
write its = xui + yu 2 and solve for x and y. 


а) Mi — (1, 3, —3) , 

б ) mi - ( 3 , 1 , 2 ) , 

c) Mi — (-1, -1,1), 

d) Mi - ( 1 , - 2 , - 1 , 3 ), 

e) mi = (1, 1, 1, -1), 

/) Mi - (1, -1, 1, -1), 


u 2 = (2, 5, -2) , 
u 2 = (4, 1, 3) , 
m 2 = (3, 2, 1 ) , 
u 2 = ( 2 , - 1 , 1 , 6 ), 
m 2 — (2, 2, 2, —2) , 
M2 - (1, 1, 1,0), 


Us = (1, 1, 6) 

m 3 = (-1, -1,0) 

M3 - (7, 3, 9) 
hi - (0, -1, -3, 1) 
m 3 = (1,2, 3, 4) 
u 3 — (o, 1, 5, 3) 


2. Show that the space L = { ui, m 2 , m 3 ) spanned by the following sets of vectors 

is V ». Hint: In every case one of the vectors, say mi, has first coordinate not zero, 
and L contains U 2 — x 2 Ui = (0, b 2t c 2 ), u 3 — x 3 ui = (0, 6 3 , c 3 ). Some linear combina¬ 
tion of these two quantities has the form (0, 0, d) and L contains e 3 = (0, 0, 1). It 
is easy to show then that L contains e 2 = (0,1, 0) and ei = (1, 0, 0), L = V 3 . 

а) Mi — (1, 2, 3) , Ui = (2, 3, 1) , m 3 = ( 1, 3, 2) 

б) mi = (0, 1, 8) , m 2 = (1, -3, 6) , m 3 = (1, -1, 23) 

c) mi = (0, 3, 2) , Ui = (0, 2, 1) , m 3 = (1, 5, 4) 

3. Determine whether or not the spaces spanned by the following sets of vectors 
Mi, ut coincide with those spanned by the corresponding v x , v 2 . Hint: If L\ = 
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(fci, v*} = L 2 = {t>i, t> 2 }, there must exist solutions xi, x 2 , x 3 , x 4 of the equations 
Vi = XiWi + XiU 2 t v 2 = x 3 wi + z 4 u 2 such that the determinant x t x 4 — xax 3 7 ^ 0 . 

a) u»- (3, -1,2,1), u 2 = (1,4, 6,1), 

ri= (7, -11, -6,1), r 2 = (7,2,10,3) 

fe) lit - (1,2, -1,2), u 2 = (1,2, 3, 4), 

(1,2, -13, -4), *- (0,0,4, 1 ) 

c) Ui - ( 1 , - 1 , 2, —3) , u 2 = (2, -2, 4, - 6 ) , 

t>i = (3, — 3, 6 , —9) , v 2 = (4, —4, 8 , —9) 

d) m- (1, -1,0,3), u 2 = (2,1, 0,1), 

*1 = (-1, “2, 1,2), *- (2, -2,0,6) 

e) m-(1,-2,0,1), in - (2, -1,1,0), 

t’i = (3, -3,1,1), ** (-1, -1, -1,1) 

/) m - (1,0, 1, —2) , tn - (0, 1, -1,2), 

— (2, -2,4, -8), V 2 — (1, 1,0,0) 

4. The row and column spaces of a matrix. We shall obtain the principal 
theorems on the elementary properties of linear spaces by connecting this 
theory with certain properties of matrices which we have already derived. 
Let us consider a set of m vectors, 

(9) Ui = (an, . . . , ciin) (i = 1, . . . , m) , 

of V n over F. Then we may regard u, as being the ith row of the correspond¬ 
ing rn by n matrix A — (a» ; ) and the space L = [u iy ... } Um} as being 
what we shall call the row space of A. Thus every m by n matrix defines a 
linear subspace of F n > every subspace of V n spanned by m vectors defines 
a corresponding m by n matrix. 

If P = ( p k i ) is any q by m matrix, the product PA is a q by n matrix. 
The jth coordinate of the vector 

(10) Wk = PkiUi + . . . + PkmUm 

is pk\dij + . . . + PkmUmj, that is, the element in the kth row and jth col¬ 
umn of PA. Hence the k th row of PA is that linear combination of the rows 
of A whose coefficients form the kth row of P. 

We have now shown that every row of PA is in the row space of A. It 
follows that the row space of PA is contained in that of A. If P is nonsingu- 
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lar, then the result just derived implies that the row space of A = P~ l (PA) 
is contained in the row space of PA and therefore that these two linear 
spaces are the same. Thus we have 

Lemma 1. Let P be an m -rowed nonsingular matrix. Then the row spaces 
of PA and A coincide. 

In the proof of Theorem 3.6 we used the matrix product equivalent of 
the elementary transformation Theorem 2.4, and we shall now state this 
theorem as the useful 

Lemma 2. Let Abe an m by n matrix of rank r. Then there exist nonsingular 
matrices P and Qo/m and n rows , respectively , such that 

(11) PA = , AQ - (H 0) , 

where G and H have rank r, G is an r by n matrix , H is an m by r matrix. 

We use (11) and note that the rows of G differ from those of PA only in 
zero rows. Then the rows of G span the row space of PA. By Lemma 1 
we have 

Lemma 3. The row spaces of G and A coincide. 

We shall use this result in the proof of 

Theorem 2. The r rows of G form a basis of the row space of A. Any 
r + 1 vectors of the row space of A are linearly dependent in F. 

For we may designate the rows of G by Vi, . . . , v r . Our definition of G 
implies that there is no loss of generality if we permute its rows in any de¬ 
sired fashion. Thus, if b&i + ... + b r v r = 0 for 6; in F not all zero, we 
may assume for convenience that hi 0. The determinant of the r-rowed 
square matrix 

(12) P = (£) , P l = (b h ...,b r ), Pi = (0, I r -l) 

is clearly b iy and hence P is nonsingular. Then PG is r-rowed and of rank r. 
But this is impossible since b\Vi + . . . + b r v r = 0 is the first row of PG. 
Hence the rows of G are linearly independent in F y they span the row space 
of G and of A , they are a basis of the row space of A. 

Assume, now, that w k = b k \V\ + . . . + b kr v T for k = 1, . . . , r + 1, so 
that w k are any r + 1 vectors of the row space L = v\F ■+...+ v r F of A. 
Define B = (b k %) and obtain w k as the kth row of the r + 1 by n matrix 
BG. By Theorem 3.6 the rank of BG is at most r, and by Lemma 2 there 
exists a nonsingular (r + l)-rowed matrix D = ( d gk ) for g y k = 1, . . . , 
r + 1, such that the (r + l)st row of D{BG) is the zero vector. But then 
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d r + 1 , 1 W 1 + ... + rfr+ 1 , r+ittV+i = 0, the d r +i, * form a row of the nonsingu¬ 
lar matrix D and cannot all be zero, the vectors w i} . . . , w r+ 1 are linearly 
dependent in F. This proves Theorem 2. 

We now apply Theorem 1 to obtain 

Theorem 3. The matrix G of Lemma 2 may he taken to he a suhmatrix of A. 
The integer r of Theorem 1 is in fact the rank of the m by n matrix whose i th 
row is Ui. 

For let A be an m by n matrix whose ith row is Ui and let r be the integer 
of Theorem 1, the rank of A be s. After a permutation of the rows of A, if 
necessary, we may assume that u h . . . , u r are a basis of the row space of A. 
Then u k = bk\Ui + . . . + h kr u r for b k i in F, and k = r + 1, . . . , m. The 
matrix P given by 

(,3) Mb' ?.-.)• 

is nonsingular, and it is clear that if 


(14) 


G - 


Ul 


n r 


then PA is given by (11). It follows that s is the rank of G, s < r. li s < r 
we apply Lemma 2 to obtain a nonsingular r-rowed matrix D = ( d gk ) such 
that the rth row of DG = 0, the rth row of D is not zero, d n ui + . . . + 
d rr Ur = 0 contrary to our hypothesis that ui, . . . , u r are linearly inde¬ 
pendent. This completes our proof. 

We may now obtain the principal result on linear subspaces of V n . 

Theorem 4. Every linear subspace L over F of V n over F has a basis. Any 
two bases of L have the same number r < n of vectors , and we shall call r the 
order of L over F. 

For V n — Ci F + . • . + e n F, and by Theorem 2 any n + 1 vectors of 
V n are linearly dependent in F. Thus any linear subspace L over F of V n 
contains at most n linearly independent vectors. It follows that there 
exists a maximum number r < n of linearly independent vectors u if . . . , u r 
in L, and that u , Wi, . . . , u r are linearly dependent in F for every u of L. 
By the proof of Theorem 1 the vector u is in \u lf . . . , u r }, L = uiF + 
. . . + u r F. But if also L = v x F + ... + v t F , then Theorem 2 implies that 
s < r and similarly that r < s, r - s is unique. 
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In closing this section we note that the rows of the transpose A' oi A 
are uniquely determined by the columns of A and are, indeed, their trans¬ 
poses. Thus we shall call the row space of A' the column space of A. It is 
a linear subspace of V m . We call the order of the row and column spaces 
of A, respectively, the row and column ranks of A. By Theorem 3 the rank 
of A is its row rank. Also A and A' have the same rank and we have proved 
Theorem 5. The row and column ranks of a matrix are equal to its rank. 


EXERCISES 

1 . Solve Ex. 1 and 2 of Section 3 by the use of elementary transformations to 
compute the rank of the matrices whose rows are the given vectors. 

2 . Form the four-rowed matrices whose rows are the vectors wj, ui , v if v 2 of Ex. 3, 
Section 3. Show thus that L\ = {iq, 1 / 2 }= L* = {t>i, t’ 2 } if and only if the ranks 
of the corresponding matrices A are equal to the order of the subspace L h that is, 
the rank of the matrices formed from the first two rows of each A. 

3. Find a basis of the row space of each of the following matrices, the basis to 
consist actually of rows of the corresponding matrix. 


a) 


1-2 10 
2-4 2 0 

2-3-1 1 

4-7 11 

2-3-1 1 


b) 


2 3 0 1 

-2 -1 -1 1 
-2-10 3 

0 4-12 

2 12-2 


c) 


-10 0 1 
0-11 1 
2 3 0 -1 

5 10 0 

12 3 4 


d) 


12-1 3 4 

3 5 1 0 2 

11 3-6 -6 

4 7 0 3 6 

10 7-15-16 


4. Let A be a rectangular matrix of the form 


( 


A 

A 


1 

3 



where Ai is nonsingular. Show that then there exists a nonsingular matrix P such 
that 



Give also a simple form for P. Hint: Show that the rows of At are in the row space 
of A,. 
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5. Let the matrix of Ex. 4 be either symmetric or skew. Show that the choice 
of P then implies that 



Show also that, if the order of A\ is the rank of A, the matrix A* = 0. 

5. The concept of equivalence. In discussing the properties of mathemati¬ 
cal systems such as fields and linear spaces over F it becomes desirable 
quite frequently to identify in some fashion those systems behaving exactly 
the same with respect to the given set of definitive properties under con¬ 
sideration. We shall call such systems equivalent and shall now proceed to 
define this concept in terms of that of junction. 

Let G and G' be two systems and define a single-valued function f on G 
to G'. In elementary mathematics it is customary to call G the range of 
the independent variable and G' the range of the dependent variable. How¬ 
ever, it is more convenient in general situations to say that / is a function 
on G to G' or that / is a correspondence on G to G'. Then/ is given by 

(15) g -> g’ = Kg) 

(read g corresponds to g') such that every element g of G determines a unique 
corresponding element g' of G'. In elementary algebra and analysis the sys¬ 
tems G and G' are usually taken to be the field of all real or all complex 
numbers and (15) is then given by a formula y = f(x). But the basic idea 
there is that given above of a correspondence (15) on G to G'. This concept 
may be seen to be sufficiently general as to permit its extension in many 
directions. 

Suppose now that (15) is a correspondence such that every element g f of 
G' is the corresponding element f(g) of one and only one g of G. Then we 
call (15) a one-to-one correspondence on G to G'. It is clear that (15) then 
defines a second one-to-one correspondence 

(16) fig) = g'-* 9 > 

which is now on G' to G, and we may thus call (15) a one-to-one correspond¬ 
ence between G and G' and indicate this by writing 

(17) 

Note, however, that, if G and G' are the same system, the functions (15) 
and (16) are, in general, distinct. Thus we may let G be the field of all real 
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numbers and (15) be the function x —> 2x , so that (16) is the function 
x —> \x. Of course, if G and G ' are distinct it is not particularly important 
whether we use (15) or (16) to define our correspondence. 

We proceed to use the concept just given in constructing the fundamental 
definition of this section, that of equivalence. Let us consider two mathe¬ 
matical systems <?, G' of the same kind such as two fields or two linear 
spaces over a fixed field F. These systems consist of sets of elements g , A, 

. . . closed with respect to certain operations. (Thus, for example, we might 
have g + h in G for every g and h in G and also g' + h' in G' for every g' 
and h f in O'.) We then call G and G' equivalent if there exists a one-to-one 
correspondence between them which is preserved under the operations of 
their definition. We now see that we have defined two fields F and F* to 
be equivalent if there exists a one-to-one correspondence (15) between them 
such that ( g + h)' = g' + h ( gh )' = g'h' for every g and h of F. Let us 
then pass to the second case which we require for our further discussion of 
linear spaces. 

Let F be a fixed field and V consist of a set of elements such that u + v 
and au are unique elements of V for every u and v of V and a of F. Then 
we shall call V a general linear space over F. If V 0 is a second* such space 

and there is a one-to-one correspondence u *-> w 0 between V ami F 0 such 

that 

(u + v) 0 = u 0 + v 0 , (aw)o = au 0 

for every u and v of V and a of F, then we shall say that V and V 0 are equiva¬ 
lent over F. We have thus introduced two instances of what is a very im¬ 
portant concept in all algebra. 

The reader should observe that under our definition every mathematical 
system G is equivalent to itself and that if G is equivalent to a system G' f 
then G' is equivalent to G. Finally, if G ' is equivalent to (7", then G is 
equivalent to C". 


EXERCISES 

1 . Verify the statement that the field of all rational functions with rational coeffi¬ 
cients of the complex number ^2 is equivalent to the field of all matrices 



for rational a and b under the correspondence A <-> a + 6^2. 

* We use the notation Vo instead of V' to avoid confusion in the consequent usage 
of v! for the arbitrary vector of V' as well as for the transpose of u. 
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2. Verify the statement that the field of complex numbers is equivalent to the 
field of matrices 

C "«) 

for real a and b. 

3. Verify the statement that the set of all scalar matrices with elements in a field 
F forms a field equivalent to F. 

4. Verify the statement that the mathematical system consisting of all two-rowed 
square matrices with elements in F and defined with respect to addition and multi¬ 
plication is equivalent to the set of all matrices 


All 

0 

an 

0 ' 

0 

All 

0 

<*12 

an 

0 

<*22 

0 

k 0 

&21 

0 

<*22, 


under the correspondence indicated by the notation. 


(<*»/ in F) 


6. Linear spaces of finite order. We shall restrict all further study of 
linear spaces to linear spaces of order n over F. Then we define V to be a 
linear space of order n over a given field F if V is equivalent over F to V n 
over F. Clearly, our definition implies that every two linear spaces of the 
same order n over F are equivalent over F. 

Our definition also implies that, if V is a linear space of order n over F, 
then the properties (2), (3), (4), and (5) hold for every w, v , w , of V and 
a, b in F. Moreover, every quantity of V is uniquely expressible in the form 

(18) ci^i + . . . + c n v n , 
where the equivalence between V and V n is given by 

(19) CiVy + . . . + c n v n < -> (ci, . . . , c n ) . 


Conversely, define an = ua and u + v to be in V for every u, v of V and 
a of F. Then it may be shown very easily that, if (2), (3), (4), and (5) hold 
and every u of V is uniquely expressible in the form (18) for c» in F, then V 
is equivalent over F to F n . However, we prefer instead to define V by its 
equivalence to V n . This preference then requires the (somewhat trivial) 
proof of 

Theorem 6. Every linear subspace L of order r over F of V n is equivalent 
over F to V r . 

Thus we justify the use of the term linear subspace L of order r over F by 
proving that L is indeed what we have called a linear space of order r over F 
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contained in the space V n . For proof we merely observe that every vector, 
of L = UiF + . . . + u r F is uniquely expressible in the form 

(20) U = CiUi + . . . + CrU r 

for d in F. Thus u in L uniquely determines the c» in F and conversely. 
It follows that 

(21) u (ci, . . . , Cr) 

is a one-to-one correspondence between V and V T and it is trivial to verify 
that it defines an equivalence of L and V r . 

We have now seen that every linear space L of order n over F may be 
regarded as a linear subspace of a space M of order m over F for any 
m > n. Moreover , L = M if and only if m — it. 

Theorem 3 should now be interpreted for arbitrary linear spaces of order 
n, and we have a result which we state as 

Theorem 7. Let L = UiF + . . . + u n F and Vi, . . . , v m be in L so that 
there exist quantities a,*/ in F for which 


Vi = aiiUi + . . . + a in u n , 

and the coefficient matrix A = (aij) is defined. Then the number of the Vk 
which are linearly independent in F is the rank of the matrix A. Moreover, 
Vi, . . . , v m form a basis of L over F if and only if m = n and A is nonsingular. 

EXERCISES 

1. Verify the statement that the following sets of matrices are linear spaces of 
finite order over F and find a basis for each. 

a) The set of all rn by n matrices with elements in F. 

b) The set of all m by n matrices whose elements not in the first row are zero. 

c) The set of all n-rowed scalar matrices. 

d) The set of all n-rowed diagonal matrices. 

2. Find bases for the following linear spaces of polynomials with coefficients in F. 

a) All polynomials in x of degree at most three. 

b) All polynomials in independent variables x and y and degree at most two 
(in x and y together). 

c) All polynomials in x = t 2 + t and y = l z + t 2 and degree at most two in 
x and y. 

d) All polynomials in ^2, with F the field of all rational numbers. 

e) All polynomials in a primitive cube root of unity with F the field of all 
rational numbers. 
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/) All polynomials in w = u(i + 1) with F the field of all rational numbers 
and i 2 = — 1, u 2 = 2. Hint: Prove that 1, u, i, ui are a basis. 
g) The polynomials of (/) but with F the field of all real numbers. 

3. Let A = diag {1, —1, 2}. Show that /, A, A 2 form a basis of the set of all 
three-rowed diagonal matrices. 

4. Show that I, A, B, AB form a basis of the set of all two-rowed square mat¬ 
rices if 





5. Show that, if 


0 0 0 
0-10 
0 0 1 , 


/° 1 °\ 

B= 0 0 1 , 

\1 0 0 / 


then I, A, A 2 , B, AB, A 2 B, B 2 , AB 2 , A 2 B 2 form a basis of the set of all three-rowed 
square matrices. 


7. Addition of linear subspaces. If Li = {vi, . . . , v m ] and L 2 = 
{u>i, . . . , w q \ are linear subspaces over F of a space L of order n over F, 
the subspace L 0 = {Vi, . . . , v m , W\, .,. . , w q ] of L will be called the sum of 
Li and L 2 and will be designated generally by 

(22) L 0 = \L h U] . 

If the only vector which is in both L x and L 2 is the zero vector, we shall say 
that Li and L 2 are complementary subspaces of their sum and write 

(23) L 0 = L\ + Li . 

In this case the order of L 0 is the su^ of the orders of L\ and L 2 and in fact 
we shall show that, if L\ = v x F v m F } L 2 = w\F + . . . + w q F , 

then L 0 = v x F + . . . + v m F + w x F + . . . + w q F. 

For it is clear that a<V\ + . . . + OmV m + biWi + . . . + b q w q = 0 if and 
only if v = aiv x + . . . + a m v m = ( — bjwi + . . . + (— b q )w q is in both L x 
and L a . Thus v 0 implies that the a< and h, are not all zero and therefore 
that the vectors v < and Wj spanning L 0 are linearly dependent in F and do 
not form a basis of L 0 . Conversely, if necessarily v = 0, then the Vi and Wj 
do form a basis of L 0 , and L 0 has order m + q. 

If Li and L 0 are linear subspaces of L and if L 0 contains L h we may ask 
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whether a linear subspace L 2 of L 0 exists such that L 0 = Li + L 2 . The 
existence of such a space is clearly a corollary of 

Theorem 8. Let Li of order m over F be a linear subspace of L of order n 
over F. Then there exists a linear subspace L 2 of L such that L\ and L 2 are 
complementary subspaces of L. 

The result above may be proved by the method we used to prove Theo¬ 
rem 1 , where we apply this method to the set v Xj . . . , v m , u x , . . . , u n , in 
L = u\F + . . . + u n F. However, let us give, instead, a proof using matrix 
theory. We put L = F n , let G be the m by n matrix whose rows are the basal 
vectors v Xy ... t v m of L\. Then G has rank m , and there exists a nonsingular 
matrix Q such that the columns of GQ are a permutation of those of G, 

GQ = (Gx G*), 

where G x is nonsingular. Then the matrix 



is nonsingular, and A = A 0 Q~ l is obtained by permuting the columns of A 0 . 
But then 

a -(h) 

is nonsingular, the rows of A span F n , and the rows of H span the space L 2 
which we have been seeking. Moreover, it is clear that the rows of // are 
certain of the vectors which we defined in Section 2. 


EXERCISES 


1. Let L\ be the row space of each of the m by n matrices of Section 4, Ex. 3. Use 
the method above to find a basis of the corresponding V\» consisting of a basis of 
L\ and of a complementary space L 2 . 

2. Let the following vectors Ui span L\ 1 span L 2 . Find a complement in { L Xf L 2 \ 
to Li and to L 2 . 


ju t = ( 1 , - 1 , 1 , 1 ), 
a) l »1 = (1, 2, 1, 0) , 



u, = (1, 2, 0, 0), 
t-, = (4, 3, 1, 1), 

= ( 1 , 0 , 2 , - 1 ), 

», = (3, -2, 2,-1), 


— ( 2 , — 2 , 1 , 2 ) , 

t> 2 = (4, -1,4,0) 

ui = ( 1 , - 1 , 1 , 0 ), 

= (4, -3, 3, 2), 

( 0 , 1 , 2 , - 1 ), 

v t = (-5, 3, -4, 2) , 


= (1, -1,0, 1) 


M» = (1, 0, 0, 1) 

v 3 = (1, 4, 2, -3) 

u 3 = (2, 1, (>, -3) 
v> = (-2, 1, -2, 1) 
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8. Systems of linear equations. The set of all linear forms in xi, . . . , x n 
with coefficients in F is a linear space 

(24) L = X\F + . . . + x n F 
of order n over F . The left members 

(25) fi = /*(xi, . . . , X n ) = a»iXi + . . . + a in x n 
of a system 

(26) ^il^l “1“ • • • “1“ O'inXn = C{ (z = lj . . . , wi) , 

of in linear equations in n unknowns, are such forms and are in L. Then, 
if r is the rank of the m by n matrix A = (a,,) of coefficients of (26), we 
set 4 by Theorem 7 that eertnin r of the forms/, are linearly independent in 
F and the remaining m — r forms are linear combinations of these r. 

We may assume without loss of generality that the equations (26) have 
been labeled so that/i, . . . , f r are linearly independent in F , and 

(27) fk = bkifi + . . . + r/r (k = r + 1, . . . , m) 


for b k] in F. Then 


Ml 


( 28 ) A = 


\ U m 


(i = 1, . . . , m) , 


where u t , . . . , u r are linearly independent in F, and 

(29) Uk = bnui + . . . + b k rUr (fc = r + 1, . . . , vi ). 

If the system (2(j) is consistent, there exist quantities d h . . . , d n in F such 
that/i(rfi, . . . , d n ) = Ci. Then (27) implies that 

Ck = fk(d t , . . . , dn) = bklCi + • . . + bkrCr 

(k = r + 1 , . . . , m) . 


(30) 
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Define the augmented matrix A* of the system (26) to be the m by n + 1 
matrix 


(31) 


A* = 


t*i 


< 


W? = (u»> Ci) = (a iu . . . , a,*, Ci) 

(i = 1, . . . , m), 


and see that (29) and (30) imply that 

(32) uj = b kl uZ + . . . + bkrut (k = r + 1, . . . , m) , 


so that the rank of A* is at most r. But A* has A as a submatrix; A* has 
rank r. 

Conversely, if A* has the same rank r as A and we choose Ui, . . . , u r to 
be linearly independent, then uf, . . . , u* are clearly linearly independent. 
We then have (29) and may apply elementary row transformations of type 
1 to A* which add — ( b kl u\ + . . . + b kr u*) to u\ for k = r + 1, . . . , m. 
These replace the submatrix A of A * by 


(33) 



/ 


G = 


Ui 


Ur 


and replace A* by 


(34) 



Ci 


C x = 


\ ^ 


where the (m — r)-rowed and one-columned matrix has c* 0 = c k — 
(bk\Ci + . . . + bkrCr) as its elements. But, clearly, if any c* 0 ^ 0, the 
matrix A* has a nonzero (r + l)-rowed minor. This is impossible if A* is 
of rank r. 

We now see that the system (26) is consistent if and only if A* has the 
same rank r as A. Moreover, we have already shown that, if A * does have 
the same rank as A, then m — r of the equations (26) may be regarded as 
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linear combinations of the remaining r equations and are satisfied by the 
solutions of these equations. It thus remain^ to show that any system of r 
equations in Xi, . . . , x n with matrix of rank r has solutions. We may write 
x = (xi, . . v, x n ) and see that such a system may be regarded as a matrix 
equation 

(35) Gx' = v ’, v = (ci, . . . , c r ) . 

Before solving (35) we prove 

Lemma 4. Let G be an r by n matrix of rank r. TTien 

(36) G = (I r 0)Q 
for a nonsingular n-rowed matrix Q. 

For Lemma 2 states that G = (H 0)Qi, where Qi is n-rowed and non¬ 
singular and H is an r by r matrix of rank r. Then H is nonsingular, so is 
Q 2 = diag {//, / n -r}, (Ir 0)Q a = (H 0). Thus we have (36) for Q = 
QzQi* 

The system (35) may now be written as 

(37) (I r 0)y' = v', y = xQ f = (y u . . . , y n ) . 

But then y — xQ' is a nonsingular linear transformation and the are 
linearly independent linear forms in xi, . . . , x n . Evidently, (37) has the 
solution yi = c< for i = 1, . . . , r; the solution of (35) is then given by 
x = y(Q')~ l for yi = Ci(i = 1 , . . . , r) and y T + y . . . , y n arbitrary. Ob¬ 
serve that, if we choose the notation of the x* so that G = (C?i, G 2 ) with 
G\ an r-rowed nonsingular matrix, then 

(38) Q = (J ) , y = (Y u Yt) , x = (X u Xi) , 
where l r i and A'i have r columns. From this we obtain 

Q' = (JS J , X Q' = CM + X,G'i, X t ) , 


so that 


X 2 — Fa — (X r +i> • • • ; ^n) 7 

= -XjCrJ T" = (Cf» • • • t *V) f 
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and our solution of (35) is given by 

(39) = (I\ - X&)(G[)-' . 

But then we have solved for xi, . . . , x r as linear functions of x r +\, . . . , x n . 
For exercises on linear equations we refer the reader to the First Course in 
the Theory of Equations. 

9. Linear mappings and linear transformations. The system of equations 
(3.1) of a linear mapping was expressed in (3.32) as a matrix equation 
x f = y'A'. Let us interchange the roles of m and n, A and A 1 in this equa¬ 
tion. Then we see that a linear mapping may be expressed as a matrix 
equation 

(40) v = uA , 
where A is an m by n matrix and 

(41) u = (j/i, . . . ,y m ) , v = (i,, . . . , x n ) . 

Clearly, u is a vector of V m over F, v is a vector of V n over F, and (40) may 
be regarded as a correspondence u —> uA defined by A whereby every u 
of V m determines a unique vector uA of V n . We now proceed to formulate 
this concept more abstractly. 

Let L and M be linear spaces of respective orders m and n over F and 
consider a correspondence on L to M . Designate the correspondence by the 
symbol S y so that S is the function 

S : u —> u& 


(read u goes to u upper S) wherein every vector u of L determines a unique 
u s in M. Suppose also that 

(42) (au + &U 0) 5 = au s + bv$ 


for every a and b of F f u and w 0 of L. Then we shall call S a linear mapping 
of the space L on the space M and describe (42) as the property that S is 
linear. 

Suppose now that L = u\F + . . . + u m F and M = v\F + . . . + v n F y 
so that we are given not only the spaces L and M but fixed bases as well. 
Then a linear mapping S uniquely determines uf in M y and hence 


(43) 


uf = duV, + . . . + a in v n 


(i = 1, . . . , m) , 
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for an in F. Thus S determines also an m by n matrix A = (a i; ). But, con¬ 
versely, if A is an m by n matrix and we define u$ by (43) for the given ele¬ 
ments cii] of Ay then the property that S is linear uniquely determines the 
mapping S. This is true since every u of L is uniquely expressible in the 
form 

(44) u = yiUi + ... + y m u m , 
for 2 /< in F and (42) implies that 

( 45 ) u s = + . . . + y m ui . 

It follows that to every linear mapping S of L on M and given bases of L 
and My there corresponds a unique m by n matrix A and conversely. We 
shall call A the matrix of S with respect to the given bases of L and M. 

We now observe that 


m m / n \ n 

u s = S^y iU f =^PyA ) = 2 x ' t '' » 

t-i »= i i / ; = i 

where 

TO 

x,- = ^ a (j = 1, . . . , n). 

i * 1 

But then, if we assume temporarily that L = V m and M = V n and put 
v = u s in (40) and (41), we see that S is the linear mapping 

(46) u—>u s = uA 

for the given matrix A. Thus every linear mapping which is a change of 
variable as in (3.1) may be regarded as a linear mapping of the space V m 
on the space V n . 

Let us next observe the effect on the matrix defined by a linear mapping 
of a change of bases of the linear spaces. Define new bases of L and M, 
respectively, by 


4 0) =^p* , 

t- l 


”i (0) =]£?<,«',» 

y-x 


(47) 
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for k = 1, . . . , m and l = 1, . . . , n. Then P = ( p ki ) and Q = (qn) are 
nonsingular and, as we saw in (3.33), we may also write the second set of 
equations of (47) in the form 

n 

(48) ", 

2-1 

where R = (r,i) = Q”" 1 . We apply the linearity of S to (47) and obtain 

m 

(ul 0) ) s = Substituting (43) and (48), we have 

»« i 

n 

(49) (i4 0) ) s -£M 0) > 

2 » 1 

where 6±/ = Hence the matrix B = (6a-/) of the linear mapping 

S with respect to our new bases is given by 

(50) B = PAQr 1 . 

Since P and Q -1 are arbitrary nonsingular matrices of m and n rows, re¬ 
spectively, we see that changes of basis in L and M replace A by an equiva¬ 
lent matrix. Thus any two equivalent m by n matrices define the same 
mapping of L on M. 

If L and M are the same space we shall henceforth call a linear mapping 
S of L on L a linear transformation of L. Since we are now considering only 
a single space, the only possible meaning of the u, and Vj in (43) can be that 
of a fixed basis u lf . . . , u n of L of order n over F and of a second basis 
v if . . . , v n of L. Let us restrict our attention to the case where Vi = iq. 
Then we define the matrix A of a linear transformation S on L with respect 
to a fixed basis ui, . . . , u n of L to be the matrix A determined by (43) 
with Ui = We have defined thereby a one-to-one correspondence be¬ 
tween the set of all n-rowed square matrices with elements in F and the set 
of all linear transformations on L of order n over F. If L = V n , such a 
correspondence is given by 

u—>u s = uA. 

Clearly, we should and do call S a nonsingular linear transformation if A is 
nonsingular. Since we may then solve for u = u 5 A -1 , we see that a non¬ 
singular linear transformation defines a one-to-one correspondence of L 
and itself. 
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We now observe the effect on A of a change of basis of L. We use (47) 
but note that u 4 = v { and u[ a) = i'[ 0) , so that we have P — Q. Hence a 
change of basis* of L with matrix P replaces the matrix A of a linear trans¬ 
formation by 

(51) B = PAP~ X . 


It is now clear that in order to study those properties of a linear trans¬ 
formation on L which do not depend on the basis of L over F we need only 
study those properties of square matrices A which are unchanged when we 
pass to PAP~ l . We shall call two matrices A and B similar if (51) holds 
for a nonsingular matrix P and shall obtain necessary and sufficient condi¬ 
tions that A and B be similar in our next chapter. 


EXERCISES 


1 . Let S be the linear mapping (46) of F 3 on F« defined for the following mat¬ 
rices A. Find the vectors of V, into which ( — 2, 3, 4), ( 1 , 0, 0), (0, 1 , 0) of F a are 
mapped by S. 


/ 1 -3 0 1\ 

a) A = ( 2-6-1 2 ) 

\-l 3 1 -1/ 


/—2 3 4 0\ 

b) A = ( -l 2 0 4) 

V 0 0 2 -3/ 


/2 -1 2 — 5 \ 

C )A= 0-2 3-2 

\1 1 -1 - 1 / 


/—2 0 0 0 \ 
d) A = I 0 -1 0 0) 
\-l 1 1 1 / 


2. Show that the linear transformations (46) of F s defined for the following mat¬ 
rices .4 are nonsingular and find their inverse transformations. Apply both S and 
iS ’ -1 to the vectors of Ex. 1. 


/— 4 3 5 \ 

a) A = ( 3 0 20 

\ 5 -2 11/ 


/ 3 -1 "4\ 

b) A = 5 -2 -1 
\2 -1 2 / 


3. Define S for the matrices of Ex. 2 and let C be one or the other of the curves 
of all vectors (points) u = (jtj, x* t 2 * 3 ) whose coordinates satisfy the following equa¬ 
tions. Find the equation of the curves C s into which each S carries each C. 

а ) 3 j? - 2 x ] + 2 xl + 4*1X2 ~ 2xix 8 ~ 2*2X3 = 1 

б ) —4*i + llx§ = - 6 * 1 X 2 - 10*1*3 — 18*2*3 + 1 


* Observe that a change of bases (47) of L defines a linear transformation of L when 
we put uf * u\ 0) . Thus we may regard a change of basis as being induced by a nonsingu¬ 
lar linear transformation. 
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10. Orthogonal linear transformations. The final topic of our study of 
linear spaces will be a brief introduction to those linear transformations per¬ 
mitted in what is called Euclidean geometry. Let then V n be the n-dimen- 
sional linear space of vectors u = (ci, . . . , c„) over a field F. We define 
the norm of u (square of the length of u) to be the value f{u) = f(ci, . . . , c„) 
of the quadratic form 

( 52 ) /(*„ + 4 - 

We propose to study those linear transformations S on V n which are said 
to be length preserving and define this concept as the property that/(?/) = 
f(u s ) for every u of V n . Such transformations S will be called orthogonal. 

We may define S by u s = uA for an n-rowed square matrix A. Then, 
clearly, f(u) = uu', f(u s ) = u s (u s )' = uAA'u'. Write B = AA' = (6„) 
and see that f(u s ) — ZcibijC,, f(u) = 2 c\. Put c p = 1, cy = 0 for j 9 ^ p and 
have f(u) = f(u s ) only if bn = 1 for every i. Then take c v = c q = 1, all 
other Cj — 0, from which f(u) = 2, f(u s ) = 2 + 26 Pfl = f{u) only if 
b pq = 0 for a p q } B = I is the identity matrix. We call matrices A 
satisfying 

(53) AA' = I 

orthogonal matrices and have shown that S is orthogonal if and only if its 
matrix is an orthogonal matrix. 

We have seen that S determines A uniquely only in terms of a fixed 
basis of V n . Now in Euclidean geometry the only changes of basis allowed 
are those obtained by an orthogonal linear transformation, that is, those 
for which the matrix P of (51) is orthogonal. But then PP' = /, P' = P~ l , 
P'P = /, so that if S is a linear transformation with orthogonal matrix A 
and we replace A by PAR' 1 = B } then 

BB' = ( PAP'){PA'P ') = PAA'P' = I , 
and B is also orthogonal. 

EXERCISES 

1. What are the possible values of the determinant of an orthogonal matrix? 

2. Let I be the identity matrix, A be a skew matrix such that I + A is nonsingu¬ 
lar. Show that (7 + AR^I — A) is orthogonal. 

3. Show that every orthogonal two-rowed matrix A has the form 

(±5 To) . 
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where a = s(s 2 + t 2 )~ 1/2 f b = t(s 2 + s and t range over all quantities in F 

such that s 2 + t 2 0 . Hint: The result is trivial if A is a diagonal matrix. Now 
afi + a ?2 = 1 so that s = an, t = an are of the form above, while, conversely, 
a 2 + b 2 = ( s 2 + t 2 )(s 2 + t 2 )~ l — 1. The values 021 = ±6 ,022 = + a are derived from 
AA ' « J. 

4. The equations for rotation of axes in a real Euclidean plane through an angle 
h are x = x 0 cos h — y 0 sin h, y = x 0 sin h + y 0 cos h . Show that the corresponding 
matrix is orthogonal and that every real orthogonal two-rowed matrix is either the 
matrix of a rotation or its product by the matrix 



of a reflection x = Xq, y = — y 0 . 

5. Find a real orthogonal matrix P for each of the following symmetric matrices 
A such that PAP ' is a diagonal matrix. Hint: Compute PAP ' = D = (<2»,) and 
put dn = 0. 

»> C l) *> (1 -J) 0 (! !) « G 3 

11. Orthogonal spaces. We shall call two vectors u and v of V n orthog¬ 
onal (that is, in a geometric sense, perpendicular) if uv f = 0. Then, if A is 
any orthogonal matrix and we define a linear transformation S of V n by 
u s = a A, we have 

u s = uA , V s = vA , n s (v s ) f = uAA'v' = uv r = 0 

if and only if uv' = 0. Thus orthogonal transformations on V n preserve 
orthogonality of vectors of F n . 

If L = ViF + . . . + v m F is a linear subspace of V n , we define the set 
O(L) to be the set of all vectors w in V n such that vw r = 0 for every v of L. 
Then O(L) is a linear subspace of V n which we shall call the space orthogonal 
in V n to L. For, clearly, if u\ and w 2 are in 0(L) and a and b are in F, we 
have v(au\ + biv 2 )' = avw[ + bvw f 2 = 0, aw x 4 bw 2 is in 0(L). We now 
prove 

Theorem 9. Let L be a linear subspace of order m of V n . Then the order 
of O(L) is n — m. 

In fact we shall prove 

Theorem 10. Let L be the row space of the m by n matrix G of rank m so 
that G = (I m 0)Q for a nonsingular n-rowed mairix Q. Then O(L) is the 
row space of H = (0 I n -m)(Q') _1 of rank n - m. 

For proof we first note that the elements of GH f are the products i\w] 
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of the rows !>,• of G by the transposes of the rows w>, of H. But GH' — 
(I m 0)[(0 /»_*)], and, clearly, then GH' = 0. Let then O(L) = y x F + 

. . . + y q F of order g over F so that O(L) contains the row space of H. 
Evidently H has rank n — m, and we must have q > n — m. We let H a be 
the g by n matrix whose fcth row is ijk, and we have GH' 0 = 0, 0 = 
(/„ 0)QH' 0 . The rank of the n by g matrix 

d -««; - (d|) 

is q since Q is nonsingular. But 

(/„ o )z> = (£'), 

so that Di = 0, D has at most n — m nonzero rows, and D must have rank 
q < n — m. This proves that q = n — m , the row space of H is O(L). 

Note that the row space of H is the set of all solutions x — (xi, . . . , x n ) 
of the homogeneous linear system Gx f = 0. 

The sum of the orders of L and O(L) is n, and it is natural to ask whether 
or not they are complementary in V n , that is, whether V n = L + O(L). 
This is not true in general, since if F is the field of all complex numbers and 
L = vF, i 2 = — 1, v = (1, i), then w' = 1 + i 2 = 0. Hence O(L) is con¬ 
tained in L and O(L) = L. However, we may prove 

Theorem 11. Let F be a field whose quantities are real numbers . Then 
V n = L + O(L). 

For by Theorem 9 it suffices to prove that the only vector w in L and 
O(L) is the zero vector. Hence let w be in both L and O(L) so that w = dG } 
where d = (di, . . . , d m ) has elements in F and G is an m by n matrix of 
rank m whose row space is L. Then Gw’ — 0 while also Gw' = GG'd' . By 
Theorem 3.17 the matrix GG' has rank m and is nonsingular, GG'd ' = 0 
only if d' = 0, d = 0, w = 0 as desired. 


EXERCISE 

Let L be the space over the field of all real numbers spanned by the following 
vectors U{, Find a basis of the space O(L) in the corresponding V n . 


а) u\ = (1, 2, — 1, 0) f 

б) ix! - (1, 0, 1, 1) , 

c) ui - (1, -1, 2, 1) , 

d) txi = (1, 2, -1) , 


u, = (0, 1, 2, 1) 

U2 = ( 0 , 1 , 0 , 1 ) , 

IX* = (2, -1, 2, 3) , 

f-1. 1,0), 


ix 3 = (-1,2, 1,0) 
ix 3 = (1, —2, 4, 0) 
w 3 = (3, 3, 3) 



CHAPTER V 

POLYNOMIALS WITH MATRIC COEFFICIENTS 


1. Matrices with polynomial elements. Let F be a field and designate by 
0) F[x 1 

(read: F bracked x) the set of all polynomials in x with coefficients in F. 
We shall consider m by n matrices with elements in F[x ] and define elemen¬ 
tary transformations of throe types on such matrices as in Section 2.4. As 
we stated in that section, we assume that in the elementary transformations 
of type 2 the quantities e art' permitted to be any quantities of F[x]. Rut 
those of type 3 are restricted so that the quantity a in F[x] shall have an 
inverse in F[x]. Then a ^ 0 must be a constant polynomial, that is, a may 
be any nonzero quantity of F. 

We now let A and B be m by n matrices with elements in F[x] and call 
A and B equivalent in F[x\ if there exists a sequence of elementary transfor¬ 
mations carrying A into B. The field Fix) of all rational functions of x with 
coefficients in F contains F[x], and it is thus clear that if A and B are equiva¬ 
lent in F[x | they are also equivalent in Fix). Hence we see that A and B 
arc 1 equivalent in F[x] only if they have the same rank. We may then 
prove* 

Lkmma 1. Fvery nonzero matrix A of rank r with elements in F[x] is equiva¬ 
lent in F[xl to 



where (L = diag jfi, . . . , f r | for manic polynomials U = L(x) such that fj di¬ 
vides f^ i. 

For the elements of all matrices equivalent in F[x] to A are polynomials 
in x, and in the set of all such polynomials there is a nonzero polynomial 
fi = f\{x) of lowest degree. Using elementary transformations of types 3 
and 1, we may assume that/i is monic and is the element in the first row 
and column of a matrix C = (a,) equivalent in F[x] to A. By the Division 
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Algorithm for polynomials we may write c,i = qj\ + r x for q x and r» in 
F[j] and r, of degree less than the degree of f\. But if w'e add —< 7 , times the 
first row of C to its ith row, we pass to a matrix equivalent in F[x] to A 
with r. as the element in its ith row and first column. Our definition of/i 
thus implies that r 4 is zero. Moreover, we have now shown A equivalent 
in F[x] to a matrix D = (tf„) with d lX = 0 for it* 1, dn = / 1 . Similarly, 
we see that every d u is divisible by fi and hence that A is equivalent in 
F[j*] to a matrix 



where A 2 has m — 1 rows and n — 1 columns. Then either A 2 = 0 , or we 
may apply the same process to A 2 . After a finite number of such steps we 
ultimately show that A is equivalent in F[x] to a matrix (2) of our lemma 
such that every/* = f x (x) is monic and is a polynomial of least degree in 
the set of all elements of all matrices equivalent in /'V] to 

(4) -1. = 0 ), — fliag !/i. fA 


Write / 1+ 1 = f x Si + t t , where 6 *, and /, are in F[x] and the degree of t t is less 
than the degree of /,. Then we add the first row of A x to its second row so 
that the submatrix in the first two rows and columns of the result is 



We then add —a* times the first column to the second column to obtain 
a matrix equivalent in /'[x] to Ai and with corresponding subrnatrix 



Our definition of f t thus implies that t x = 0,/» divides f xvx as described. 

We observe now that the elements of every /-rowed minor of a matrix A 
with elements in F[x ] are polynomials in x, so that these minors are also 
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polynomials in x. By Section 2.7 if B is obtained from A by an elementary 
transformation, then every /-rowed minor of B is either a /-rowed minor of 
A, the product of a /-rowed minor of A by a quantity a of F, or the sum 
M\ + /M 2 where Mi and M 2 are /-rowed minors of A and/ is a polynomial 
of F[x]. If d in F[x] then divides every /-rowed minor of A, it also divides 
every /-rowed minor of all m by n matrices B equivalent in F[x] to A. But 
A is also equivalent in F[x] to B and therefore the /-rowed minors of equiva¬ 
lent matrices have the same common divisors. We may state this result as 

Lemma 2. Let A be an m by n matrix with elements in F[x] and d t be the 
greatest common divisor of all t -rowed minors of A. Then d t is also the greatest 
common divisor of the t-rowed minors of every matrix B equivalent in F[x] to A. 

Every (k + 1)-rowed minor M k + X of A may be expanded according to a 
row and is then a linear combination, with coefficients in F[x], of Arrowed 
minors of A. Hence d k divides every M k+ i so that d k divides d k + 1 . We ob¬ 
serve also that in (2) the only nonzero Arrowed minors are the Arrowed 
minors |diag [f iv ... ,/<J | for U<i 2 < . . . <i k , and since clearly every f { 
divides /,+»we see that the g.c.d. of all Arrowed minors of (2) is fi . . . f k . 
We thus have d k = fi . . . /*, whence 


(7) 


d k -i 


= /* 


It is customary to relabel the polynomials /< and thus to write f r = g i, 
f r _ x = g 2f fi = g r . We call gr, the jth invariant factor of A and see that if we 
define d 0 = 1 we have the formula 


( 8 ) 


9i = 


dr - 


7+1 


dr—j 


(j = 1 , . . . , r) 


Moreover, g , is now divisible by 0,+1 for j = 1, . . . , r — 1. We apply ele¬ 
mentary transformations of type 1 to ( 2 ) and see that if A has invariant 
factors pi, . . . , 0 r , then A is equivalent in F[x] to 

(9) (o o) ’ G = diag {9x ’ • * • ’ 9r} ' 

If, then, B is equivalent in F[x] to A, it has the same greatest common di¬ 
visors dt and hence the same gj of ( 8 ), while, if the converse holds, B is 
equivalent in F[i] to (9) and to A. We have proved 
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Theorem 1. Two m by n matrices with elements in F[x] are equivalent in 
F[x] if and only if they have the same invariant factors. Every m by n matrix 
of rank r with invariant factors gi, . . . , g r is equivalent in F[x] to ( 9 ). 

As in our theory of the equivalence of matrices on a field we may obtain 
a theory of equivalence in F[x ] using matrix products instead of elementary 
transformations. We thus define a matrix P to be an elementary matrix if 
both P and P " 1 have elements in F[x\. Then, if a = \P\ and b = I/*" 1 !* 
the quantities a and b are in F[x], ab = 1 1 | = 1, so that a and b are nonzero 
quantities of F. But if P has elements in F[x] so does adj P. Hence so will 
P- 1 = |P|~ l adj P if |P| is in F. Thus we have proved that a square ma¬ 
trix P with polynomial elements is elementary if and only if its determinant 
is a constant (that is, in F ) and not zero. 

We now observe that, in particular, the determinant of any elementary 
transformation matrix is in F. Hence, if A and B are square matrices with 
elements in F[x] and are equivalent in F[x], their determinants differ only 
by a factor in F. Moreover, if |A | and \B\ are monic polynomials, then 
the equivalence of A and B in F[x] implies that | A | = \B\ and, in fact, 
that when A is a nonsingular matrix with <71, . . . , g n as invariant factors its 
determinant is the product g\ ... g n . 

It is clear now that a square matrix P with elements in F[x] is elementary 
if and only if P is equivalent in F[x] to the identity matrix, and thus the 
invariant factors of P are all unity. We may now redefine equivalence. We 
call two m by n matrices A and B with elements in F[x] equivalent in F[x] if 
there exist elementary matrices P and Q such that PAQ = B. Then we again 
have the result that A and B are equivalent in F[x] if and only if they have 
the same invariant factors. For under our first definition P and Q are equiva¬ 
lent in F[x] to identity matrices and hence may be expressed as products of 
elementary transformation matrices. But, if P 0 and Q 0 are elementary trans¬ 
formation matrices, the products PqA, AQ 0 are the matrices resulting from 
the application of the corresponding elementary transformations to A. 
Hence, PAQ must have the same invariant factors as A. The converse is 
proved similarly, and we have the result desired. 

In closing let us note a rather simple polynomial property of invariant 
factors. The invariant factors of a matrix A with elements in F[x\ and rank 
r are certain monic polynomials 0,(2) such that <7»+i(x) divides g^x) for 
i = 1, . . . , r — 1. If g k {x) = 1 , then gr,(x) = 1 for larger j = A; + 1 , 
. . . , r. Let us then call those Qi(x) ^ 1 the nontrivial invariant factors 
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of A, the remaining g<(x) = 1 the trivial invariant factors of A. Thus there 
exists an integer t such that g,(x) has positive degree for i = 1, . . . , t, 
9t+i(x) = ... = g T {x) = 1. 


EXERCISES 

1. Express the following matrices as polynomials in x whose coefficients are 
matrices with constant elements. 


jl + 2x x 3 + 4x 2 + x + 2 
a) I 0 x 2 + x 

\1 — 2x x 3 + 3x 2 — 3a; — 1 


x 3 + 4x + 2 
x 2 

x 3 — x 2 + 4x — 2 


b) 


' x 2 
x 2 - 1 

<x 3 -h 2a; 2 - 2 


x 

1 + x 

x 2 + 2a; + 2 



/x 3 + x 2 — 2a; 
c) ( X 3 

\2x 3 + x 1 - 2x + 1 


X 2 - X 3 x 3 + x\ 

X -X 3 X 3 ) 

x 2 — x — 1 X + 1/ 


( X 2 + X X 

X 2 + a; 1+1 

a; 2 - 1 2a; 


z 2 + l\ 
x 2 + x 1 
2x 2 + 2/ 


2. Let -A be an m by n matrix whose elements are in F[x]. Describe a process by 
means of which we may use elementary row transformations involving only the 
ith and Arth rows of A to replace A by a matrix with the g.c.d. of a» 7 * and auj in its 
ith row and jth column. Hint: Use the g.c.d. process of Section 1.6 with / = a t7 , 
g == cikj. 

3. Use elementary transformations to carry the following matrices into the 
form (9). 




x(x - 1) 
0 


0 

x(x + 1) 


) 


* + x - 2 
0 


0 

x* + 2x 


-a) 


4. Describe a process for reducing a matrix A with elements in F[x] to an equiva¬ 
lent diagonal matrix. How, then, may we use the process of Ex. 3 to carry this 
preliminary diagonal matrix into the form (9)? 

5. Reduce the matrices of Ex. 1 to the form (9) by the use of elementary trans¬ 
formations. 

6. Determine the invariant factors of the matrices of Ex. 1 by the use of (8). 
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7. Use elementary transformations to reduce the following matrices to the 
form (9). 



2 -h x — x 3 — x 4 
2 + x + x 2 — 2x 4 


x - 2 
1 + 2x - x 2 - x 3 
-1 + x - 2a; 3 


0 


2x — 2 
— 3x 
1 — x 

— 2x 2 + x — 2 


x 2 - 2 
4x + 6 
x+l 
4x 2 + 5x - 


^2x 2 + 4x + 2 
x 2 + x — 3 
3x + 6 
x + 2 


3x 2 + 3x 

2x + 3 

x + 3' 

2x 2 

2x x 

3x 2 -f x 

2x + 1 

x + 1 

x — 2x 2 

1 - 2x 

1 — X 

2x + 1 

x — 1 

x — 1 

x + 2 

3x -}- 4 

X 

x 2 4- x -b 2 

X 2 + 3x - 1 

— x 

2x 2 + 3x + 1 

5x 2 + 2x 

x 2 - 1 

-2x 2 + x 

3x + 3 

-x 2 

x + 1 



6x 2 + 5x + 2 
-2x 2 + 2x + 1 
—x 2 


2. Elementary divisors. Let K be the field of all complex numbers so that 


9 i(x) = (x - ci) e » . . . (x - c,) e *, 


where Ci, . . . , c, are the distinct complex roots of the first invariant factor 
of a matrix A with elements in K[x]. Since g i+ i(x) divides < 7 t(x), it is clear 
that every g»(x) divides fifi(x), and thus 

(10) £<(x) = (x - . . . (x - c.) e " (i = 1, . . . , r) . 

Here r is the number of invariant factors of A, €\,• = 6„ e »,• ^ 0 for i = 2, 
. . . , r. We shall call the rs polynomials 

/.7 = (x - Cy)**7 

the elementary divisors of A. Those for which > 0 will be called the 
nontrivial elementary divisors of A. 
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The invariant factors of A clearly determine its elementary divisors 
uniquely as certain rs powers / i7 of linear functions x — c, where r is the 
rank of A, and s is the number of distinct roots of the invariant factors 
of A. Conversely, the elementary divisors of A uniquely determine its in¬ 
variant factors. In fact, let us consider a set of q polynomials each a power 
with positive integral exponent of a monic linear factor with a complex 
root. The distinct roots in our set may then be labeled Ci, . . . , c, and our 
polynomials have the form hkj = (x — c ,)"0 for n*, >0. For each c, let tj be 
the number of hkj in our set, t be the maximum tj. Then, clearly, q = ti + 
...+*« ^ ts, and our set of polynomials may be extended to a set of 
exactly ts polynomials by adjoining t — tj polynomials (x — c,-) n< * with 
n»; = 0. Let us then order the t exponents n^ to be integers e t y satisfying 
e X j e 2 j ^ ^ e ti ^ 0. Define < 7 t(x) as in (10) for i = 1, . . . , t and 

obtain a set of polynomials Qi{x) such that ^ l+ i(x) divides < 7 »(x) for i — 1, 
. . . , t — 1 , the < 7 i(x) are the nontrivial invariant factors of a matrix A 
whose nontrivial elementary divisors are the given hkj. If A has rank r, we 
have r ^ t } and we adjoin (r — t)s new trivial elementary divisors /„ to 
obtain the complete set of r invariant factors g\(x) of A. 

It is now evident that two m by n matrices with elements in K[x] are 
equivalent in K[x] if and only if they have the same elementary divisors. 

The matrix (9) has prescribed invariant factors and hence prescribed 
elementary divisors. However, it is desirable to obtain a matrix of a form 
exhibiting the elementary divisors explicitly. We shall do this. Let us prove 
first 

Theorem 2. Let fi, . . . , f g be monic polynomials of F[x] which are rela¬ 
tively prime in pairs. Then the only nontrivial invariant factor of the matrix 
A = diag {fi, . . . , f s } is its determinant g = fi . . . f». 

The result is trivial for s = 1. If 5 = 2, the g.c.d. of the elements f\ and 
/ 2 of A 2 = diag {/i, / 2 ) is unity, d i = l,/i/ 2 is the only nontrivial invariant 
factor of A 2 , and A 2 is equivalent in F[x] to diag {fif 2 , 1}. Assume, then, 
that A,_i = diag \f h . . . ,ij is equivalent in F[x] to i = diag 
1, . . . , 1}, where g = fi . . . /— i. Then A = diag {/i, ...,/,) is equiv¬ 
alent in F[x] to diag /., 1, . . . , 1}. But g a -i is prime to /„ diag 

[g t -i, /,} is equivalent in F[x] to diag {< 7 , lj. Hence A is equivalent to 
diag {< 7 , 1, . . . , 1}. Then g is the only nontrivial invariant factor of A, 
and our theorem is proved. 

We see now that if gi(x) is defined by (10) for distinct complex numbers 
Cj f the corresponding elementary divisors /,i, . . . , /»* are relatively prime 
in pairs. By Theorem 2 the matrix Ai = diag {/n, . . . ,/»} is equivalent 
in F[x] to diag {gi , But then A = diag {A 1 , . . . , A t } is equiv- 
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alent in E[x] to diag {ffi, . . . , Qt, 1, . . . , 1}, and hence the nontrivial in¬ 
variant factors of A are gi,, Qt . We examine the form of A to see that 
we have proved 

Theorem 3. Let Ci, . . . , c r be complex numbers and f, = (x — Ci)“i for in¬ 
tegers ni Si 0. Then the matrix 

A = diag {f l( . . . , f,} 

has the fi as its elementary divisors. 

EXERCISES 

1. The following polynomials are the nontrivial invariant factors of a matrix. 
What are its nontrivial elementary divisors? 

a) x* + 2x 4 + x 3 , x 3 + x 1 , x 3 + x 

b) X* + x 3 + 2x 4 + 2x 3 + x 3 + x, x 3 + x, x 

c) x(x - l) 2 (x 2 - 2x + 1), x 3 - 2x 2 + x, x - 1 

d) x 4 - 5x 3 + 9x 2 -7x + 2, x 3 - 4x 2 + 5x - 2, i* - 3x + 2 

e) (x> - l) 3 , (x* - 1)*, (x* - 1), x - 1 

/) (x*+l) 4 , (x’+l) 5 , x*+l 

2. The following polynomials are the nontrivial elementary divisors of a matrix 
whose rank is six. What are its invariant factors? 

a) (x - 1)’, (x - 1)*, (x - 1), (x + 1)*, (x + 1) 

b) (x - 2)*, (x - 2) 3 , (x - 2), x, (x + 1)*, x 2 

c) (x — 3), (x — 3) 3 , (x — 3) 5 , x 2 , x 4 , x 4 

d) x, (x - 1), (x - 2), (x - 3), (x - 4), (x - 5) 

e) (x + 1)», (x + l) 2 , (x - 1), x, x 2 , x 2 , x 3 

/) x, (x - 1), x 2 , (x - l) 2 , x 3 , (x + l) 4 , x 2 

3. Find elementary transformations which carry the following matrices into the 

form (9). 

/(x - l) 3 0 0 \ /x 3 0 0 \ 

a) I 0 x — 2 O) i) (0 x + l 

\ 0 0 x - 1/ \0 0 x + 2/ 

/x 3 0 ON 

c) (0 (x - l) 2 0 I 

\0 0 x + l/ 
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3. Matric polynomials. Let the elements a it of an m by n matrix A be 
in F[x] and let 5 be a positive integer not less than the degree of any of the 
an. Then every a i3 has s as a virtual degree and we may write 

(11) an = a[fx * + a[?x-' + . . . + a<‘> 

for a\f in F. Define A k — (a[f) and obtain an expression of A in the form 

(12) /Or) = A 0 x* + . . . + A, 

for m by n matrices A*. Thus f(x) is a polynomial in x of virtual degree 
s with m by n matrix coefficients A* and virtual leading coefficient A 0 . 
Moreover, we say that/(x) has degree s and leading coefficient A 0 if A 0 9 * 0. 

In order to be able to multiply as well as to add our matrices we shall 
henceforth restrict our attention to n-rowed square matrices with elements 
in F[x\ and thus to the set of all polynomials in x with coefficients n-rowed 
square matrices having elements in F. Let us call these polynomials n- 
rowed matric polynomials. If all the A k in (11) are zero matrices, the poly¬ 
nomial/^) is the zero polynomial, and we again designate it by 0. Evident¬ 
ly we have 

Lemma 3. The degree of f(x) + g(x) is not greater than the degree of f(x) 
or g(x). 

Lemma 4. Let f(x) have degree n and leading coefficient A 0 , g(x) have de¬ 
gree m and leading coefficient B 0 such that AoB 0 ^ 0. Then the degree of 
f(x)g(x) is m + n and the leading coefficient of f(x)g(x) is AoB 0 . 

As in Chapter I we use Lemma 4 in the derivation of the Division Algo¬ 
rithm for matric polynomials which we state as 

Theorem 4. Let f(x) and g(x) be n-rowed matric polynomials of respective 
degrees s and t such that the leading coefficient of g(x) is nonsingular. Then 
there exist unique polynomials q(x), Q(x), r(x), R(x) such that r(x) and R(x) 
have virtual degree t — 1 and 

(13) f(x) = q(x)g(x) + r(x) = g(x)Q(x) + R(x) . 

Moreover , if s < t, then q(x) = Q(x) = 0, while if s ^ t, then q(x) and 
Q(x) have degree s — t. 

While the proof is very much the same as that in Theorem 1 . 1 , let us 
give it in some detail. We assume first that s ^ t and let A 0 5 ^ 0 and B 0 be 
the respective leading coefficients of f(x) and g(x). Then Bq 1 exists, and if 
q x (x) = A 0 Bq 1 x 9 ~* the polynomial f\(x) = f(x) — qi(x)g(x) has virtual de- 
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gree s — 1. This implies a finite process by means of which we begin with 
a polynomial /»(x) of virtual degree s — i^t and leading coefficient A 
and then form /»+i(x) = /*(x) — q%+\{x)g{x) of virtual degree s — i — 1 
for g<+i(x) = A^B^x 4- *” 1 . The process terminates when we obtain an 
fi(x) of virtual degree t — 1. Thus we have the first equation of (13) with 
q(x) = 0, r(x) = /(x) if s < t, and otherwise with q{x) of degree 8 — J and 
leading coefficient AqBq 1 and r(x) the/,(x) above. Ifalso/(x) = qo{x)g{x) + 
r 0 (x) for q 0 {x) — q(x) 0, the degree of this polynomial is h ^ 0, 4 and its 
leading coefficient is C 0 0. Also C 0 B 0 s* 0 since B 0 is nonsingular. But 
then by Lemma 4 the degree of [g 0 (x) — q{x)]g{x) = r(x) — r 0 (x) is t + h, 
whereas its virtual degree is t — 1, a contradiction. This proves the unique¬ 
ness of g(x) and r(x). The existence and uniqueness of Q(x) and R{x) is 
proved in exactly the same way except that we begin by forming f{x) — 
g(x)BQ l A 0 x 9 ~ t . 

Let us regard the first equation of (13) as the right division of f{x) by 
g{x) f the second as the left division of f{x) by g{x). Then we shall speak 
correspondingly of q{x) and r(x) as right quotient and remainder, of Q{x) 
and R{x) as left quotient and remainder. If r(x) = 0, we have /(x) = 
q{x)g{x) y and g(x) is a right divisor of /(x). Similarly, we call g{x) a left 
divisor of g{x) if /(x) = g(x)Q(x), so that R{x) = 0 in (13). 

It is natural now to try to prove a Remainder Theorem for matric poly¬ 
nomials. However, the theorem in usual form is ambiguous since, for ex¬ 
ample, if C is an n-rowed square matrix and/(x) = A 0 x 2 = x 2 A 0 = xAox, 
the polynomial f(C) might mean any one of A 0 C 2 , C 2 A 0 , CA 0 C, and these 
matrices might all be different. Thus we must first define what we shall 
mean by /(C). We shall do this and obtain a Remainder Theorem which we 
state as 

Theorem 5. Let f(x) be an n-rowed square matric polynomial (12), C be 
an n-rowed square matrix with elements in F. Define f R (C) {read: i right of C), 
and fi,(C) {read: f left of C) by 

(14) f R (C) = AoC- + A X C- 1 + . . .+A i _ 1 C+A 1 , 
and 

(15) fi/C) = C"Ao +C*” 1 Ai + . . . +CAa_i+ A,. 

Then the right and left remainders on division of f(x) by xl — C are f R (C) 
and fi,(C), respectively. 
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To make our proof as in the case of polynomials with coefficients in a 
field we use Theorem 4 with g(x) = zl — C and have /(x) = q(x)g(x) + 
r(x) where q{x) has degree s — 1 and r(x) = B has elements in F. We then 
wish to draw our conclusion from/«(C) — ^/e(C)(C — C ) + B. This state¬ 
ment is correct, but let us examine it more closely. We write q(x) = 
Cox*” 1 + . . . + C,_ i and have 

h(x) = q(x)(xl - C) 

- (Cox* + Cix-* + . . . + C^x) - (CoCx- 1 + . . . + C.- X C) . 
Then, if D is any n-rowed square matrix with elements in F } we have 
h R (D) = CoZ> + (Ci - C 0 C)Z>~ 1 + . . . + (C*_i - C._*C)D - C.^C, 
while 

q R (D)(DI - C) - C 0 D‘ + ( C X D> - CoD-'C) + . . . 

+ (C f _xD - C._*DC) - C,_xC , 


and these matrices are equal in general* if and only if D and C are commu¬ 
tative. They are equal if D = C, and thus /(x) = h(x) + B implies that 
f R (C) = Aft(C) + 5 = q R (C)(C — C) + B = B. The second part of our 
theorem is proved similarly. 

As a consequence of the result just proved we have the Factor Theorem 
for matric polynomials which we state as 

Theorem 6. The matric 'polynomial f(x) has xl — C as a right divisor if 
and only if f R (C) = 0; it has xl — C as a left divisor if and only if fj,(C) = 0. 

For Theorem 5 implies that, if f R (C) = 0, then in (12) the polynomial 
r(x) — 0, /(x) = q(x)(xl — C). Conversely, if /(x) = q(x)(xl — C), we 
have seen that/a(C) = q R (C)(CI — C) = 0. The results on the left follow 
similarly. 

Our principal use of the result above is precisely what is usually called 
the trivial part of Theorem 5, that is, if /(x) has x — C as a factor, then 
C is a root of /(x). However, it is nontrivial that f R (C) = 0 and follows 
only from the study above where we showed that if D is any square matrix 
such that DC = CD, then h(x) = q{x)(xl — C) implies that h R (D) = 
q R (D)(D - C). 

* E.g., if q(x) - iso that h(x) - x* — Cx , then h R (D) — D* — CD, q R (D)(D — C) « 
D 1 — DC & D* — CD unless DC - CD. 
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EXERCISES 

1. Express the following matrices as polynomials f(x) and g(x) with matric coeffi¬ 
cients and compute q(x), Q(z), r(x), R(x) of (13). 


.... / + 5x + 1 3x> + x — 1 \ , , _ /2x’ 

a) Six) - \2x' + x* +2 4x s + 2x + 2/ ' “ Us 2 


b) fix) 


e 


+ 3x s — 1 

+ 1 


t + l) > ?(x) " (- 


3x 2 + x 
x 2 + 2x 


9(z) = 



- 1 x*\ 

2x 2 ) 

-2x 2 + x + 1> 
x 2 + x ) 


/2x 4 — x 2 + 2 —x 8 + x — 1 1 — x 2 

c) /(s) = I x s — x + 1 —x 4 + x 2 — 2 1 + x 

Vx 2 -1 —1 — x x 4 + x 2 -l/ 

/x 2 + 2x + 1 

1 — X 

2x 


d) fix) = 


2 

4 

3x 2 

2x 2 — x 


x 2 - 1 
2x 2 + x 
x 2 + 1 


2x 2 - 
x 2 - 
1 


x 2 + 1 

X 2 - 1 
X 2 + 2X/ 


3x 2 4- n 
x 2 + x 


2. Use/(x) of Ex. 1(c) and find//*(£) and /l(C) by the use of the division process 
as well as by substitution if 


a) C 


0 1 0 

0 0 1 

2 0 0 . 


b) C 


/0 2 0 \ 
" 1 0 0 
Vo o i/ 


(l 2 
c) C = 2-1 

Vo 0 


4. The characteristic matrix and function. If fix) is a matric polynomial 
(12) with n-rowed scalar matrix coefficients A k , then we shall call/(x) a 
scalar 'polynomial. Thus A * = a*J for the a* in F , and 

(16) fix) = (a 0 x* + . . . + a,)/ , 

where I is the n-rowed identity matrix. We call fix) monic if a 0 = 1. If 
now gix) is also a scalar polynomial, the quotients qix) and Q(x) in (13) 
are the same scalar polynomials and also r(x) = R(x) is scalar. For obvi¬ 
ously this case of Theorem 4 is now the result of multiplying all the poly¬ 
nomials of Theorem 1.1 by the n-rowed identity matrix. 
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If A is any n-rowed square matrix, the polynomials fn(A) and / L (A) are 
equal for every scalar polynomial/(x), and we shall designate their com¬ 
mon value by 

(17) f(A) = aoA 9 + . . . + + aj . 

We now say that either the polynomial/(x) or the equation/(x) = 0 has A 
as a root if /(A) = 0. By Theorem 6 the matrix A is a root of /(x) if and 
only if /(x) has x7 — A as either a right- or a left-hand factor. 

We shall call the matrix xl — A the characteristic matrix of A, the de¬ 
terminant | xl — A | the characteristic determinant of A, the scalar poly¬ 
nomial 

(18) fix) - |x7- A| •/ 

the characteristic function of A, and the corresponding equation/(x) = 0 
the characteristic equation of A. We now apply (3.27) to xl — A to obtain 

(19) 0 xl - A)[adj {xl - A)] = [adj {xl - A)](x7 - A) 

= |x7- A | •/. 

Then the elements of adj {xl — A) are the cofactors of the elements of 
xl — A, and adj {xl — A) is a matric polynomial (in general, nonscalar). 
By the argument above we have 

Theorem 7. Every square matrix is a root of its characteristic equation. 
The g.c.d. of the elements of adj {xl — A) is clearly the polynomial d n -i(z) 
defined for xl — A and d n {x) = |x7 — A|. But then adj (x7 — A) =d n ^i{x)B{x) f 
where B{x) is an n-rowed square matrix with elements in F[x]. Hence B{x) 
is a matric polynomial. The invariant factors of x7 — A were defined in (8), 
and by (19) |x7 — A|7 = gi{x)d n -\{x)l = {xl — A)B{x)d n - i(x). By the 
uniqueness of quotient in Theorem 4 we have 

(20) g{x) = <7i(x)7 = {xl - A)B{x) . 

Hence, clearly, g{A) = 0. Observe also that the g.c.d. of the elements of 
B{x) is unity so that if B{x) = Q{x)q{x) for a monic scalar polynomial 
q{x), then q{x) = 1. 

We now define the minimum function of a square matrix A to be the 
monic scalar polynomial of least degree with A as a root. The remark just 
made above then implies 

Theorem 8. Let gi(x) he the first invariant factor of the characteristic ma¬ 
trix of A. Theii g(x) = gi(x)I is the minimum function of A. 
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For if h(x) is the minimum function of A we may write g(x) = h(x)q(x) + 
r(x) for scalar polynomials h{x) and r(x) such that the degree of r(x ) is less 
than that of h (x). But/ 1 (A) = 0, and gr(A) = 0 so that r(A) = 0, r(x) = 0. 
Hence g{x) = A(x)g(x), and since g{x) and h(x) are monic so is q(x). By 
Theorem 6 we have h(x) = (xl — A)Q(x ), and by ( 20 ) we have g(x) — 
(x7 — A)Q(x)q(x) = (x7 — A)B(x). The uniqueness in Theorem 4 then 
states that B(x) = Q(x)q(x), q(x) = 1 , from which g(x) = h(x) as desired. 

We see now that \xl — A | is a monic polynomial of degree n and is not 
zero, r = m = n in (9), and 

(21) P{x) • (x7 - A) • Q(x) = diag {gi, . . . , 0 n } 

for elementary matrices P(x) and Q(x). But then, as we have already ob¬ 
served in an earlier discussion, 

( 22 ) c • | x/ — A | • d = gi. . . g n 

for c = | P(x) | and d = |Q(x) | in F. Hence cd — 1 , and we have proved 
Theorem 9. The characteristic function of a square matrix A is the product 
by I of the product of the invariant factors of the characteristic matrix of A. 

This result implies that |x7 — A | is the product of £i(x) by divisors 
</t(x) of g i(x). It follows that every root of |x7 — A | in any field K con¬ 
taining F is a root of gi(x). But in fact we have already seen that if F is 
the field of all complex numbers, the elementary divisors of x7 — A are 
polynomials (x — c,)* 0 ' whose product is |x7 — A|. Then the c, are the 
distinct roots of < 7 i(x) as well as of |x7 — A |. They are called the charac¬ 
teristic roots of A. 

In closing this section we note that if we write /(x) = |x7 — A | • 7 = 
(x n + aix n_1 + . . . + a*) 7 , then/( 0 ) = | —A | • 7 = (— l) n | A | • 7, so that 
|A | = ( —l) n a n . It follows that, if |A | 0 , then 

(23) A " 1 = ~a“ 1 (A n " 1 + a,A "' 2 + . . . + a n _ a • 7) , 

and hence it is obvious that AA~ l = A _ 1 A. Moreover, if |A| = 0 , then 
A -1 does not exist. If, then, A ^0 and g\{x) == x m + 6 ix m_1 + . . . + b m , 
the polynomial ^(x) = gi(x)I can be the minimum function of A only if 
b m = 0. Since A 5 * 0 , we have m > 1, and G = A m_1 + &iA m ~ 2 + . . . + 
b m -il is a nonzero matrix with the property 


(24) 


AG = GA = 0 . 
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EXERCISES 

1. When is the minimum function of a matrix linear? 

2. What, then, are the minimum functions of the following matrices? 

■»C4) «Ci) «G!) 

3. Let A = diag [A\ t A 2 \ where A\ and A 2 are square matrices of m and n rows 
respectively. Show that if /(x) is any polynomial in F[x] and we define /*(x) = 
f{x)l t for any t , then fm+n(A) = diag [f m (Ai), f n (A 2 )\. Hint: Prove first by induc¬ 
tion that A k = diag {Af, A%), 

4. Let A have the form of Ex. 3 and let g(x)Im+ n , 0i(x)L», and g 2 (x)I n be the 
respective minimum functions of A , Ai, A 2 . Prove that g(x) is the least common 
multiple of gi(x) and g 2 (x). 

5. Apply Ex. 4 in the case where A\ is nonsingular and A 2 = 0. 

6. Compute the characteristic functions of the following matrices. 

/I 2 3\ (2 3 1\ 

a) 2-14 b) 1 2 1 

\3 11/ \0 0 -4/ 


7. It may be shown that the characteristic function f(x) = x n — aix n_1 + 

. . . + (~l) i aiX n_< + . . . + ( — l) n a« of an n-rowed square matrix A has the 
property that a» is the sum of all t-rowed principal minors of A. Verify this for the 
matrices of Ex. 6. 

6. Similarity of square matrices. We have defined two n-rowed square 
matrices A and B with elements in a field F to be similar in F if there exists 
a nonsingular matrix P with elements in F such that 

PAP- 1 = B . 

The principal result on the similarity of square matrices is then given by 
Theorem 10. Two matrices are similar in F if and only if their charac¬ 
teristic matrices have the same invariant fadbrs. 

For if PAP- 1 = B, then P(xl - A)P~ l = xPIP~ l - B = xl - B. 
But P has elements in F, | P | ^ 0 is in F, P is elementary. Hence xl — A 



I 4 3 2 —1\ 

3 2 1 0\ 

l 2 1-1-2 

\-l 0 -2 -5/ 
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and x7 — B are equivalent in F[x ] and have the same invariant factors. 
Conversely, let xl — A and xl — B have the same invariant factors, so that 

(25) P(x)[xl - A]Q(x) = xl - B 
for elementary matrices P(x) and Q(x). We define 

(26) P = P l (B) , Q = Q r (B) 
as in Theorem 5 and have 

(27) P(x) = (x7 - B)P 0 (x) + P , Q(x) = Qo(x) (x7 — B) + Q. 

Then 

P(x)[(xJ - A)]Q(x) - (x7 - B)Po(x)(x7 - A)Q(x) + 

P • (x7 - A) • Q 0 (x)(x7 - B) + P • (x7 - A) • Q = x7 - B. 

We now use (25) and the fact that P(x) and Q(x) are elementary to write 
[P(x)] _1 = C(x), Q(x) -1 = 7)(x) for matrix polynomials C(x) and D(x) 
such that 

(28) (xl - A)Q(x) = C (x)(x7 - B), P(x)(x7 — A) = (x7 - B)D(x ). 

But then from (27) and (28) 

P • (x7 — A) = (x7 - B)[D(x) - P 0 (x)(x7 - A)] 

and thus 

(29) (x7 — B) — P • (x7 — A) • Q = (x7 - B)B(x7 - B) , 

where B = B(x) = P 0 (x)C(x) + 7>(x)Q 0 (x) — P 0 (x)(x7 — A)Qo(z). By Lem¬ 
ma 4 the degree in x of the right member of (29) is at least two unless 
R(x) — 0. But the degree in x of the left member of (29) is at most one, 
R(x) = 0, 

(30) P • (x7 - A) • Q = x7 - B . 

It follows that PAQ = B, P7Q = 7, Q = P _1 , PAP~ l = B as desired. 

Observe that the degree of | x7 — A | is n and hence that if n* is the de¬ 
gree of the ith nontrivial invariant factor gi(x) the property |x7 — A | = 
g i(x) . . . g t (x) implies that 

(31) n\ + n 2 + . * . + n t = n , n\ ^ ^ ... g; n* = 0 . 


Obviously, this is an important restriction on the possible degrees of the in¬ 
variant factors of the characteristic matrix of an n-rowed square matrix. 
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EXERCISES 

1. What are all possible types of invariant factors of the characteristic matrices 
of square matrices having 1, 2, 3, or 4 rows? 

2. Give the possible elementary divisors of such matrices. 

3. Use the proof of Theorem 10 to show that if A Xf A 2 , B h and B 2 are n-rowed 
square matrices such that A x and B x are nonsingular, then A x x + A 2 and B x x -f B 2 
are equivalent in F[x] if and only if there exist nonsingular matrices P and Q with 
elements in F such that PA X Q = B x and PA 2 Q = B 2 . Hint: Take A = —At x A 2 , 
B = -Rr 1 ^ in (25). 

4. Show that the hypothesis that A x is nonsingular in Ex. 3 is essential by proving 
that A\x — I and B& — I are equivalent in F[x\, yet P and Q do not exist, if 

/0 1 0\ /0 1 0\ 

A x = 0 1 0 , J5i = (0 0 l) . 

\0 0 0 / Vo 0 1 / 

6. Characteristic matrices with prescribed invariant factors. If g i(x), 
. . . , g t {x) are the nontrivial invariant factors of a matrix xI—A and n» is 
the degree of gi(x), then by Theorem 1 the n-rowed square matrix 

B = diag {B h . . . , B t ] 


will be similar in F to A if is an w»-rowed square matrix such that the 
only nontrivial invariant factor of the characteristic matrix of is g»(x). 
For then xl — B = diag \xl nx — B i, . . . , xl nt — B t ) is equivalent in 
F[x] to diag {Gh, . . . , G t \, where Gi = diag {g i} 1, . . . , 1}, and we con¬ 
clude that xl — B has the same invariant factors as xl — A . Thus the problem 
of constructing an n-rowed square matrix A whose characteristic matrix has 
prescribed invariant factors is completely solved by the result we state as 
Theorem 9. Let g(x) = x n — (bix n ~ l + . . . + b a ) and 


(32) 


1 0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 


ba—i 

bn—2 • • • 

bi , 


Then g(x)I is both the characteristic and the minimum function of A, g(x) is 
the only nontrivial invariant factor of xl — A. 

For 


(33) 


xl - A = 


1 z 

-1 

0 

. . 0 

0 ' 

0 

X 

-l 

0 

0 

0 

0 

0 

X 

-1 

~bn 

“bn-1 

bn— 2 • 

.. -b t 

x — bi | 
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The complementary submatrix of the element —b» of (33) is an (n — 1)- 
rowed square triangular matrix with diagonal elements all — 1 and its de¬ 
terminant is (—l) n_1 . Thus the cofactor of —b n is (—l) n+1 (—l) n_1 = 1 
and hence d*_i(x) = 1. It remains to prove that d„(x) = \xl — A\ — g(x). 
This is true if n = 1, since then A = A\ = (&i), and \xl — A \ = x — bi. 
Let it be true for matrices A n -i of the form (33) and of n — 1 rows, so that 
| xl — 4,_i | = x n ~ l — (&i£"“ J + . .. + b n - i) is the cofactor of the element 
x in the first row and column of (33). We now expand (33) according to its 
first column and obtain \xl — A \ = x[x n-1 — (6ix n- * + ... + i>„_i)] — 
b n = g(x) as desired. This proves our theorem. 

The construction of square matrices A with complex number elements 
whose characteristic matrices have prescribed elementary divisors has a 
simple solution, and we shall see that the argument preceding Theorem 9 
reduces the solution to the proof of 

Theorem 10. Let c be a complex number, A be the n-rowed square matrix 

c 1 0 ... 0 0 

0 c 1 ... 0 0 

(34) A = . 

0 0 0 ... c 1 

i 0 0 0 ... 0 c 

Then the only nontrivial invariant factor of x I — A is (x — c)“. 

For xl — A is a triangular matrix with diagonal elements all x — c, 
\xl — A\ = (x — c) n . The complementary minor of the element in the 
nth row and first column of xl — A is a triangular matrix with diagonal 
elements all unity, d n ~i(x) = 1, and d n (x) - (x — c)” is the only nontrivial 
invariant factor of A. 

Thus if Ci,... , c t are complex numbers and n\,... ,n t are positive in¬ 
tegers, we construct matrices Aj of the form (34) for c = c, and with n, 
rows. The matrix A = diag {Ai,. .. , A,) then has n rows, and its char¬ 
acteristic matrix xl — A = diag {B\,. .. , B t \, where B t — xl ni — Ai is 
equivalent in F[x] to diag {/<, 1,..., 1J such that /< = (x — c<) n< . But 
then xl — A is equivalent in F[x\ to diag \fi,... ,f h 1,... ,1], and by 
Theorem 3 the nontrivial elementary divisors of x/ — A are /i, ...,/< as 
desired. 
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EXERCISES 

1. Compute the invariant factors and elementary divisors of the characteristic 
matrices of the following matrices. 




d) 


1 -1 0 

2 -1 -1 

6 -2 — 2 , 


e) - 


4 6 3 

3 5 4 

4 —5 3 



o) 





'-6 2 -5 —19\ 

2 0 1 5 

-2 1 0-5 

3-1 2 9/ 


2. Find a matrix B = diag {B\, . . . , B t \ similar in F to each of the matrices of 
Ex. 1, respectively, where Bi has the form (32), and the characteristic function of 

is the /th nontrivial invariant factor of A. 

3. Solve Ex. 2 with the characteristic function of Bi now the ith nontrivial ele¬ 
mentary divisor of A , Bi of the form (34). 


7. Additional topics. There are many important topics of the theory of 
matrices other than those we have discussed, and we leave their exposition 
to more advanced texts. Let us mention some of these topics here, however. 
The quantities of the field K of all complex numbers have the form 

c = a + bi (a, b in R, i 2 = — 1) , 


where R is the field of all real numbers, a subfield of K. The complex con¬ 
jugate of c is 

c = a — bi , 


and the correspondence c <-> c defines a self-equivalence or automorphism 

of K. a fact verified in Section 6.9. This automorphism of K leaves the 
elements of its subfield R unaltered, that is, a = a for every real a. We 
then may call it an automorphism over R. 
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If A is any m by n matrix with elements in K , we define A to be the 
m by n matrix whose element in its ith row and jth column is a t y. It is then 
a simple matter to verify that 

AB = AB , (AY = A 7 , (F l = (U )~ l , 

for every A, /?, and nonsingular C. Moreover, then (AB)' = WA 7 . We 
now define a matrix A to be Hermitian if JP = A, skew-Hermilian if A 7 = 
— A. Two matrices A and B are said to be conjunctive in K if there exists a 
nonsingular matrix P with elements in K such that 

PAF = B . 


The results of this theory are almost exactly the same as those of the theory 
of congruent matrices, and it is, in fact, possible to obtain a general theory 
including both of the theories above as special cases. 

Two symmetric matrices A and B with elements in a field F are said to 
be orthogonally equivalent in F if there exists an orthogonal matrix P with 
elements in F such that PAP' = B. But PP' = P'P - I so that B = 
PAP “ 1 is both similar and congruent to A. Analogously, we call a matrix P 
with complex elements such that PP' = 7 17 P = /, a unitary matrix. Then 
we say that two Hermitian matrices A and B are unitary equivalent if 
PAP 1 = B , where P is a unitary matrix. Both of these concepts may also 
be shown to be special cases of a more general concept. 

Finally, let us mention the topic of the equivalence and congruence of 
pairs of matrices. Let A, /i, C, I) be matrices of the same numbers of rows 
and columns. Then we call the pairs A, B and C, B equivalent * pairs if 
there exist nonsingular square matrices P and Q such that simultaneously 
PAQ = C and PBQ = I). Similarly, if A, B ) C, I) are n-rowed square 
matrices, we call A, B and C, I) congruent pairs if simultaneously PAP' = 
C and PBP' = D for a nonsingular matrix P. 

References to treatments of the topics mentioned above as well as others 
will be found in the final bibliographical section of Chapter VI. We shall 
not state any of the results here. 

* In this connection see Exs. 3 and 4 of Section 5. 
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FUNDAMENTAL CONCEPTS 

1. Groups. These pages were written in order to bridge the gap between 
the intuitive function-theoretic study of algebra, as presented in the usual 
course on the theory of equations, and the abstract approach of the author’s 
Modern Higher Algebra. The-objective of our exposition has now been at¬ 
tained. For our study of matrices with constant elements led us naturally 
to introduce the concepts of field, linear space, correspondence, and equiva¬ 
lence, and we are ready now to begin the study of abstract algebra. How¬ 
ever, we believe it desirable to precede the serious study of material such 
as that of the first two chapters of the Modern Higher Algebra by a brief 
discussion of this subject matter, without proofs (or exercises). We shall 
give this discussion here and shall therewith not only leave our readers with 
an acquaintance with the basic concepts of algebraic theory but with a 
knowledge of how these concepts may lead into those branches of mathe¬ 
matics called the Theory of Numbers and the Theory of Algebraic Numbers. 

Our first new concept is that of a set G of elements closed with respect 
to a single operation, and we wish, to define the concept that G forms a 
group with respect to this operation. It should be clear that if we do not 
state the nature either of the elements of G or of the operation, it will not 
matter if we indicate the operation as multiplication. If we wish later to 
consider special sets of elements with specified operations we shall then re¬ 
place “product” in our definition by the operation desired. Thus we shall 
make the 

Definition. A set G of elements a, b, c, . . . is said to form a group with 
respect to multiplication if for every a, b, c of G 

I. The product ab is in G; 

II. The associative law a(bc) = (ab)c holds; 

III. There exist solutions x and y in G of the equations 

ax = b , ya = b . 

The reader is already familiar with the groups (with respect to ordinary 
multiplication) of all nonzero rational numbers, all nonzero real numbers, 
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and, indeed, all nonzero elements of any field. These are all examples of 
groups G such that for every a and b of G we have 

IV. The products ab = ba. 

Such groups are called commutative or abelian groups. An example of a 
nonabelian (noncommutative ) group is the group, with respect to matrix mul¬ 
tiplication, of all nonsingular n-rowed square matrices with elements in a 
field. 

Every group G contains a unique element e called its identity element , 
such that for every a of G 

(1) ae = ea = a . 

Moreover, every element a of G has a unique inverse a -1 in G such that 

(2) aa~ l = a~ l a = e . 

Then the solutions of the equations of Axiom III are the unique elements 

(3) x = a~ l b , y = ba ~ l . 

A set H of elements of a group G is called a subgroup of (7 if the product 
of any two elements of H is in H f II contains the identity element of G 
and the inverse element h~ l of every h of II. Then H forms a group with 
respect to the same operation as does G. 

The equivalence of two groups is defined as an instance of the general 
definition of equivalence which we gave in Section 4.5. The concept of 
equivalence of two mathematical systems of the same kind as well as the 
concept of subsystem (e.g., subgroup of a group, subfield of a field, linear 
subspace over F of a linear space over F ) are two concepts of evident funda¬ 
mental importance which are given in algebraic theory whenever any new 
mathematical system is defined. 

The number of elements in a group G is called its order. This number is 
either infinity, and we call G an infinite group; or it is a finite number n , and 
we call G a finite group of order n. For finite groups we have the important 
result which we shall state without proof.* 

Lemma 1. Let H be a subgroup of a finite group G. Then the order of II 
divides the order of G. 

A simple example of a finite abelian group is the set of all nth roots of 
unity. The reader may verify that an example of a finite nonabelian group 

* The proof is given in Chapter VI of my Modern Higher Algebra. 
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is given by the quaternion group of order 8 whose elements are the two- 
rowed matrices with complex elements, 



The set of all powers 

(5) a 0 = 6, a, a~ l ) a 2 , a -2 , . . . 

of an element a of a group G forms a subgroup of G which we shall desig¬ 
nate by 

(6) [a] 

and shall call the cyclic group generated by a. Its order is called the order of 
the element a and it can be shown that either all the powers (5) are distinct 
and a has infinite order, or a has finite order m, and [a] consists of the m 
distinct powers 

(7) e y a, a\ ... , a— 1 , 

where e is the identity element of [a], a m = e. Then the order m of a is the 
least integer t such that a* = e. Moreover, it can be shown that a 1 = e if 
and only if m divides t. 

The order m of an element of a finite group G divides the order of G f 
since m is the order of the subgroup [a], and we may apply Lemma 1. Thus 
n — mq , a n = ( a m ) q — e q — e. We therefore have 

Lemma 2. Let e be the identity element of a group G of order n. Then 

(8) a n = e 
for every a o/ G. 

2. Additive groups. In any field and in the set of all n-rowed square mat¬ 
rices there are two operations. Thus we have said that the set of all nonzero 
elements of a field and the set of all nonsingular n-rowed square matrices 
form multiplicative groups. But the elements of any field, the set of all 
m by n matrices with elements in a field, the elements of any linear space, 
all form additive groups, that is, groups with respect to the operation of 
addition as defined for each of these mathematical systems. The reader will 
observe that the axioms for an additive abelian group G are those axioms 
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for addition which we gave in Section 3.12 for a field. Additive groups are 
normally assumed to be abelian, that is, the use of addition to designate the 
operation with respect to which a group G is defined is usually taken to 
connote the fact that G is abelian. 

The identity element of an additive group is usually called its zero 
element, that is, the element 0 such that a + 0 = a = 0 + a. The in¬ 
verse with respect to addition of a is designated as —a and is such that 
a + ( — a) = ( — a) + a = 0. Thus the solutions of the additive formula¬ 
tion 

(9) a + x = b , y + a = b , 
of the equations of our group Axiom III are 

(10) x = (-a) + b , y = b + (-a ). 

When G is abelian, we have x — y and designate their common value by 
b — a. Thus we define the operation of subtraction in terms of that of addi¬ 
tion. 

In a cyclic additive group [a] the elements are always designated by 

(11) 0, a, -a, 2 • a, — (2 • a), . . . , 

where, clearly, if m is any positive integer — (w • a) = m • ( — a), and we 
define ( — m) • a = — (m • a). Here m • a does not mean the 'product of a 
by the positive integer m but means the sum a + . . . + a with m sum¬ 
mands. If [a] is a finite group of order m f the elements of [a] are 0, a, 2a, 
. . . , (ra — 1) • a, and m is least positive integer such that the sum of m 
summands all equal to a is zero. However, if [a] is infinite, then it may be 
seen that n • a = q • a for any integers n and q if and only if n and q are 
equal. 

3. Rings. The set consisting of all n-rowed square matrices with elements 
in a field F is an instance of certain type of mathematical system called a 
ring. Many other systems which are known to the reader are rings and we 
shall make the 

Definition. A ring is an additive abelian group of at least two distinct 
elements such ihat f for every a, b, c, of R, 

I. The product ab is in R; 

II. The associative law a(bc) = (ab)c holds; 

III. The distributive laws a(b + c) = ab + ac, (b + c)a = ba + ca hold. 
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We leave to the reader the explicit formulation of the definitions of sub¬ 
ring and equivalence of rings. They may be found in the first chapter of the 
Modern Higher Algebra. The reader should also verify that all nonmodular 
fields are rings, the set of all ordinary integers is a ring. 

The zero element of a ring R is its identity element with respect to addi¬ 
tion. Observe that by making the hypothesis that R contains at least two 
elements we exclude the mathematical system consisting of zero alone from 
the systems we have called rings. 

Rings may now be seen to be mathematical systems R whose elements 
have all the ordinary properties of numbers except that possibly the prod¬ 
ucts ah and ba might be different elements of R, the equations ax = b of 
ya — b might not have solutions in R if a ^ 0, b are in R. A ring may also 
contain divisors of zero , that is, elements a ^ 0, c 0 such that ac = 0. 
In particular, the ring of all n-rowed square matrices has already been seen 
to have such elements as well as the other properties just mentioned. 

A ring is said to possess a unity element e if e in R has the property 
ea = ae = a for every a of R. The element e then has the properties of the 
ordinary number 1 and is usually designated by that symbol. The unity 
element of the set of all n-rowed square matrices is the n-rowed identity 
matrix, and the unity element was always the number 1 in the other rings 
we have studied. However, the set of all two-rowed square matrices of the 
form 

< i2 > (2 o) 

with r rational may easily be seen to be a ring without a unity element. In 
fact, all nonzero elements of this ring are divisors of zero. 

A ring R is said to be commutative if ab = ba for every a and b of R. The 
ring of all integers is a commutative ring, the ring of all n-rowed square mat¬ 
rices with elements in a field F is a noncommutative ring. 

4. Abstract fields. If R is any ring, we shall designate by 

R* 

the set of all nonzero elements of R. Then we shall call R a division ring if 
R* is a multiplicative group. This occurs clearly if and only if the equations 
ax = b f ya = b have solutions in R * for every a and b of R *. The set of all 
n-rowed square matrices is not a division ring. However, let c and d range 



114 


INTRODUCTION TO ALGEBRAIC THEORIES 


over all complex numbers, 8 and i be the complex conjugate of c and d. 
Then the set Q of all two-rowed matrices 

< 13 > A - (S ~i) 

is a noncommutative division ring. The reader should verify this, noting 
in particular that every A 0 is nonsingular since A ^ 0 implies that 
c 7 * 0 or d 0 and | A | = cc + dd > 0. The ring Q is a linear space of 
order 4 over the field of all real numbers and it has the matrices I, A, B, AB 
of (4) as a basis over that field. It is usually called the ring of real quater¬ 
nions. 

Until the present we have restricted the term “field” to mean a field con¬ 
taining the field of all rational numbers. We now define fields in general. 

Definition. A field is a ring F such that F* is a multiplicative abelian 
group . 

The identity element of the multiplicative group F* is then the unity 
element 1 of F . The whole set F is an additive group with identity element 
0, and 1 generates an additive cyclic subgroup [1]. If this cyclic group has 
infinite order, it may be shown to be equivalent to the set of all ordinary 
integers. But F is closed with respect to rational operations, and [1] then 
generates a subfield of F equivalent to the field of all rational numbers. We 
call all such fields nonmodular fields. 

The group [1] might, however, be a finite group. Its elements are then 
the sums 

(14) 0 • 1 = 0, 1, 2, . . . , p - 1 , 

where p is the order of this group, and we have the property that the sum 
1 + 1 + . . . + 1 with p summands is zero. It is easy to show that p is a 
prime y and it follows that if a is in F , then the sum a + a + . . . + a with 
p summands is equal to the product (l + l + ... + l)a = 0. We call 
such fields F modular fields of characteristic p. It may easily be shown that 
the characteristic of all subfields of a field F is the same as that of F } and, 
in fact, every subfield of F contains the subfield generated by the unity ele¬ 
ment of F under rational operations. 

6. Integral domains. A commutative ring with a unity element and with¬ 
out divisors of zero is called an integral domain. Any field forms a somewhat 
trivial example of an integral domain. Less trivial examples are the set 
F[x] of all polynomials in x , the set F[x i, . . . , x q ] of all polynomials in 
X\, 9 x q and coefficients (in both cases) in a field F, the set of all ordinary 
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integers. The set of all integers may be extended to the field of all rational 
numbers by adjoining quotients a/6, 6^0. In a similar fashion we adjoin 
quotients a/6 for a and 6 ^ 0 in J to imbed any integral domain J in a field. 

When we studied polynomials in Chapter I we studied questions about 
divisibility and greatest common divisor. We may also study such ques¬ 
tions about arbitrary integral domains. Let us then formulate such a study. 

Let J be an integral domain and a and 6 be in J. Then we say that a is 
divisible by b (or that b divides a) if there exists an element c in J such that 
a = be. If u in J divides its unity element 1, then u is called a unit of J. 
Thus u is a unit of J if it has an inverse in J. The inverse of a unit is clearly 
also a unit. 

Two quantities 6 and 6 0 are said to be associated if each divides the other. 
Then 6 and 6 0 are associated if and only if 6 0 = bu for a unit u. Moreover, 
if 6 divides a so does every associate of 6. Every unit of J divides every a 
of J. Thus we are led to one of the most important problems about an inte¬ 
gral domain, that of determining its units. 

A quantity p of an integral domain J is called a prime or irreducible quan¬ 
tity of J if p 9 ^ 0 is not a unit of J and the only divisors of p in J are units and 
associates of p. Every associate of a prime is a prime. A composite quantity 
of J is an a ^ 0 which is neither a prime nor a unit of J. It is natural then 
to ask whether or not every composite a of J may be written in the form 

(15) a = pi . . . p r 

for a finite number of primes pi of J. We may also ask if it is true that 
whenever also a — qi ... q t for primes qi , then necessarily s = r and the q> 
are associates of the p, in some order. When these properties hold we may 
call J a unique factorization integral domain. The reader is familiar with the 
fact that the set of all ordinary integers is such an integral domain. This 
fact, as well as the corresponding property for the set of all polynomials in 
xij . . . , x n with coefficients in a field F are derived in Chapter II of the 
author’s Modern Higher Algebra. 

The problem of determining a g.c.d. (greatest common divisor) of two 
elements a and 6 of a unique factorization domain is solvable in terms of the 
factorization of a and 6. However, we saw that in the case of the set F[z] 
the g.c.d. may be found by a Euclidean process. Let us then formulate the 
problem regarding g.c.d.’s. We define a g.c.d. of two elements a and 6 not 
both zero of J to be a common divisor d of a and 6 such that every common 
divisor of a and 6 divides d. Then all g.c.d.’s of a and 6 are associates. We 
call a and b relatively prime if their only common divisors are units of J , 
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that is, they have the unity element of J as a g.c.d. We now see that one of 
the questions regarding an integral domain J is the question as to whether 
every a and 6 of J have a g.c.d. Moreover, we must ask whether there exist 
x and y in J such that 

d = ax + by 

is a common divisor and hence a g.c.d. of a and b; and finally whether or not 
d may be found by the use of a Euclidean Division Algorithm. 

6. Ideals and residue class rings. A subset M of a ring R is called an 
ideal of R if M contains g — ft, ag y ga , for every g and h of M and a of R. 
Then M either consists of zero alone and will be called the zero ideal of 72, 
or M may be seen to be a subring of R with the property that the products 
am , ma of any element m of M and any element a of R are in M. 

If H is any set of elements of a ring /?, we may designate by [H] the set 
consisting of all finite sums of elements of the form xmy for x and y in R, 
m in H. It is easy to show that {H} is an ideal. If H consists of finitely 
many elements mi, . . . , m t of R } we write {//J = {mi, . . . , m t \, and if H 
consists of only one element m of R, we write 

(16) M = {m\ 

for the corresponding ideal. This most important type of an ideal is called 
a principal ideal. It consists of all finite sums of elements of the form amb 
for a and b in R. When R is a commutative ring, M = [m] consists of all 
products am for a in R. 

The ring R itself is an ideal of R called the unit ideal. This term is de¬ 
rived from the fact that in the case where R has a unity quantity R = {1) . 
Evidently {OJ is the zero ideal. 

Let M be an ideal of R and define 

(17) a s b (M) 

(read: a congruent b modulo M) if a — b is in M. We may then define 
what we shall call a residue class a of M for every a of R. We put into the 
class every b in R such that a = b (M). Clearly a s b (M) if and only if 
b s a (.M ). Moreover, if a — b is in M and b — c in M, then (a — b) + 
(b — c) = a — c is in M. It follows that a = b (a and b are the same resi¬ 
due class) if and only if b is in g. 

Let us now define the sum and product of residue classes by 


(18) 


g + b == a b , 


g • b = a • b. 



FUNDAMENTAL CONCEPTS 


117 


It may be verified readily that if = g, 61 = b , then Qi + bi = a + b , 
flibi = It follows that our definitions of sum and product of residue 
classes are unique. Then it is easy to show that if M is not R the set of all 
the residue classes forms a ring with respect to the operations just defined. 
We call this ring the residue class or difference ring R — M (read: R minus 
M). When M = R the residue classes are all the zero class, and we have 
not called this set a ring. 

When the residue class ring R — M is an integral domain, we call the 
ideal M a prime ideal of R. We call M a divisorless ideal of R if R — M is 
a field. These concepts coincide in the case where R — M has only a finite 
number of elements since it may be shown that any integral domain with a 
finite number of elements is a field. This coincidence occurs in most of the 
topics of mathematics (in particular the Theory of Algebraic Numbers ) where 
ideals are studied. 

7. The ring of ordinary integers. The set of all ordinary integers is a ring 
which we shall designate henceforth by E. It is easily seen to be an inte¬ 
gral domain, and we shall prove that it has the property of unique factori¬ 
zation. 

We observe first that the units of E are those ordinary integers u such 
that uv = 1 for an integer v. But then 1 and — 1 are the only units of E. 
Thus the primes of E are the ordinary positive prime integers 2 , 3, 5, etc., 
and their associates — 2 , —3, —5, etc. Every integer a is associated with 
its absolute value \a\ ^ 0. We note now that if b is any integer not zero, 
the multiples 

(19) 0,161,-161,2161,-2161,... 

are clearly a set of integers one of which exceeds* any given integer a. 
Then let (g + 1)|6| be the least multiple of |6| which exceeds a so that 
g\b\ < a, (g + 1) | b | > a, a — g\b\ = r such that 0 < r < | b |. We put 
g = g if 6 > 0, q = —g otherwise, and have g\b\ = qb } a = bq + r. If also 
a = bq i + r x with 0 < r x < |6|, then b(q — qi) = r x — r is divisible by 6, 
whereas |r — r x \ < |6|. This is possible only if n = r and qi = q. We 
have thus proved the Division Algorithm for E , a result we state as 

Theorem 1. Let a and b tA 0 be integers. Then there exist unique integers 
q and r such that 0 <r< | 6 |,a = bq + r. 

We now leave for the reader the application of the Euclidean process, 
which we used to prove Theorem 1.5, to our present case. The process yields 

* We use the concept of magnitude of integers throughout our study of the ring E. 
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Theorem 2. Let f and g be nonzero integers . Then there exist integers a 

and b such that 

(20) d = af + bg 

is a positive divisor of both i and g. Then d is the unique positive g.c.d. of i 
and g. 

The result above implies 

Theorem 3. Let f, g, h be integers such that f divides gh and is prime to g. 
Then f divides h. 

For by Theorem 2 we have af +bg — 1 } afh + bgh = h . But by hypothesis 
gh = fq y h = {ah + bq)f is divisible by /. 

We then have 

Theorem 4. Let p be a prime divisor of gh. Then p divides g or h. 

For if p does not divide g, the g.c.d. of p and g is either 1 or an associate 
of p. The latter is impossible, and therefore p is prime to g. 

We also have 

Theorem 6. Let m be an integer . Then the set of integers prime to m is 
closed with respect to multiplication . 

For let a and b be prime to m and d be the g.c.d. of ab and m. If ab is not 
prime to m we have d > 1, and by Theorem 3 if d is prime to a it divides b. 
But then a divisor c > 1 of d divides a or b as well as m contrary to hy¬ 
pothesis. 

We may now conclude our proof of what is sometimes called the Funda¬ 
mental Theorem of Arithmetic. 

Theorem 6. Every composite integer a is expressible in the form 

(21) a = ±pi. . . p r 

uniquely apart from the order of the positive prime factors pi, . . . , p r . 

For if a = be, every divisor of b or of c is a divisor of a. If a is composite, 
it has divisors b such that 1 < b < |a|, and there exists a least divisor 
Pi > 1 of a. But then p\ is a positive prime, a = p\Oi for \a^\ < |a|. If 
02 is a prime, we write o 2 = ± p 2 with p 2 a positive prime and have ( 21 ) for 
r = 2 . Otherwise o 2 is composite and has a prime divisor p 2 by the proof 
above, o 2 = p^a* and o = pip 2 a 3 for |a 3 | < |a 2 |. After a finite number of 
stages the sequence of decreasing positive integers |o| > |a 2 | > |a 3 1 > 

. . . must terminate, and we have (21). If also a = ±<71 ... q B for positive 
primes gi, , 0 „ the sign is uniquely determined by a, and p\ . . . p r = 
0 i. . . 0 ,. Then either we may arrange the 0 < so that q\ = pi, or pi ^ 0 , 
for j = 1 , . . . , s. But if the divisor pi of 01 . . . 0 . is not equal to 01 , it docs 



FUNDAMENTAL CONCEPTS 


119 


not divide q\ and by Theorem 4 must divide q 2 . . . By our hypothesis 
it does not divide q 2 and must divide <73 . . . q a . This finite process leads to 
a contradiction. Thus Pi = < 71 , p 2 . . . p r = qi . . . q t , and the proof just 
completed may be repeated, and we may take p 2 — q 2 . Proceeding simi¬ 
larly, we ultimately obtain r = s, the pi = qi for an appropriate ordering* 
of the q it 

8. The ideals of the ring of integers. Let M be a nonzero ideal of the 
set E of all integers and m be the least positive integer in M. Then M con¬ 
tains every (dement qm of the principal ideal {mj. If A is in M, we may 
use Theorem 1 to write h = mq + r where 0 < r < m. But mq and h are 
in M y h — mq — r is in M. Our definition of m implies that r — 0 and thus 
that every element of M is in {mj. We have proved 

Theorem 7. The ideals of E are principal ideals {m}, m a positive integer. 
The residue classes of E modulo { m } are now the classes 

(22) 0, 1, — , m - 1 . 

For if a is any integer, we have a = mq + r for r = 0, 1, . . . , m — 1. 
Then a — r is in \m \, a = r. Thus the elements of the residue class ring 

E — {mj 

defined for m > 1 are given by (22). Then E — {m} is a ring whose zero 
element is the class 0 of all integers divisible by m and whose unity element 
is the class 1 of all integers whose remainder in Theorem 1 on division by 
m is 1. 

If a is an integer prime to m, the elements of the residue class q are all 
prime to m. For by Theorem 2 there exist integers c and d such that 

(23) ac + md = 1 . 

If b is in g, then b = a + mg, be + m(d — qc) = ac + m(qc + d — qc) = 
ac + md = 1, and therefore 6 and m are relatively prime. But a • c = 1, 
and Theorem 7 implies 

Theorem 8. The residue classes a in E — {m} defined for a prime to m 
form a multiplicative abelian group. 

If m is a composite integer cd where c > 1, d > 1, then m > c, > d, 
and c and d are both not the zero class. But c • d = cd = m = 0, E — {mj 
has divisors of zero and is not an integral domain. If m is a prime, then 

* We may order the p i so that p x <, p 2 <^ ... p r and similarly assume that qi ^ 
q 2 £ . . . £ q v Then we obtain r = s, p i — q { for i — 1, . . . , r. 
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every a not in m is prime to rn, and Theorem 8 states that 2 ? — {m\ is a 
field. We have proved 

Theorem 9. An ideal M of E is a prime ideal if and only if M = {pj for 
a positive prime p, M is a divisorless ideal. 

We observe now that a = b (Af) for M = \m\ means that a — b = mq, 
that is, a — 6 is divisible by m. Thus it is customary in the Theory of Num¬ 
bers to write 

(24) a = b (mod m) 

0 read : a congruent b modulo m) if a — b is divisible by m. But then a = 6, 
and if we also have c = d, we will have a + c = b + d as well as a • c = 
b • d. Hence if (24) and 

(25) c = d (mod m) 
hold, we have 

(26) a + c = 6 + d (mod m) , ac = bd (mod m) . 

Thus the rules (26) for combining congruences are equivalent to the defini¬ 
tions of addition and multiplication in E — {m}. 

We next state the number-theoretic consequence of Theorem 8 and Lem¬ 
ma 2 which is called Euler's Theorem and which we state as 

Theorem 10. Let f(m) be the number of positive integers not greater than 
m > 0 and prime to m. Then if a is prime to m we have 

(27) a f(m) = 1 (mod m) . 

For/(m) is clearly the order of the multiplicative group defined in Theo¬ 
rem 8. Our result then follows from Lemma 2. 

We next have the Fermat Theorem. 

Theorem 11. Let p be a prime. Then 

(28) a p = a (mod p) 
for every integer a. 

For (28) holds if a is divisible by p. Otherwise a is prime to p> and a is 
one of the residue classes 1,2,..., p — \ . Thus f(p) = p — 1 and by 
Theorem 10 a p ~ l — 1 is divisible by p } a(a p “ 1 — 1) = a v — a is also di¬ 
visible by p. 

The ring E — {p} defined by a positive prime p is a field* P whose 

* This field is equivalent to the sub field generated by its unity quantity of any field 
of characteristic p. 
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nonzero elements form an abelian group P* of order p — 1. The elements 
of P* are the distinct roots of the equation x p_1 = 1 and are not all roots 
of any equation of lower degree. Then it may be shown that P* is a cyclic 
multiplicative group [r] where r is an integer such that 1 , r, rf , . . . , are 
a complete set of nonzero residue classes modulo {p\ . Such an integer r is 
called a primitive root modulo p. W e then have 

Theorem 12. Let p be a prime of the form 4n + 1. Then there exists an 
integer t such that 

(29) t 2 + 1 ss 0 (mod p) . 

For p — 1 = 4n, and we let r be a primitive root modulo p, t = r\ Then 
t* = r p -\ P - 1 = (P + 1)(^ 2 - 1), (t 2 + !)(£? — 1) — 0 in the field E - 
{p\. However, t 2 = 1 since r is primitive, t? + 1 = 0 as desired. 

There is a number of other results on congruences which are corollaries 
of theorems on rings and fields. However, we shall not mention them here. 

9. Quadratic extensions of a field. If F is a subfield of a field K which 
is a linear space U\F + . . . u n F over F, we say that K is a field of degree 
n over F. The theory of linear spaces of Chapter IV implies that U\ may be 
taken to be any nonzero element of K. Hence we may take U\ to be the 
unity element 1 of F. Then if n = 1, the field K is F. We call K a quad¬ 
ratic, cubic, quartic, or quintic field over F according as n = 2, 3, 4, or 5. 

Let n = 2 so that K has a basis u\ = 1, u 2 over F. The quantities of K 
are then uniquely expressible in the form k = Ci + c 2 u 2 for Ci and c 2 in F, 
and k is in F, k = k • 1 + 0u 2 , if and only if c 2 = 0. Clearly, if k is in F, 
then 1, A; do not form a basis of K over F. We now say that a quantity u 
in K generaies Iv over F if 1, u are a basis of K over F. Then u generates K 
over F if and only if u is not in F. For if 1, u are linearly dependent in F we 
have a\ + a 2 u — 0 for a 2 5 * 0, u= —a^ l a\ is in F. 

The elements fc of a quadratic field have the property that 1, fc, k 2 are 
linearly dependent in F, c 0 fc 2 + Cik + c 2 = 0 for c 0 , c h c 2 not all zero and 
in F. If Co = 0, then c\ cannot be zero, and k = —C& 1 is in F, k is a root 
of the monic polynomial {x — k) 2 with coefficients in F. If k is not in F, 
then co ^ 0, fc is a root of a monic polynomial of degree two. Thus every 
element fc of a quadratic field is a root of an equation 

(30) /(x, fc) - x 2 - T(k)x + iV(fc) = 0 

with T(k) and N (fc) in F. In particular, if u generates K over F we have 

(31) u 2 — bu + c = 0 
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for b and c in F, and we propose to find the value of T(k) and N(k) as poly¬ 
nomials in fc, c and the coordinates ki and k 2 in F of k = ki + k 2 u . 

Define a correspondence S on K to K by 

(32) k = fci + k 2 u — > k s = fci + fc 2 ^ s , 

for all fci and k 2 in F where u s is in K. Then (32) is a one-to-one corre¬ 
spondence of K and itself if and only if u s generates K over F. If ko = 
k 2 + km for kz and fc 4 in F we have ftf = kz + fc 4 u 5 , fc + A; 0 = (fci + fc 3 ) + 
(fc 2 + fc 4 )w, so that (A: + k 0 ) s = fci + kz + (fc 2 + A; 4 )u* s , and we have 

(33) (Jfe + fc 0 ) 5 = + A: 0 5 . 


Also kko = A?iA: 3 + (fcifc 4 + k 2 kz)u + k 2 k^u 2 = (A;i/c 3 — A: 2 fc 4 c) + (&ifc 4 + fc 2 fc 3 
+ k 2 kj^i)u ) while k s k% = fcifc 3 + (fcifc 4 + k 2 kz)u s + k 2 k A (u s ) 2 . But then 

(34) (kk 0 ) s = 

if and only if (w s ) 2 = bu s — c, that is, u s is a root in K of the quadratic 
equation x 2 — bx + c = 0. But the quadratic equation can have only two 
distinct roots in a field K, the sum of the roots is b } 

(35) u s = u or b — u . 

In the former case S is the identity correspondence k <- > k. In either case 

S defines a self-equivalence of K leaving the elements of F unaltered and 
is called an automorphism over F of K. 

If K were any field of degree n over F and if S and T were automorphisms 
over F of K, we would define ST as the result k —> k ST = ( k s ) T of apply¬ 
ing first k-+k s and then k s —* ( k s ) T . It is easy to show that the set of all 
automorphisms over F of K is a group G with respect to the operation just 
defined. In case G has order equal to the degree n of K over F it is called 
the Galois group of K over F, and that branch of algebra called the Galois 
Theory is concerned with relative properties of K and G. In our present 
case u s * = ( u s ) s = (5 — u) s = b — (b — u) = u so that S 2 is the identity 
automorphism. But b — u = u if and only if 2u = 6, which is not possible 
since u is not in F unless K is a modular field in which u + u = 0. Hence 
if K is a nonmodular quadratic field over F, the automorphism group of K 
over F is the cyclic group [S] of order 2 and is the Galois group of K. 

We see now that if k = ki + k 2 u f then kk H = (ki + A: 2 u)[A;i + k 2 (b — 
u )] = k\ + k\kjb + k\u{b — u). But bu — u 2 = c. Hence if 

T(k) = k + k s , N(k) = kk s , 


(36) 
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we have 

(37) T(k) = 2ki + k 2 b , N(k) = k\ + , 

and the polynomial of (30) is 

(38) /(x, k) = (x — fc)(x — fc s ) . 

The function 7\fc) is called the froce of A:, and N(k) is called the norm of k . 
Moreover, we may show that 

(39) T(aik + a 2 A: 0 ) = a^fc) + a 2 T(k 0 ) , tf(Adb 0 ) = tf(fc) • AT(Jfc 0 ) 

for every a \and a 2 of F and k and fc 0 of K. For T{a\k + a* fc 0 ) = {ajc + a 2 k 0 ) + 
(ciik + a 2 k 0 ) s = ciik + a 2 k 0 + d\ k s + avkfi — a\{k + k s ) + a 2 (k 0 + k$) as 
desired. Similarly,* N(kk 0 ) = kk 0 (kk 0 ) s = kk 0 k s k$ = (kk s )(k 0 k$). Note that 
if A: is in F, then N(k ) = A; 2 . 

Let us now assume that the field K is nonmodular so that K = F + uF 
where u satisfies (31). Then K is also generated by the root u — \b of the 
equation 

(40) x 2 = a ( a in F) , 

if a = (u — b/2) 2 = u 2 — 6u + 6 2 /4 = — c + 6 2 /4. Let us then assume with¬ 
out loss of generality that u is a root of (40) so that in (31) we have 6 = 0, 
c = — a. Then (37) has the simplified form given by 

(41) T(k) = 2 ki , A (AO = k\- k\a . 

Now if a = d 2 for d in F, we have u 2 = d 2 , (u + d)(u — d) = 0, whereas u is 
not in F and u + d 0, u — d 0. This is impossible in a field. 

The quantities of K now consist of all polynomials in u with coefficients 
in F. For if k(x) is any polynomial in F[x], we may write k(x) = A;i(x 2 ) + 
k 2 (x 2 )x y k(u) = ki(a) + k 2 {a)u. Thus every nonmodular quadratic field is the 
ring F[u] of all polynomials with coefficients in F in an algebraic root u of an 
equation x 2 = a where a is in F, a is not the square root of any quantity of F, 
and F is nonmodular . Since K is a field, it is actually the field F(u) of all 
rational functions in u with coefficients in F. 

Conversely, if the nonmodular field F and the equation x 2 = a in F are 
given, then K is defined. For we may take 



* The trace function is thus called a linear function and the norm function a multipli¬ 
cative function. 
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and identify F with the set of all two-rowed scalar matrices with elements 
in F (so as to make K contain F). Then every polynomial in u has the form 
k = ki + k 2 u and N(k) — k\ — k\a = 0 if and only if ki = k 2 = k = 0. 
For otherwise k 2 0, a = (A^A^ 1 ) 2 contrary to hypothesis. But then we 
take K = F + uF and have Ar 1 = ( k\ — k\a)~ l (k x — k 2 u) for every k of 
K. Hence K is the field F(u). That K is nonmodular is evident since the 
unity element of K is that of F . We have thus constructed all quadratic 
fields K over a nonmodular field F in terms of equations x 2 = a for a in 
F, a 9 * <P for any d of F. 

Observe in closing that if K = F(u), then K = F(v ) for every v = bu 
such that b 0 is in F. But v is a root of 

(43) z 2 = b 2 a . 

Thus we may replace the defining quantity a in F by any multiple b 2 a for 
6 0 in F. It is also shown easily that if K is defined by a and K 0 by a 0 , 

then K and Kq are fields equivalent over F if and only if a 0 = b 2 a. 

10. Integers of quadratic fields. The Theory of Algebraic Numbers is 
concerned principally with the integral domain consisting of the elements 
called algebraic integers in a field K of degree n over the field of all rational 
numbers. We shall discuss somewhat briefly the case n = 2. 

Let, then, K be a quadratic field over the rational field so that K is 
generated by root u of u 2 = a where a is rational and not a rational square. 
By the use of (43) we may multiply a by an integer and hence take a inte¬ 
gral. Write a = <?d where d has no square factor and c and d are ordinary 
integers. If we take b = c~ l in (43) we replace a by d. Hence every quadrat¬ 
ic field is generated by a root u of the quadratic equation 

(44) x 2 = a = ±pi . . . p r , 

where the are distinct positive primes and r > 1 for a > 0, while if 
a = — 1 we interpret (44) as the case a negative and r = 0. 

The quantities k of K have the form k = k\ + k 2 u where ki and k 2 are 
ordinary rational numbers. We call k an integer of K if the coefficients 
T(k), N(Jc ) of (30) are integers. Thus k is an integer of K if and only if 
2 k x and k\ — k\a are both ordinary integers. We shall determine all integers 
of K stating our final results as 

Theorem 13. The integers of a quadratic field Kform an integral domain J 
containing the ring E of all ordinary integers . Then J consists of all linear 
combinations 


(45) 


Cl + C 2 W 
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for Ci and c 2 in E, where w = u if a = 2 (mod 4) or a = 3 (mod 4) bvA 


(46) 


w 


1 + U 
2 


if a = 1 (mod 4). 

Note that a ^ 0 (mod 4) since a has no square factor. We write 


ki = 


bi 

b Q > 



k = 


bi + b 2 u 
bo 


for b 0 , bi, b 2 in E and b 0 the positive least common denominator of k\ and k 2 . 
Then 1 is the g.c.d. of b 0 , bi , b 2 . Now 2/ci is an integer, b 0 divides 2bi, 
k\ — k\a = bo 2 (b\ — bfa), so that bg divides b\ — b|o. If p were an odd prime 
factor of bo, it would divide 2bi only if it divided b%. But then p 2 would 
divide 6?, and p 2 would divide bf — b\a as well as — bfa. Since a has no 
square factors this is possible only if p divides b 2 , a contradiction. Hence 
b 0 is a power of 2. If 6 0 > 4 divides 2bi, then 2 divides bi, 4 divides b\ and 
— bfa, and thus 2 divides b 2 . We have a contradiction and have proved that 
b 0 = 1, 2. If b 0 = 2 and bi is even, then 4 divides — bfa, and hence 2 
divides b 2 contrary to the definition of b 0 . Similarly, if b 2 were even, then 
4 would divide bf, a contradiction. Hence bi = 2mi + 1, b 2 = 2m 2 + 1 for 
ordinary integers m x and ra 2 , and b\ — b\a = 4 [m\ + m x — a(m£ + m 2 )] + 
1 — a is divisible by 4 if and only if a = 1 (mod 4). Thus we have proved 
that J consists of the elements (45) with w = u if a = 2, 3 (mod 4). But 
if a s 1 (mod 4), then we have shown that either k = bi + b 2 u with b t 
and b 2 in E, b 0 = 1, or k = mi + ??? 2 7 ' + w for ty in (46). However, u = 
2iy — 1, and in either case fc has the form (45) with Ci and c 2 in E. 

It remains to show that J is an integral domain. The elements of J are 
in the field K , and thus it suffices to show that k + h y kh are in J for every 
k and h of J. But the sum of k and h of the form (45) is clearly of that 
form, the product 

(ci + c 2 w)(di + d 2 w) = Cidi + (cid 2 + c 2 d x )w + c^ 2 w 2 

is of the form (45) if w 2 is of that form. But this condition holds since if 
w = u then w 2 — a + Ow, while otherwise u 2 = a = 4 m + 1 for an inte¬ 
ger m, T(w) = i + h = 1, N(w) = J(1 - a) = -m, u; 2 — w — m = 0, 


(47) 

This completes the proof. 


w 2 = m + io . 
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The units of J are elements k such that kh =' 1 for h also in J . Then 
N(kh) = N(k)N(h) = 1, N(k) is an ordinary integer which divides unity, 

(48) N(k) = ± 1 . 

Conversely, if (48) holds, k has the property kk s = 1 or — 1. But k s is in 
J when & is in J since T(k s ) = T(k) and N(k s ) = N(k). Hence k s = Ar 1 or 
— k s = Ar 1 . Thus (48) is a necessary and sufficient condition that A; be a 
unit of J . 

If a s= 2, 3 (mod 4) so that k = Ci + c 2 u, then (48) is equivalent to 

(49) c\ — c\a — ± 1 . 

However, if a = 1 (mod 4), then (48) becomes N(A;) = cj + C]C 2 — c\m = 
±1 so that 4 N(k) = 4 c\ + 4cic 2 + c\ — (4 m + 1 )c% = ±4. But this is 
equivalent to 

(50) (2 Cl + c 2 y - c\a = ±4. 

We may determine the units of J completely and simply in case a is 
negative. For both (49) and (50) have the form x\ + x\g = ±1, ±4 for 
g = — a > 0, and this is possible for ordinary integers Zi and x 2 and g > 4 
only if x 2 = 0, C\ = 1, — 1. Now a has no square factors, and hence <7^4, 
the only possible remaining cases are g = 1, 2, 3. If g = 2 we have xj + 
2xf = 1 only if x 2 = 0. Hence we have proved that the units of J are 1,-1 
for every a < 0 save only a = —1, —3. 

Now let a = — 1 so that (49) becomes c? + ci = 1. Then one of Ci and c 2 
is zero, the other is 1 or — 1, and the units of J are 

(51) 1, u, —u,u 2 = —1 . 

In the remaining case a = —3 we use (50) and have (2ci + c 2 ) 2 + 3 c\ = 
4, and if c 2 ^ 0 we must have c 2 = ± 1, 2ci + c 2 = ±1 with any choice 
of signs. Then c 2 = 1 gives C\ = 0 or —1, while c% = — 1 gives 0 or 1 
for Ci. Clearly, the units of J are 

(52) 1, —1, w, —w, w s , —w s . 

The units of J form a multiplicative group, and we have shown that if a is 
negative this group is a finite group. 

If a = 2, then h = 3 + 2u has norm 9 — 4 u 2 = 9 — 8=1. Hence h is 
a unit of J, and so is h* for every integer t. If h l = h B for $?* t, we may take 
t > s and have h l ~ 9 = 1. But we may regard u as the ordinary positive \^2 
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and have h > 5, h l ~ 8 > 5 £ ~* > 1. Hence the multiplicative group of the 
units of J is an infinite group. It may similarly be shown to be infinite for 
every 'positive a. 

We shall not study the units of quadratic fields further but shall pass on 
to some results on primes and prime ideals in two special cases. 

11. Gauss numbers. The complex numbers of the form x + yi with ra¬ 
tional x and y are called the Gauss complex numbers , and those for which 
x and y are integers are called the Gauss integers. Then the Gauss complex 
numbers are the elements of a quadratic field K of our study with 

u 2 = -1 , 

and the Gauss integers comprise its integral domain J. We have deter¬ 
mined the units ±1, ± u of J and shall now study its divisibility properties. 
Our first result is the Division Algorithm which we state as 

Theorem 14. Let i and g ^ 0 be in J. Then there exist elements h and r 
in J such that 

(53) f = gh + r 
and 0 < N(r) < N(g). 

For fg~ l is in K , fg~ l = fci + k z u with rational fci and fc 2 . Every rational 
number t lies in an interval s < t < s + 1 for an ordinary integer s. If 
s < t < s + \ then (t — s) < while if s + ?<t<s + 1 then \t — 
(s + 1) | < Hence there exist ordinary integers hi and h 2 such that 

(54) Si = ki — hi , s 2 = k 2 — h 2 , |*i| < 4 , |«a| < 4 • 

Put h = hi + h 2 u t s = si + s 2 u y r = sg so that fg~ l = h + s y f = gh + 
sg = gh + r. Then N(s) = s? + si < 4, N(r) = \(s)N(g) < N(g) as de¬ 
sired. 

We observe that the quotient h and the remainder r need not be unique, in 
our present case. For example, if / = 2 + u and g = 1 + w, we have 
N(g) = 2. Then / = l'g + l = 2- g — u with N( 1) = A r (—u) ~ 1. 

We shall use the Division Algorithm to prove the existence of a g.c.d. 
Our proof will be different and simpler than that we gave in the case of 
polynomials and indicated in the case of integers but has the defect of being 
merely existential and not constructive. We first prove 
Theorem 16. Every ideal M of J is a principal ideal . 

For the norms of the nonzero elements of M are positive integers, and 
there exists an m ^ 0 in M such that N(m) < N(f) for every/of Af. By 
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Theorem 14 we have / = mh + r for h and r in J and N(r) < N(m). But 
/, m and mh are in M, so is r = / — mh , and it follows that r must be zero. 
Thus/ = mh, M = [m\. 

The above result then implies 

Theorem 16. Let f and g be nonzero elements of J. Then there exist b and c 
in J such that 

(55) d = bf + eg 

is a common divisor of i and g. Thus d is a g.c.d. of f and g. 

For the set of all elements of the form xf + yg with x and y in J is an 
ideal M of J. By Theorem 15 we have M = {d J, d in ilf has the form (55) 
and divides / = 1 • / + 0<7, and g = 0 • / + 1 • g in M. 

The above result may now be seen to imply that if / divides gh and is 
prime to g, then / divides h, and also if p is a prime of J dividing gh, then p 
divides g or h. Moreover, if / is a composite integer of J, then / = gh for 
nonunit g and h, N(g) < N(f). Then if p is a divisor of / of least positive 
norm it is a prime divisor of f, and we continue the proof of Theorem 6 
to obtain 

Theorem 17. Every composite Gauss integer is expressible as a product 

(56) f = pi . . . p r 

of primes pi which are determined uniquely by f apart from their order and 
unit multipliers. 

We observe that Theorem 17 implies that an ideal M of J is a prime ideal 
(and, in fact, a divisorless ideal) if and only if M = {d| for a prime d of J. 
Let us then determine the prime quantities d = di + d 2 u of J . We see first 
that if d is a prime of J, so is d s — di — d 2 u. For if d s = gh with N(g) j* 1, 
N(h) 5 * 1, we have d * g s h s for N{g s ) = N(g), N(h s ) = N{h), and d is 
composite. We now prove 

Theorem 18. A positive prime po/E is either a prime of J or the product 
N(d) = dd s , where d is a prime of 5. Every prime d of J is either associated 
with a positive prime of E or arises from the factorization in J of a positive 
prime p = dd s of E which is composite in J. 

For if p is a positive prime of E and is composite in J, then p = dk for 
N(d) > 1 ,N(k) > 1 ,N(p) = p 2 = N(d)N(k). But then p = N(d). lid = 
gh with N(g) > 1, then JV(d) = p = N(g)N(h) so that N(h) = 1, h is a 
unit, d is a prime. Conversely, let d be a prime of J. Then c = N(d) is a 
positive integer of E and is either a prime or is composite. But c = dd* s can 
have at most two prime factors in E, c = pp 0 for positive primes p and p Q 
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of E. We may assume that p is associated with d and po with d s so that 
Po = Po is associated with d and with p. Since p and p 0 are both in E , we 
have p 0 = ± p. But p and p 0 are positive and must be equal, N(d) = p 2 . 
Therefore d is associated with the prime p of E which is prime in J. 

We now clearly complete our study of the primes of J by proving 

Theorem 19. A positive prime p of E is a prime of J if and only if p has 
the form 4m + 3. 

To prove this result we observe first that if t is an odd integer of E we 
have t = 2s + 1, t 2 = 4s 2 + 4s + 1 = 4 s(s + 1) + 1. One of s and $ + 1 
is even, t 2 = 8r + 1 1 (mod 8). If t is even, we have t 2 ss 0 (mod 4) 

and t 2 3= 0, 4 (mod 8). Thus a sum of two squares is congruent to 0, 1, 2, 4, 
or 5 modulo 8 while 4n + 3 is congruent 3 or 7 modulo 8. It follows that 
p = 4n + 3 t x 2 + y 2 . We now assume that p = gh for g and h in J and 
have p 2 = N(p) = N(g)N(h). If neither g nor h were a unit, we would 
have N(g) > 1, and both of these integers would be divisors of p 2 . But 
then N(g) = N(h ) = p = x 2 + y 2 , which is impossible. Hence, p = An + 3 
is prime in J. We note that 2=1 + 1 = (1 + u)(l — u) is composite in 
J and that it remains to show that p = 4n + 1 is composite in J. We know 
by Theorem 12 that there exists an integer h in E such that h 2 + 1 is divisible 
by p. If p divides 6 + u or h - u, then h ± u = p(fci + fc 2 w), and ± 1 = 
pfc 2 , which is impossible. But a prime p of J cannot divide the product 
(6 + u)(b — u) — b 2 + 1 without dividing one of its factors b + u> b — u. 
Hence p is not a prime of J. This completes the proof. 

We use the result above to derive an interesting theorem of the Theory of 
Nurtibers. We call a positive integer c a sum of two squares if c = x 2 + y 2 
for x and y in E. Then we have 

Theorem 20. Write c = f 2 g where f and g are positive integers and g has 
no square factors. Then c is a sum of two squares if and only if no prime factor 
of g has the form 4n + 3. 

For if c = x 2 + y 2 = (x + yu)(x — yu ), we may write x + yu = d\ 
... d r for primes d» in J, c — N{d\) . . . N(d r ). Then N(di) = p< is a prime 
of E if and only if p* it 4n + 3 and otherwise N(di) = pt Thus the prime 
factors of c of the form 4n + 3 occur to even powers and are not factors of g. 
Conversely, if g = pi . . . p r for positive primes p< of E not of the form 4n + 
3, we have p< = N(di) for d< in J, g = iV,(di) . . . N(d r ) = N(di . . . d r ), 
and c = N(fdi . . . d r ) = N{k ) = x 2 + y 2 for fc = x + yu in J. 

Note in closing that a positive prime p of the form 4n + 3 divides x 2 + y 2 
if and only if p divides both x and y. For p is a prime of J and divides 
( x + yu)(x — yu) if and only if p divides either x + yu or x — pu. Then 
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x ± yu = p{k\ ± k 2 u), x = pfci, y = pk 2 for integers fci and of #. It 
follows, then, that if x , y, z are ordinary integers such that 

(57) x 2 + y* = z 2 , 

every prime divisor p = 4n + 3 of z divides both x and y . Then let x, ?/, z, 
have g.c.d. unity. It follows that the odd prime divisors of z have the form 
4 n + 1. We may show readily also that x and y cannot both be even or be 
odd and that z must be odd.* 

12. An integral domain with nonprincipal ideals. We shall close our ex¬ 
position with a discussion of some properties of the ring J of integers of 
the field K defined by a = u 2 — —5. Since — 5 3 (mod 4), the elements 

of J have the form fci + k 2 u for fci and k 2 in the set E of all integers. Ob¬ 
serve that if b is an integer of E which is a divisor in J of fci + k 2 u , then b 
must divide both ki and k 2 . For k i + k 2 u = b(hi + h 2 u ), fci — bhi, k 2 = 
bh 2 . We now prove 

Theorem 22. The elements 3, 7, 1 + 2u, 1 — 2u are primes of J no two 
of which are associated. 

For if A; is a composite of J we have k = gh, N(k) = N(g)N(h) for ordi¬ 
nary integral proper divisors A r (^) and N(h) of N(k). The norms of the 
integers of our theorem are 9, 49, 21, 21, respectively, and the only positive 
proper divisors of these norms are 3, 7. But if g = g x + g 2 u } we have 
N[g) = g\ + > 0, g\ + 5 g\ = 3,7. Evidently g 2 ^ 0, g\ < 1. But g\ = 1 

is impossible since g\j± —2, 2. Thus 3, 7, 1 + 2u, 1 — 2u are primes of J. 
The units of J are 1,-1, and clearly no two of them are associated. 

We see now that 21 = 3 • 7 = (1 + 2u)(l — 2u), and we have factored 
21 into prime factors in J in two distinct ways. Moreover, the principal 
ideal {3} of J defined for the prime 3 is not a prime ideal. For 3 does not 
divide 1 + 2u and 1 — 2u in J yet does divide their product, and therefore 
the residue class ring J — {3J contains 1 + 2u and 1 — 2 u as divisors of 
zero. 

The ring J contains nonprincipal ideals one of which we shall exhibit. 
We let M be the ideal of J consisting of all elements of J of the form 

(58) 3x + 2/(1 + 2 u) (x, y in J) . 

If this ideal were a principal ideal {d}, there would exist a common divisor 
d of 3 and 1 + 2u. Since these are nonassociated primes, d must be a unit 

* For farther results see L. E. Dickson's Introduction to the Theory of Numbers , pp. 
40 - 42 . 
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1 or —1 of J, and {d j = {1). Then 1 = 36 + (1 + 2 u)c for 6 and c in J , 

7 = 216 + (1 + 2u)7c = (1 + 2u)[6(l — 2 u) + 7c], But 1 + 2u does not 
divide 7, a contradiction. We now proceed to prove that M is a divisorless 
ideal of J and hence is a prime ideal of J. 

We have shown that M is not a principal ideal and does not contain 1. 
Since M contains 3 but not 1 it cannot contain 2. For otherwise 3 — 2 = 1 
would be in M. Let c be an integer of E in M so that c = 3# + r where 
r = 0, 1, 2. Since M contains 3q it contains c — 3g = r, r = 0. Hence 
the ordinary integers in M are divisible by 3. 

The ideal M contains 3 and 1 + 2u and so contains 

(59) Wi = 3 , w 2 = 3u — (1 + 2 u) = u — 1. 

If h is in M, then 6, = 6 0 + h 2 u for h 0 and h 2 in E, h — h 2 w 2 = ho + h 2 in 
E and in M. By the proof above ho + h 2 = 36i for hi in 2?, 

(60) h = /wi + h 2 w 2 . 

We have thus proved that M consists of all quantities of the form (60) for 
hi and h 2 ordinary integers. Thus we may call w 2 a basis of M over E. 
Note that {1} has a basis 1, u over E in this sense. It may be shown that 
every ideal of J has a basis of two elements over E . 

Every integer of J has the form k = ki + k 2 u = ki + k 2 w 2 + k 2 . Write 
ki + k 2 — 3 c + r for c in E and r = 0, 1, 2. Then k = ewi + k 2 w 2 + r = 
r (Af), the elements of J — M are the residue classes 0,1, 2 such that 3 = 0. 
We have proved that J — M is equivalent to E — {3} and is a field, M is 
a divisorless ideal of J. 

This completes our discussion of ideals and of quadratic fields. We shall 
conclude our text with the following brief bibliographical summary. 

Let us begin with references to standard topics on matrices not covered 
in our introductory exposition. The theory of orthogonal and unitary equiv¬ 
alence of symmetric matrices is contained with generalizations in Chapter V 
of the Modern Higher Algebra and is further generalized and connected with 
the theory of equivalence of pairs of matrices in the author’s paper entitled 
“Symmetric and Alternate Matrices in an Arbitrary Field,” in Transactions 
of the American Mathematical Society , XLIII (1938), 386-436. See also pages 
74-76 and Chapter VI of L. E. Dickson’s Modern Algebraic Theories , and 
Chapter VI of J. M. H. Wedderburn’s Lectures on Matrices. Both of these 
texts as well as Chapters III, IV, and V of the Modern Higher Algebra in¬ 
clude, of course, all the material of our Chapters II-V. For a discussion 
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without details but with complete references of many other topics on ma¬ 
trices see C. C. MacDuffee, The Theory of Matrices (“Ergebnisse der Mathe- 
matik und ihrer Grenzgebiete,” Vol. II, No. 5 [pp. 110]). 

The theory of rings as given here is contained in the detailed discussion 
of Chapters I and II of the Modern Higher Algebra, and the theory of ideals 
in Chapter XI. See also the much more extensive treatment in Van der 
Waerden’s Modeme Algebra. The units of the ring of integers of a quadratic 
field are discussed on page 233 of R. Fricke’s Algebra, Volume III, and its 
ideals on pages 106-10 of H. Hecke’s Theorie der algebraischen Zahlen. We 
close with a reference to the only recent book in English on algebraic num¬ 
ber theory, H. Weyl’s Algebraic Theory of Numbers (Princeton, 1940). 
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J-rowed, 27 

Diag {Ai, .. . , Ak }, 31 
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Element —continued 
unity, 113 
zero, 112, 113 
Elementary: 
divisors, 94 
matrix, 92 
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cogredient, 52 
matrices of, 43 
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Fermat theorem, 120 
Field, 56, 114 
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independence, 67 
space, 66, 74 
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cogredient, 51 
inverse of, 17, 46 
matrix of, 36, 83 
nonsingular, 17 
product of, 36 
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Linear space, G6, 71, 74 
basis of, 67 
order of, 71 

Linear transformations, 84 
orthogonal, 86 

Mappings; see Linear mappings 

Matrices: 

addition of, 59 
commutative, 37 
congruent, 51 
congruent pairs of, 108 
conjunctive, 108 
equivalent, 47, 89, 92 
equivalent pairs of, 108 
orthogonal, 86 
orthogonally equivalent, 108 
products of, 36 

rationally equivalent, 25, 34, 47 
similar, 85, 103 
unitary equivalent, 108 

Matrix: 

adjoint of, 29 

augmented, 80 

of a bilinear form, 49 

characteristic, 101 

characteristic function of a, 101 

of coefficients, 19 

columns of, 20 

determinant of, 27 

diagonal, 29 

diagonal of, 23 

diagonal elements of, 23 

elementary, 92 

of an elementary transformation, 

elements of, 20 

equal, 21 

Hermit ian, 108 

identity, 30 

invariant factors of, 91 

inverse of, 45 

of a mapping, 36, 83 

minimum function of, 101 

nonsingular, 34 

partitioning of, 21 

principal diagonal of, 23 

rank of, 32 

rectangular, 19 


rows of, 20 
scalar, 30 
skew, 52 

skew-Hermitian, 108 
square, 20 

symmetric, 52, 53, 59, 64 
transpose of, 22 
triangular, 29 
unitary, 108 
zero, 30 

Minimum function, 101 
Minors, 27 

Monic polynomial, 3, 100 
Multiple roots, 5 

n-ary g-ic form, 13 
/^-dimensional space, 66 
Nonsingular: 

linear transformation, 84 
mapping, 84 
matrix, 34 
quadratic form, 55 
Norm in a field, 123 
Norm of a vector, 86 

Order of: 

a group, 110 
a group element, 111 
a linear space, 71 
Orthogonal: 
matrices, 86 
subspaces, 87 
transformations, 86 
vectors, 87 

43 Point, 66 

Polynomials: 
associated, 5 
coefficients of, 6 
constant, 1 
degree of, 2, 7, 97 
divisibility of, 5 
equal, 2 
factors of, 5 
g.e.d. of, 9 
homogeneous, 7 
identically equal, 0 
monic, 3, 100 
n-ic, 13 
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Polynomials —continued 
products of, 3 
9 -ary, 13 

rationally irreducible, 12 
relatively prime, 11 
roots of, 5 
scalar, 100 

in several variables, 6 
virtual degree of, 2, 6, 97 
zero, 1 

Positive form, 64 
Prime to a polynomial, 11 
Prime quantity, 115 
Primitive root, 121 
Principal diagonal, 23 

Quadratic field, 121 
integers of, 124 
Quadratic form, 54 
definite, 64 
index of, 62 
negative, 64 
nonsingular, 55 
real, 62 

Quadratic forms, equivalent, 59 
Quantities: 

associated, 115 
composite, 115 
congruent, 6, 11 
irreducible, 115 
prime, 115 
relatively prime, 115 
Quartic field, 121 
Quartic form, 13 
Quasi-field, 57 
Quaternary form, 13 
Quaternions, 114 
Quinary form, 13 
Quotient, 57 
right, left, 98 

Rank of: 

an adjoint matrix, 46 
a bilinear form, 49 
a matrix, 32 
a product, 44 
a quadratic form, 54 
a row space, 72 


Rational: 

equivalence, 25 
functions, 8 
operations, 8 
Real: 

quadratic form, 62 
quaternions, 114 
symmetric matrix, 64 
Reflection, 87 
Relatively prime: 
polynomials, 11 
quantities, 115 
Remainder: 
right, left, 98 
theorem, 4, 98 
Residue classes, 116 
Right division, 98 
Ring, 112 

commutative, 113 
difference, 117 
division, 113 
residue class, 117 
Roots, 5 
multiple, 5 
Rotation of axes, 87 
Row, 69 
rank, 72 
space, 69, 72 
Row by column rule, 37 

Scalar: 

matrix, 30 
polynomial, 100 
product, 16, 41 
Scalars, 19 
Self-equivalence, 107 
Semidefinite forms, matrices, 64 
Sequence: 

elements of, 16 
zero, 16 

Similar matrices, 85, 103 
Skew forms, 53 
Skew matrices, 52 
Space; see Linear space 
Spanning a space, 67 
Subfield, subgroup, 110 
Submatrix, 21 

complementary, 21 
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Subring, 113 
Subspaces, 66, 75, 77 
complementary, 77 
sum of, 77 
Subsystems, 110 
Subtraction, 112 
Symmetric bilinear forms, 54 
Symmetric matrices, 52, 53 
congruence of, 59 
definiteness of, 64 

Ternary forms, 13 
Trace, 123 
Transformations, 84 
Transpose, 22 
of a product, 37 
of a sum, 62 

Unary form, 13 

Unique factorization, 115, 118, 127 


Unit ideal, 116 
Units, 115 
Unity element, 113 

Variables, change of, 38 
Vectors, 66 

linearly independent, 67 
norm of, 86 
orthogonal, 87 
zero, 66 

Virtual degree, 2, 6, 97 
Virtual leading coefficient, 3, 97 

Zero: 

divisors of, 113 
element, 112, 113 
ideal, 116 
matrix, 30 
polynomial, 1 
sequence, 16 
space, 67 
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